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1. INTRODUCTION 


The work presented herein represents the results of the 
first phase of activities sponsored by the NASA Lewis Research Center 
under Contract No. NAS 3-22530. The activities were directed to the "Con- 
tinued Development of a Detailed Model of Arc Discharge Dynamics". This 
work is a continuation of previously documented results on the subject 
[Beers, et a1., 1979]. The purpose of these studies is to develop a model 
of the catastrophic breakdown and discharge of spacecraft dielectrics which 
have been charged by the geosynchronous substorm plasma environment. In pre- 
vious studies, a series of codes (SEMC, CASCAD, ACORN) were developed to 
describe the development of a single electron-initiated avalanche into a 
negative tip streamer. A conceptual model was also presented for a discharge 
mechanism which proceeds by a stochastic (random) succession of this basic 
process. The results presented herein are specifically directed toward a 
more fundamental understanding of this model. A major portion of the effort 
was directed towards characterizing the numerical solutions computed with 
the code models. 

A primary goal of the studies was to explore the parameter space 
associated with the numerical solutions of the code models. Preparatory to 
initiating these parametric studies, it was necessary to modify the code pack- 
ages in order to use them. The required modifications are described in Sec- 
tion 2 . Details of these changes and a user's manual are described in a sepa 
rate document entitled Operational Aspects of the Computer Codes: SEMC, 

CASCAD, and ACO RN , Beers Associates, Inc., Report No. 1-82-16-04. 

Section 3 of the report details the results of an extensive 
series of computations over a wide range of problen parameters. The re- 
sults are presented on a code-by-code basis. In Section 3.1 the results 
of SEMC canputations of the electron distribution function are presented. 
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Computational results from the CASCAD model of the undistorted avalanche 
are given in Section 3.2. A total of 60 separate calculations were per- 
formed for each of the above 2 sections. Section 3.3 describes the de- 
tails of ACORN calculations of the evolution of the self-consistent electron 
avalanche. Eleven distinct computational configurations were considered here. 

Section 4 gives a very brief discussion of the work on discharge 
initiation in space. The reader is referred to several appendices where 
published papers on the subject have been reproduced. The published results 
describe the activities performed during the contract period. Work of a more 
speculative nature is presented in Section 5 on discharge quenching and blow 
off mechanisms. The results of all work are briefly summarized in Section 6. 

Six appendices to the body of the text present reproductions of 
published papers which were generated either directly or indirectly by the 
work described in this report. Appendices 1 and 2 describe seni -analytic 
calculations of the electron distribution function. These results support 
the interpretation of the SEMC results of Section 3.1. ^ paper describing 

the preliminary results of Section 3.3 is reproduced as Appendix 3. Appen- 
dix 4 gives detailed calculations supporting an internal discharge initiation 
mechanism which was developed. In Appendix 5, calculations of thermal break- 
down in the presence of field distortion are given. Finally, further de- 
velopments in the classical approach to the discharge problem are given in 
Appendix 6. 
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2. CODE MODIFICATIONS 


An extensive set of modifications to the SEMC-CASC AD- ACORN family 
of computer codes previously developed has been implmented. The set of modi 
fications may conveniently be split into two areas: upgrades to the physical 
models; and installation of graphics. These modifications will be discussed 
in the following sections on a code-by*code basis. Operational features of 
the updated codes are separately documented in a user's guide report. The 
modified codes were delivered and installed on the Lewis Research Center 
UNIVAC 1108 system. 


2.1 SEMC MODIFICATIONS 


Three ma^Jor changes were made to the SEMC code: 

, • Modification of the mean free path sampling algorithm; 

t Computation of an equilibrium velocity- space distribution; and, 

• Automatic linkage to the CASCAD code. 

These modifications are discussed in turn below. 

The first change modified the mean free path sampling routine. 

In the usual Monte Carlo transport calculation, the energy of the particle 
does not change between collisions. For the required calculations however, 
Ue presence of large electric fields requires a new sampling approach. 
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The probability P(x) that a particle will have a collision along 
a path between x and x ^ dx Is: 


P(x)dx 




dx 


0 


where X Is the mean free path. In order to sample from this distribution, we 
equate the probability that a random number Is between C and C dC (denoted 
G(C)dC) to the above probability that a collision occurs between x and x * dx. 

By noting that G(c) * 1 Tor random numbers uniformly distributed on the interval 
[0.1], we obtain 



Or by the Integration 


( 2 . 1 . 1 ) 


/de -y i jexp f .y } d* • (2.1.2) 

0 0 0 

We observe however, that the Integrand on the right Is an exact differential 
and obtain: 



0 


tn(i - c) = - y* ^ . 

0 


(2.1.3) 


(2.1.4) 
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For those situations In which the particle is in force free flight 
between collisions, the total macroscopic cross section remains constant be- 
tween collisions. The mean free path, related to the total macroscopic cross 
section by: 


A 



(2.1.5) 


Is also a constant. Integration is easily performed to yield: 

X - - A tnd - 0 (2.1.6) 

Since 1 - C has the same distribution as K* we have the familiar result: 

X ■ - A i6n C . (2.1.7) 


In order to sample in time, one merely replaces x by t (time) and 
A by T (mean free time) in the above expressions. 

When there is an electric field present however, a charged particle 
will experience an acceleration between collisions, resulting In a constantly 
changing particle energy over the particle path. The mean free path (mean free 
time) is no longer constant in general (E^ is a function of energy) and the in- 
tegration of Equation (2.1.4) depends on the form of A. For a tabular form of 
A(E), the integration must be performed numerically. The schene presently im- 
plemented Is as follows: 

(1) Begin with a particle energy Eq. 

(2) Select a random number C*. 

(3) Numerically integrate 1/t from 0 to t' until the integration 
is equal (to a preprescribed accuracy) to tn(C'). 
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(4) The time to the next collision is then t' . 

(5) The coordinates of the next coMision site are given by an 
integration of the e^jations of motion. 

While the above discussion could easily be given in terms of dis- 
tance and mean free path, the time to the next collision site must be determined 
from a numerical inversion of a rather complicated expression previously de- 
veloped [Beers, et al., 1979]. An overall reduction in the complexity of the 
calculation is achieved by use of the mean free time formulation. 

The second major modification was the incorporation of an algorithm 
for scoring the equilibrium velocity distribution gQ(v). Scoring the equilibrium 
velocity distribution allows for the computation of the macroscopic transport co- 
efficients directly from appropriate velocity space integrations. Routines for 
performing these computations have also been included in the code. 

A method of scoring the equilibrium velocity distribution 9 q(v) has 
been devised. An understanding of this distribution may be obtained through the 
following argument. One may obtain 9 q(v) by visualizing a snapshot In time of 
the electron dynamics In phase space. The distribution 9 q{v) may then be ob- 
tained by tabulating the Instantaneous velocity of all particles In the snapshot. 
With the Monte Carlo methods however, we are able to follow only one sample 
particle at a time. One may view each trajectory as a contribution to the phase 
space snapshot by sampling the velocity of each trajectory at a single ran- 
domly selected time. By sampling a trajectory at many points In time, for ex- 
ample at a random time between each collision, the statistical uncertainty In the 
calculated value of 9 q(v) for a fixed number of particle trajectories will be 
greatly reduced. This is the method Implemented In the code. 

The transport coefficients are obtained by velocity space Integra- 
tions. The average velocity <v> Is simply 


-► 

<v> = 



( 2 . 1 . 8 ) 
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T^e drift velocity then* is the component of <v> in the direction of the field* 
and components perpendicular to the field should vanish with sufficiently good 
statistics. The mobility is obtained in the usual fashion. 

Now, introduce the random velocity u by the equation 

V » <v> + u (2.1.9) 

The mean random energy T (mean energy minus drift energy) is given by 

T ■ J* ^ m* (u*u)gQ(v)dv . (2.1.10) 

The energy may be split into coripi.nents parallel and perpendicular to the field 
in the obvious fashion. 


The diffusion coefficients are obtained from the integrals previously 
derived using the Chapman-Enskog theory [Beers, et al . , 1979]. Using the notation 
given in that reference, let K(v,Vq) be the scattering kernel for the Boltzmann 
equation. Let t be the following vector integral; 

T(vq) = y vK(v,VQ)dv . (2.1.11) 

It follows from symmetry arguments that t has the form 

T(vo) = Vi<’o> • 


where Kj is a scalar function of the variable (vq*Vq). The diffusion tensor B 
may be written as the following integral: 




dv 


9q(v) 

K,(v^) 


V u 


(2.1.13) 


7 



ORIGINAL 
OF POOR 

where u was Introduced in Equation (2.1.9). The tensor 6 may be diagonalized* 
with diagonal components being the coefficients along tie electric field and 
perpendicular to the electric field. 

Finally, the Ionization rote 3 Is given by 

3 • y* dv QqCv) , (2.1.14) 

where (v) Is the fundamental Ionization rate for an electron with velocity 
V. Routines for performing all the above Integrations have been Included In SEMC. 

The final major modification to SEMC was thn Installation of an auto- 
mated link between SEMC and CASCAD. The link was Implemented by Installing In 
SEMC a routine to write the distribution function of Ionization sites on disk. 

This file Is subsequently read by CASCAD with no user Intervention required. 


PAGE IS 

quality 
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2.2 CASCAD MODIFICATIONS 


In addition to the automated linkage described In Section 2.1, 
two other changes were made to CASCAD. The previous version of CASCAD was not 
able to follow the calculations for mc.e than 13 generations. This limit Is im- 
posed by the amount of central memory required to store the location of each 
electron In the cascade. This limit was extended beyond 24 generations by the 
use of a dual file flip-flop buffering scheme. At the end of each generation, 
the locations of the electrons are stored on one disk file. These locations 
are then read, updated and stored on the second disk file. At the end of each 
generation, the read/write pointers are Interchanged and the process repeated. 

The disk storage available was sufficient to run all cases of Interest. 

The final modification to CASCAD Involved the Installation of graph- 
ics display routines. At the end of each CASCAD run, the following four (4) plots 
are produced: 

• charge density contour plot; 

f electric potential contour plot; 

• electric field magnitude contour plot; and, 

a electric field vector plot. 
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2.3 ACORN MODIFICATIONS 


Modifications to ACORN were made primarily In the development of a 
complete graphics output system and for the replacement of the routine formerly 
used to print the arrays. A vector plot routine was developed to provide a 
clear concise display of the electric field. At many points within the calcu- 
lated grid an arrow Is drawn. The length of the arrow Is proportional to the 
magnitude of the electric field at that point. The direction Is that of the 
electric field vector at that point. The routine chooses the density of arrows 
to be as large as possible, without producing a cluttered plot. 

The contour routine developed for ACORN was designed to provide the 
option of either user-specified contour levels or contour levels chosen by the 
routine Itself. Past experience has Indicated that many automated graphics 
schemes frequently chose scaling levels which do not reasonably span the calcu- 
latlonal regimes of interest. For this reason a more flexible scheme seemed 
desirable. For the ACORN results presented here, the contour levels chosen by 
the routine have thus far been entirely satisfactory. The routines as Implemented 
in ACORN therefore use automated scaling as a default, with user override only 
when desired. 

The remaining plots, such as the electron number density versus 
one spatial variable with the other spatial variable fixed, were also generated 
with a more generalized plot routine. This routine can also be used on a stand 
alone basis for other purposes. Again the automated scaling within this routine 
has proven satisfactory for the case of Interest. This routine Is also used to 
produce time history plots of the problem variables. 

All three of the main plotting routines, CONTOUR, VECPLT and GRAPH 
were designed as general purpose stand alone routines. These routines are now 
available for Incorporation Into future codes. A canplete description of how 
to use these graphics packages Is Included In the appendix of the companion 
operational procedures report. 
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3. COMPUTATIONAL STUDIES 


An extensive series of computations was performed using the SEMC/ 
CASCAD/ACORN family of codes. The computations were performed on site at LEWIS 
RESEARCH CENTER. These computations are discussed on a code-by*code basis In 
the following text. 


3.1 SEMC STUDIES 


The computer code SEMC solves the linear Boltzmann transport equa- 
tion for the electron distribution function In the presence of an applied electric 
field. No electron self-scattering effects are Included In the simulation. The 
electrons undergo colllslonal events with a background set of scattering centers. 

A specification of the scattering events determines a model for the electron Inter- 
action with the background material. 

A crude model which was developed earlier [Beers, et al., 1979] 
represents the Interaction of conduction band (quasi -free) electrons In a 
solid material representing Teflon. It was the purpose of the SEMC compu- 
tational studies to explore the dependence of the electron distribution function 
(and resulting transport coefficients) on the parameters which define the scat- 
tering model. As all subsequent calculations In the CASCAD/ACORN series depend 
on the coefficients determined from the SEMC calculations, a determination of 
the sensitivity of these parameters on the material model also provides Infor- 
mation for subsequent CASCAD/ACORN calculations. 

The choice of scattering models for the SEMC computations, results 
the subsequent calculations, and a discussion of these results are presented 
below. 
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3.1.1 SCATTERING PARAMETER FORMULATION: SCALING 


As currently formulated, the single scattering Monte Carlo Code 
SEMC requires the definition of three distinct processes: elastic scattering. 

Inelastic energy loss (non-lonlzating) , and Ionization. The scale for these 
processes is set by some overall dimensional quantity having units of length 
(I.e., mean free path). Thus, the ^nput depends on three dimensional parameters 
^eV ^0* ^I paths for the respective processes. 

In addition, the Inelastic and Ionization processes are characterized by their 
energy thresholds (denoted fiu) and Ej respectively). These two dimensional 
parameters characterize the energy loss of the electron. The only other di- 
mensional parameters available are the electronic effective mass m*, the 
electronic charge e, and the electric field F. One might also Introduce the 
density of the solid N. Within the present formulation however. It will be 
assumed that N only enters the problem through the mean free paths (MFP), 
with all MFP's scaling Inversely with density. 

From this list of eight dimensional variables aq, Xj, Ej, 

'flu), m*, e and F, three independent variables may be chosen to represent the 
fundamental units of mass, length, and time. All other dimensional variables 
may then be expressed in terms of this basic set. A convenient choice is the 
electronic effective mass m*, the elastic mean free path X^^^, and the ionization 
energy Ej (a quantity conjugate to time). Denoting the fundamental dimensional 
units by a subscript "naught", we then have the following expressions for the 
units of mass, length, time, velocity, force and energy: 
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Mq ■ m* • 


■■0 ■ ■ 


/ m* 

■ 


*et ’ 






J 


3b 


and 


^0 = ^1 


( 3 . 1 . 1 ) 

( 3 . 1 . 2 ) 

( 3 . 1 . 3 ) 

( 3 . 1 . 4 ) 

( 3 . 1 . 5 ) 

( 3 . 1 . 6 ) 


All other dimensional quantities may be expressed in terms of the fundamental 
units of mass, length and time. For example, the units of the diffusion co- 
efficient Dq(LqTq^) are 



(It is implicitly assumed that the tanperature of the lattice is fixed in the 
above discussion, and that the tanperature dependence of the solution enters 
only through the parameters already given). 

The ranaining problem variables may be expressed in terms of the 
three basic variables, together with dimensionless parameters formed from 
various combinations of the problem variables. A convenient set of dimension- 
less parameters is: 
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E 

* iFX 


L 

el 


and 


fi 


FX 


e 

r" 

eil 


(3.1.8) 

(3.1.9) 

(3.1.10) 

(3.1.11) 

(3.1.12) 


where Aj and Aq are the length ratios and R, X, and are energy ratios. The 
quantity X is the ratio of ionization energy to energy gained in an elastic mean 
free path, while is the ratio of electrostatic Coulomb energy in a mean free 
path to the energy gain in a mean free path. The two quantities e and F enter 
the linear Boltzmann equation as a product e»F (in the external force term). 
Dimensional analysis of this equation gives rise to quantities which depend on 
this product, rather than on e and F individually. It follows that the equation, 
and hence the solution* is independent of the variable because n varies as 
the ratio of e and F. The variable enters only when the electron can act on 
itself (i.e., in collective problems). 


With the above information, it is an elementary task to rewrite 
all the problem variables in terms of the three dimensional variables, and 
the four dimensionless parameters Aj, Aq, R and X. Distinct problans (not 
obtainable by scaling) are therefore specified by the four dimensionless 
variables noted above. Of these, only the quantity X depends on the external 
field strength, so that X may be chosen as the variable specifying field 
strength. The remaining three parameters Aj, Aq, and R characterize the 
scattering parameter set. 
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For the purposes of this Investigation, we will now suppress 
one of these variables (Aj). Shown In Figure 3.1.1 Is a typical ionization 
cross-section (determining Xj). The major feature to note Is that the cross- 
section rises very sharply out of threshold, reaches a peak substantially 
above threshold, and then begins to trail off slowly for higher energies. 

For solid densities, cross-sections with this order of magnitude correspond 

o 

to rather short mean free paths (a few Angstroms). Except for enormous 
fields (greater than lO^^V/M), the energy gained In an Ionization mean free 
path Is completely negligible compared to the energy required to pass the 
peak In the cross- section. For moderate fields the electron will lose es- 
sentially all of Its energy In the collision Ionization. This Implies that 
the results will be very Insensitive to the exact value of the cross-sectiun 
above threshold. In conformity with all other Investigators, we will make 
the assumption that this parameter Is not of concern for the field strengths 
of Interest here; only If a study of runaway electrons were required would 
it make sense to investigate this parameter. For solids these runaways occur 
at field strengths substantially higher than those of Interest. 

The parameter variations for the present study are those In- 
volving R and Aq. For the model of Teflon chosen during the first phase of 
this work [Beers, et al., 1979], the optical phonon energy 'flu) was taken to be 
0.11 eV, while the band gap Ej was taken as 6.5 eV to a corresponding 
R-value of 0.017. Reasonable band gaps for insulators may be chosen In the 
range from 2 to 10 eV. Optical phonon energies in polymer type materials 
typically vary In the range of 0.04 to 0.4 eV. This range of energy values 
produces a spread of R-values from 0.004 to 0.2. We have chosen three R- 
values to span this range [0.2, 0.03, 0.004]; other values may be interpolated 
as required. The dependence of ionization on this parameter may be partially 
understood from the literature and is discussed In Section 3.1.3 below. 

In the Teflon Model previously Investigated the quantity Aq had 

a value of 4.4. Because cross-sections in polymers do not vary widely. It Is 
unlikely that this quantity will vary by more than a factor of 5 from this 
nominal value. We therefore choose the values of Aq to be [1, 5, 25]. 
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Figure 3.1.1 Total Ionization Cross Section in Ng ; 

From [Strickland, et al . , 1976]. 
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All combinations of the above parameters give rise to nine 
different cross-section sets. Because the code SEi<iC requires dimenslona. In- 
puts* the fundamental units must be given specific values* and the actual pa- 
rameters determining R and Aq must be given their appropriate values. Specif- 
ically the quantities A^^* Ej* and m* have the same values as In the previous 
study* whereas Hu) and Xq are varied to give the required variations In R and 
Aq. Also* the Ionization mean path retains the same value as previously. The 
fixed parameters are: 


Xej^ ■ 6 X 10“^®m * 

Xj = 2 X 10“^°m * (3.1.13) 

E » 6.5 eV 

I 


and last* the electronic mass (chosen as the bare electron mass). 

Table 3.1.1 defines the variations of Hw and Xq used to define 
a material model (cross-section set). Each set is given a numerical label 
from 1 to 9* spanning the range of R* Aq discussed above. In addition* a 
10th model is included in the table representing the nominal model developed 
in the previous study. [Beers* et al.* 1979]. 

The field values chosen for the computations were all at the 
high end of the physical scale. Based on general theoretical considerations 
(see Section 3.1.3 below)* the high field results may readily be extrapolated 
to lower field values. Six field values were chosen, spanning a single decade 
of field strength. These values were [4 x 10^* 2 x 10®* 8 x 10®* 6 x 10®* 

4 X 10®]* where the quantities In brackets are in volts per meter. These field 
valu s correspond to values of the dimensionless variable X of [2.7* 5.4, 10.8* 
13.5* 18.1. 27.1]. 
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MODEL 

NUMBER 

OPTICAL 

PHONON 

ENERGY 

(eV) 

R 

OPTICAL 

PHONON 

MFP 

0 

A 

^0 

1 

■■ 

m 

6 

1 

2 


■Bii 

6 

1 

3 


0.004 

6 

1 

4 

1.300 


30 

5 

5 

0.195 


30 

5 

6 

0.026 


30 

5 

7 

1.300 


150 

25 

8 

0.195 


150 

25 

9 

0.026 

0.004 

150 

25 

10 

0.110 

0.017 

26 

4.3 


I ^ Table 3.1.1 

Definition of Material Models by Parameters R, Aq 
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3.1.2 SEMC RESULTS 


The updated version of SEMC was exercised on the ten models 
of Table 3.1.1. Each material model was run for the six field levels d1s> 
cussed in Section 3 .I.I 1 for a total of 60 distinct computations. The re- 
sults are simmarized below. 


3. 1.2.1 MOBILITY 


Shown 4 n Figure 3.1.2 is the mobility for three different models 
(3, 6. 9). Each of these models has a fixed value of the optical phonon energy. 
The mean free path for optical phonon emission increases with Increasing model 
number. Thus, we see that an increase in optical phonon MFP leads to a de- 
crease in mobility, but that the absolute mobility is relatively insensitive 
to this parameter. Each of the other sets (for constant R) shows the same 
behavior. 


Shown in Figure 3.1.3 is the mobility for three models (4, 5, 6) 
which have a fixed value of Aq. In this case the optical phonon energy de- 
creases with Increasing model number. The dependence on optical phonon energy 
and field level is somewhat more complicated, as evidenced by the cross-over 
in the curves. In absolute terms however, the mobility is very insensitive to 
this parameter. Other sets with fixed Aq show similar behavior. 

It may be concluded that the mobility is very insensitive to 
the value of and only slightly sensitive to the value of Xq. The mobility 
does however, show a strong dependence on the value of the electric field, 
varying by a factor of 5 over a decade of field strength. 


19 





lfoDility< (M /S) 


ORIGINAL PAGL' IS 

OF POOR QUALITY 


10 


-3 




rm f 


— r" 

4.0 X 10 


9 


Electric Field (V/M) 


Figure 3.1.3 Mobility Versus Electric Field for Three 
Scattering Sets with Fixed A^. 
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3. 1.2. 2 DIFFUSION COEFFICIENTS 


Shown In Figure 3.1.4 Is the transverse diffusion coefficient for 
three scattering sets with fixed R. Increasing the optical phonon MFP corres- 
ponds to an increasing scattering set number. We also see that an Increase in the 
optical phonon MFP leads to an increase In the diffusion coefficient and that 
the absolute change is small. Other groups of scattering sets with fixed R 
show similar variation. Illustrated in Figure 3.1.5 Is the transverse diffusion 
coefficient for three scattering sets having fixed Aq. Increasing set numbers 
correspond to decreasing optical phonon energy. The dependence on this parameter 
(optical phonon energy) is much stronger, with variations of nearly 50% between 
sets 4 and 6. The relationship is also not monotonic as Is evidenced by the 
shuffled order and the apparent crossing at low fields. The other sets (for 
fixed Aq) show the same ordering with respect to 'hw. The field dependence of 
the transverse diffusion coefficient is quite weak. (The scale on the figure 
is linear and exaggerates the variation of this quantity). 

Shown in Figure 3.1.6, family (3, 6, 9) represents the longitud- 
inal diffusion coefficient for three scattering sets which have the same value 
of R. For small values of the fiold, an increase in Xq results in an Increase 
in diffusion. For high values of the field, this ordering is not present. A 
cross-over point is evident and the dependence on ''v, is quite weak. The cross- 
section family (2, 5, 8) shows similar behavior, but the family (1, 4, 7) is 
different in that there is a uniform increase in diffusion for increasing Xq. 

It appears that the cross-over point may be at much higher fields. All the 
families show a similar weak variation as a function of Xq. 

Shown in Figure 3.1.7 is the longitudinal diffusion coefficient 
for variations in Hw, with Xq fixed. Both the longitudinal diffusion coef- 
ficient and the transverse coefficient demonstrate an identical behavior with 
fio). Other families show the same behavior as well. 
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Figure 3.1.4 Transverse Diffusion Coefficient Versus 
Electric Field for Three Scattering Sets 
With Fixed R. 
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Figure 3.1.5 Transverse Diffusion Coefficient 
Versus Electric Field for Three 
Scattering Sets With Fixed A^. 
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Electric Field (lO^V/m) 


Figure 3.1.6 Longitudinal Diffusion Coefficient Versus 
Electric Field for Three Scattering Sets 
With Fixed R. 
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Electric Field (lO^V/m) 


Figure 3.1.7 Longitudinal Diffusion Coefficient Versus 
Electric Field for Three Scattering Sets 
With Fixed A^. 
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The variation In the longitudinal diffusion coefficient with 
field Is comparable to the variation of the transverse coefficient, varying 
by IcjS than a factor of 2 over a decade of field strength. 


3. 1.2. 3 IONIZATION LENGTH 


Shown In Figure 3.1.8 Is the Ionization length versus electric 
field for three scattering sets having a fixed value of flw. For low values 
of the field, an increase in the optical phonon MFP leads to a decrease In 
the Ionization length. For higher fields, this rather sharp difference de- 
creases and the curves coalesce, cross, and the dependence on Xq Inverts, 
but Is much weaker. Other families of scattering sets show a similar behavior. 

Shown in Figure 3.1.9 is the Ionization length for three scatter- 
ing sets with fixed Xq and different values of flui. An Increase in optical 
phonon energy at low fields gives rise to a strong increase in ionization 
length. At high fields there is a coalescence of the curves, showing a 
strong suppression of the dependence on this parameter. Other families for 
fixed Ag show similar behavior. 

The most remarkable feature is the very strong dependence on 
electric field, with variation of more than a factor of 200 over one decade of 
field strength. Because this is the quantity most sensitive to field strength, 
a coiipilation of all the computed ionization lengths is given in Figure 3.1.10. 
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Figure 3.1.8 Ionization Length Versus Electric Field 
for Three Scattering Sets With Fixed R. 
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3.1.3 DISCUSSION 


A momentum transfer collision rate Vp is related to the mobility 
p by the relation 


p 



(3.1.14) 


The rate Vp is related to the mean free path for momentum transfer Xp, and the 
mean speed v by the usual relation 


V 


D 



(3.1.15) 


Because the inelastic processes do not lead to any net momentum transfer on the 
average, the quantity Xp is equal to X^^^ (a constant). Thus, all the dependence 
of p on the model (and field) is through the mean speed v. As the field in- 
creases and the electrons heat up, with a consequent rise in v, the mobility de- 
creases as seen in the calculated results. Equation (3.1.15) is implicit, be- 
cause the mean speed includes a contribution due to the drift velocity Vp. Using 
the random velocity u (with magnitude u) defined in Equation (2.1.9), we have 
the obvious relation 





(3.1.16) 


Combining relations (3.1.14-3.1.16) together with the relation of drift velocity 
and mobility, the resulting equation may be solved for Vp. The following re- 
sult is found 
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V? ( j“ 

= (f(y)u) 


e F A 


el 


{■j m*u^) 


(3.1.17) 


(3.1.18) 


Equation (3.1.17) displays explicitly the dependence of the drift velocity on 
the mean random energy of the electrons. Combining Equation (3.1.17) with the 
definition of mobility gives 


u(F) . u(0) , 


(3.1.J9) 


where f(Y) is defined in Equation (3.1.17) and u(0) is given by 


u{0) = 


m*u 


(3.1.20) 


Equation (3.1.19) is a useful form for the mobility. The quantity p(0) has 
the dimensions of y, reduces to the zero field mobility at F = 0, and is a rela- 
tively weak function of the electric field. Most of the variation of the mobility 
due to the field is contained in the dimensionless function f(Y)/Y, which varies 
with the dimensionless quantity Y. The variable Y is directly proportional to 
the field F and inversely proportional to the random energy of an electron 

(1/2 m*u^). Shown in Figure 3.1.11 is the random energy of the electron for the 
'arious models computed above. While the energy varies strongly with electric 
field, it does not vary too strongly from model to model. 
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Figure 3.1.11 Random Energy Versus Field Strength 
for Various Scattering Models. 
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It is possible to feel comfortable with the relative insensitivity 
of mobility results to the material parameters together with the relatively 
stong dependence on field. Scattering is dominated by the elastic process. 

The drift velocity must saturate at high fields. For the ACORN computations 
of Section 3.3, we will assume a single fit to the mobility of all models, an 
excellent assumption based on the computed results. 

In elementary kinetic theory, the diffusion coefficient D is 
related to the average energy and the collision frequency v for randomizing 
collisions. 


0 


1 <v^> 

' 3 “' 


(3.1.21) 


o 2 

Equating <v‘"> with v and assuming that only the elastic collisions are 
randomizing, we find 


D 


3 '"et 


(3.1.22) 


The dependence of D on v given by Equation (3.1.22) shows that heating of the 
electrons will lead to an increase in D. This expression may be rewritten in 
terms of u, to yield 


D = 3 u (l + A + 4 y 2 j . (3.1.23) 

The explicit dependence of D c.n the field is contained in the quantity under 
the radical, while the dimensional dependence is contained in the factor out- 
side the radical . 
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Although the quantity Y varies relatively rapidly as a function 
of the field, the quantity under the radical is a very weak function of Y 
for the range of values taken on by Y. In a first approximation, the de- 
pendence on the field may therefore be suppressed. As displayed by Firgure 
3.1.11, D may be seen to vary as the square root of the random energy of 
the electron, p. This dependence accounts for the substantially slower varia- 
tion of D with both field and model parameters. It is apparent however, that 
increasing the energy of the electrons does increase diffusion. In the computa- 
tional approach to ACORN of Section 3.3 below, we will ignore both the model and 
field dependence of the diffusion coefficient. This seems well justified in view 
of the calculated results. 

As a final parenthetical note. Equations (3.1.23) and (3.1.19) 
may be combined to give an expression for the Einstein relation 

D 

p 

where 


The quantity ; is the distribution 
relation [Bates, et al., 1962]. 

The quantity displaying the strongest dependence on both the 
scattering model and field strength is the ionization length. The dependence 
of this quantity on both model parameters and field may be understood from the 
work of Lin, [1979] and Lin and Beers, [1981], which is based on the work of 
Baraff, [1962]. The first two of these papers are reproduced herein as 
Appendices 1 and 2 respectively because they are the result of previously 
sponsored NASA work on this subject. To facilitate comparison with results 


1 7 

z, • m*u‘^ 


(3.1.24) 


/l 4Y^%- 1 
\/l + 4Y^ - 1 


1/2 


(3.1.25) 


dependent coefficient in the general Einstein 
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In those papers, we present In Figure 3.1,12, a plot of the inverse ionization 
length a (in dimensionless units) versus the dimensionless field variable X of 
Equation (3.1.11). The plot is specifically chosen to be log-linear. It can 
be seen that most of the curves are very nearly straight lines on this plot. 
This provides a rationale for extrapolation tc lower field values corresponding 
to a fit to V of the form 


a « On^ exp 



(3.1.26) 


where oj^j and are fitting parameters. In fact, this form of a is used in 
the ACORN computations within Section 3.3 below. 
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Figure 3.1.12 Inverse Ionization Length a Versus the 

Dimensionless Field Variable x. 
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3.2 CASCAD STUDIES 


CASCAO Is a Monte Carlo computer code which uses the distribu- 
tion of Ionization sites computed by SEMC to generate an avalanche of elec- 
trons. New electrons are created at each step, and the resulting distribu- 
tion of charge Is used to compute the self-field. When the self- field reaches 
a pre-assigned value, the code scores the final distribution of positive and 
negative charges. This code was exercised for each of the solution sets gener 
a ted by SEMC, I.e., 10 scattering models at 6 different field levels. 

It Is possible to estlnete the fields due an avalanche by In- 
voking various approximations. A discussion of these estimates Is given In 
Section 3.2.1. These estimates are used as background for Interpreting the 
results of the CASCAD calculations. The computed results are discussed In 
Section 3.2.2. 


3.2.1 ANALYTICAL ESTIMATES 


In a continuum approximation to the basic avalanching process, 
the number density of electrons n_ is determined by the equation 


on 5 

-^ + V • (n_v) = 6n, + DVX * (3.2.1) 

If the self fields are ignored, then v, 3, and D are constants. The equation 
is then linear in n^. The quantity 3 provides a unit of time, while [v] 
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provides a unit of velocity. These together provide a unit of length \ (the 
ionization length) given by A = iv|/3i or B = a jv]. Equation (3.2.1) may be 
brought to dimeiisionless form by scaling out these quantities. Choose the drift 
velocity to be in the z direction. Equation (3.2.1) reduces to 


3n_ V 


• 9 

n_ + 


(3.2.2) 


where 



(3.2.3) 


is a dimensionless diffusion coefficient! and the space and time variables are 
now dimensionless. The solution to Equation (3.2.2) is 


n ( r' ,t ' ) 


o t '/ 
n_e I 


2 


) exp 






(3.2.4) 


The boundary conditions for the problem are such that the solution reduces to 
a 5 function centered at the origin at t' = 0, and vanishes at infinity. The 
solution represents a diffusion sphere which is translating in space with unit 
velocity. The r.m.s. sphere radius satisfies the relation. 


where <> indicates a density weighted average. Larger values of Cj therefore 
correspond to more rapid diffusion. 
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The total number of electrons N_ may be shown to Increase 

exponentially 


N° e^' 


(3.2.6) 


The density of positive charge sites n^ satisfies the equation 


9n^ 


= n 


(3.2.7) 


with n_ In the form given by Equation (3.2.4), the required Integral Is not 
analytically tractable, so that the resulting charge density and electric field 
cannot be brought to simple form. 

Reverting to dimensional variables, and working In the co-mov1ng 
frame (y = x - vt), the magnitude of the electric field F due to the electrons 
alone may be shown to be 

|F| = Fj fj |y37Tl^l/ro j . (3.2.8) 

where 

F . 3 e exp(Bt) ,325, 

e Is the electronic charge, Eq Is the permitivlty of free space, 

r^ = 6 Dt ; (3.2.10) 


I 
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fj(x) 


erf (x) 
.2 


2 e 

v/tt 


-X 


(3.2.11) 


In Equation (3.2.11), erf Is the standard error function [Abramowitz, 1965]. 

The function fj is plotted in Figure 3.2.1 • It may be seen to have a maxi- 
mum at x « 1 ( |y| = \/2/3 rg). The value at the max' mum is 0.428. The 
maximum value of the electric field is therefore 

= 0.428 F, . (3.2.12) 

max 1 


A convenient form of this equation may be obtained by using Equation (3.2.10), 
and introducing the dimensionless time t' = 6t 


F 


max 




where 


F 


1 


(0.428) 


eB 

ISTTEgD 


(3.2.13) 


(3.2.14) 


Introduce another dimensionless variable Cg equation 




2 


e a 


(3.2.15) 


where Fq is the ambient field. Equation (3.2,14) may then be written 


F 


1 


(0.428) 

16 tt 



(3.2.16) 
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Equations (3.2.16), (3.2.13), together with the spatial dependence given by 
fj(x) in Equation (3.2.11) give a complete characterization of the field due 
to the diffusing electrons. The field is radial and centered about the co- 
moving origin. 

The function (e /t') has a minimum at t' = 1. Assume the solution 
is only valid for t' > 1. Then this expression for the self-electric field gives 
a first estimate of when the self-field will become important. It will not 
become important until the field due to the negative charge alone becomes com- 
parable to the external field. When this field does become comparable to the 
external, however, it does not follow that the total self-field is important. 

The field due to the positive charge which is left behind will give a canceling 
contribution over much of the volume, leading to a reduction in field. 

The electric field due to the positive charge satisfies the 

equation 


3^+(x.t) ^ 

=-6^_(x,t) (3.2.17) 

by virtue of Equation (3.2.7). The total electric field f is therefore given 
by 

/ t 

df t (x,t') . (3.2.18) 

0 

with (x,t) given by Equation (3.2.8). Unfortunately, the integral in Equa- 
tion (3.2.18) is not analytically tractable (the situation is not helped by 
working with the potential). 

Estimates of the peak field due to Equation (3.2.18) are necessary. 

It is possible to compute the multipole moments of the charge distribution 
giving rise to the electric field (3.2.18). The dipole movement p is found 
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to be 


ev 


(1 - ' 


Since the total charge Q is simply e^^, we have the expression 


(3.2.19) 


p| = 2Qa 


where 2a is the length of the dipole and 


(3.2.20) 


2a = i|i- (1 - * (3.2.21) 

After a few generatioi this tends to a constant (the ionization length). Thus, 
a reasonable approximation to the fields is provided by the field associated with 
two diffusion spheres separated by an ionization length. 

A good estimate of this field is provided by the field from two 
spheres of uniform charge density separated by an ionization length. As may 
be seen from Equations (3.2.12) and (3.2.9), a maximum field matching may be 
accomplished by putting all the charge Q inside a sphere of radius 1.25 rg, 
which we denote r^. Three cases may be distinguished: 


(I) ^ 2r^ ; (3.2.22) 

(II) < X < 2r^ ; and (3.2.23) 

(III) ^ • (3.2.24) 
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Case (I) corresponds to two non-overlapping spheres, while case (III) cor- 
responds to two spheres with a strong overlap (separaclon of centers less than 
a radius). 


0 

With this model, a rigorous value may be derived for the maximum 
field strength. The relation Is 


F 

max 


F i /F (0»428) ^2 I ( * | 

o jya* 16 tt I ^,3/2 j 



(3.2.25) 


where the function h is defined by 


h(x) 


1 | x^-2x-t-2 I 
^ ix(x-l)2 j 


X > 2 


0.618 <x<2 (3.2.26) 


1 i (x+2) I 

2 |(l^x)2) 


0<x<0.618 . 


The function h(x) Is plotted in Figure 3.2.2. It has a maximum value of unity 

for small x (large time). The quantities in curly brackets will be recognized 

as the field of a dipole, written In dimensionless form. Note that this de- 

-3/2 

pends only on the dimensionless ratio (^ 2^1 quantity h gives a measure 

of the degree of overlap between the positive and negative charge. Note that 
for long times h tends to unity, at which time there Is substantial overlap of the 
distributions. Note that for small values of (diffusion less Important) h 
gives substantial corrections to the dipole field for long periods of time. 
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In the limit that a sufficient number of electrons are involved 
In an avalanche so that a continuum representation is sensible. Equation 
(3.2.25) should provide a rather good representation of peak electric field 
in an avalanche. The spatial distribution implied by the model should also 
be rather good over most of space. 

The rather complicated considerations were investigated in an 
attempt to understand the CASCAD computations presented in the next section. 
As will be explained in that section, the analytical treatment becomes in- 
adequate when only a small number of electrons are involved in the avalanche. 
Further discussion may be found in the next section. 
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3.2.2 CASCAO RESULTS AND DISCUSSION 


CASCAD cotnpul'.es the evolution of a single-electron initiated 
avalanche by a Monte Carlo method. The avalanche proceeds in the undistorted 
applied field. The self-fielo is computed at each stage of the process. When 
the self-fields reach a preassigned percentage of the applied field, the pro- 
cess is terminated. The purpose of this code is to provide the initial con- 
ditions for the ACORN computations of the next section. Unfortunately , it 
has not been possible to automate the interface of these two codes during the 
current period of activity. 

The lack of an automatic linkage between the two codes means that 
it is necessary for an analyst to look over the CASCAD output very carefully, 
and decide on a procedure for providing an analytical fit to the data. The 
data fit may then be used as an input routine for defining the initial condi- 
tion of the ACORN routine. 

Because of the large number of computations, it was decided to 
attempt a universal fitting procedure for the CASCAD output. The analytical 
considerations of the previous section were pursued in an attempt to provide a 
rationale for a fitting procedure. It was hoped that an adequate fitting of 
the data would be provided by using the charge densities given by Equations 
(3.2.4) and (3.2.7). 

The diagnostic which was used to compare the adequacy of a 
possible fit to the computed data was a comparison of the number of genera- 
tions required to reach a preassigned self-field level as obtained in 
CASCAD versus the same quantity predicted by Equation (3.2.25). Additionally, 
the functional form of the fields, as computed by CASCAD, was compared to 
the form of the fields as given by Equation (3.2.25) and related equations. 
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Qualitatively, the fields computed by CASCAD appear to be 
describable by the model of Equation (3.2.25) and related formulae. Quantita- 
tively. however, significant disagreement was found between the predicted 
values of the maximum field using Equation (3.2.25) and those computed from 
ACORN. The reason for this discrepancy has been uncovered. The understanding 
required to elucidate the discrepancy has shed new light on the avalanche 
process. 


The understanding and resulting conclusion Is as follows. When 
the self-fields of an avalanche become Important after a relatively few genera- 
tions 10), the transition to a continuum model of the process Is not strictly 
justified. A discrete description of the avalanche will have statistical features 
which are not present in the analytic continuum description of Section 3.1. The 
transition to a continuum description required by the CASCAD/ACORN Interface must 
therefore be viewed as a transition to an ensemble averaged description. Further 
definition of the Interpretation of the ACORN results will be required. Further 
discussion of the observations which lead to these conclusions is given below. 

For each of the material models and field strengths of Section 3.1, 
a CASCAD computation of the avalanche was performed (60 calculatiun:}. The 
calculation proceeded until the self-field reached 50% of the applied fu'*ld. In 
most cases, this required very few generations (approximately 10). AnalytJc 
computations of the number of generations required to reach this same field 
level were also performed. The first computations were made using the most 
rudimentary estimates of the field as presented In Section 3.3.2 of a previous 
report [Beers, et al . , 1979]. Preliminary spot checks comparing these analytical 
results with the CASCAD results suggested good agreenent. A later complete 
comparison of all 60 cases was performed. Significant differences appeared. 

Because of these differences. It was felt that perhaps the analytic 
estimates which led to the equation referred to In the previous paragraph were 
too crude. A more detailed estinate was therefore performed. The results of 
this estimation procedure are presented In Section 3.2.1 above. The best 
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estimate of the maximum field in the avalanche is provided by Equation 
(3.2.25). It was believed that this equation should correctly predict the 
maximum field within a factor of 2. The number of generations to reach a 
given field should therefore be predicted to within differences of one or two 
from the ac tua 1 . 

Table 3.2.1 provides a comparison of the number of generations 
required to reach 50 % of the impressed field for all sixty computations. The 
number of generations predicted by Equation (3.2.25) is shown in brackets, 
while the actual number achieved in the CASCAD calculation is unbracketed. 

It is evident from this table that good agreement between the two is not 
present. 


The discrepancies of Table 3.2.1 are similar when cruder estimates 
of the self-field are made. As noted above, the first approach to explaining 
the differences was to obtain a better analytic estimate of the field. As 
evidenced by Table 3.2.1, the better estimates did not eliminate the differences. 
The second approach looked more carefully at the parameters controlling the 
CASCAD calculations. 

The most obvious quantity to check is gridding which is used to 
perform the field calculations in CASCAD. Figure 3.2.3 summarizes the effect 
of grid size. The ordinate of the figure is the ratio of the number of genera- 
tions obtained by CASCAD to number predicted, as taken from Table 3.2.1. The 
abscissa is the ratio of the axial grid size which was used in the CASCAD 
calculation divided by the characteristic ionization length for the problem 
(as obtained from SEMC). The figure is constructed as a standard scatter-plot. 

It is obvious from the figure that the degree of agreement be- 
tween the CASCAD computed result, and the analytical prediction is very 
strongly related to the normalized grid size. Large grid size corresponds 
to analytic underprediction, while small grid size corresponds to analytical 
overprediction. A linear regression of the data is also shown on the figure 
(the dashed line). The equation of this fit is shown on the figure. 


an 


/ 


ORIGINAL PAGE t3 
OF POOR QUALITY 

TABLE 3.2.1 

Comparison of Number of Generations Required for Avalanche Self-Field 
to Reach 50:<. of Impressed Field CASCAD Result is Unbracketed. Pre- 
dicted Result from Equation (3.2.25) is Given in Brackets. 


Cross 

Section 

Set 

4 X 10^ 

Ambient Electric Field (V/M) 

2 X 10^ 1 X 10^ 8 X 10® 

6 X 10® 

4 X 10® 

1 

9(8) 

18(11) 





2 

5(8) 

6(9) 

11(9) 

18(9) 

19(11) 


3 

2(8) 

5(9) 

6(9) 

8(9) 

9(9) 


4 

5(10) 

9(11) 

16(12) 

18(13) 

19(13) 


5 

5(9) 

5(9) 

7(9) 

8(9) 

10( 10) 

14(10) 

6 

1(9) 

5(9) 

7(10) 

7(10) 

7(10) 

8(10) 

7 

5(10) 

7(11) 

9(12) 

11(12) 

14(13) 

16(13) 

8 

5(8) 

6(9) 

6(9) 

7(10) 

7(10) 

8(8) 

9 

5(9) 

7(10) 

6(10) 

7(10) 

7(10) 



6(10) 

7(10) 

8(11) 

10(12) 

12(12) 

17(12) 
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Figure 3.2.3. Scatter Plot Displaying the Covariance 

Between the Ratio of Computed to Predicted 
Generations to Achieve a Specified Relative 
Field Level and the Normalized Grid Size - 
see text for explanation. 
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A reasonable interpretation of Figure 3.2.3 is that the results 
of the CASCAD computations are strongly grid-size dependent. This hypothesis 
was tested by computing the identical problem for several choices of grid 
size. These comparative computations confirmed the hypothesis. Smaller 
grid size gives rise to larger values of the calculated fields. Larger grid 
size gives rise to smaller values of the computed field. The number of 
generations required to reach 50% of the applied field increases with increas- 
ing grid size. 

The fact that the CASCAD results are sensitive to the choice 
of grid size should not be construed to mean that the results are incorrect, 
but only that the results are subject to interpretation. Some comments must 
be made about the numerical algorithm. In principle, one can easily compute 
the electric field from a system of discrete charges by using Coulomb's law 
for point charges. It was deliberately decided to not follow this procedure 
in computing the self-field of the avalanche. The rationale for the decision 
is simple. The purpose of the calculation is to provide initial conditions 
for the transition to a continuum description of the avalanche. The quantity 
which is desired is the value of the self-electric field in this same con- 
tinuum description. Only if the field is large in this smeared description 
will it be important to the subsequent evolution of the smeared description. 

If a discrete representation of the field is used, it is very difficult to 
decide what is meant by the maximum field. Getting close to a bare change 
will always lead to a large field, independent of the number of particles in 
the avalanche. Not only is it desired to make the transition to a continuum 
approximation, but it is deliberately intended to exclude this field close to 
the bare charge. The force between bare positive and negative charges (in the 
short range sense) is meant to be handled via the primary ionization pro- 
cesses. 


We are led to the conclusion that a smeared description of the 
charge is necessary. Given this choice, the difficulties of gridding become 
apparent. If only a small number of generations are required to reach Important 
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levels of the self-field, then the statistical aspects of the avalanche may 
dominate. Choosing too small a grid leads to a significant probability that 
bins will end up with a single charge in them, leading to the same problems 
with large local fields as with Coulomb's law. Choosing too large a grid 
tends to give cancellation of the fields due to separated charge if the grid 
size is larger than the separation distance. The separation of charge occurs 
over an ionization distance. The conclusion that must be reached from this 
discussion is that the best one can do in gridding the problem is to choose 
the grid size to be very nearly equal to the ionization length. Choosing 
the grid either larger, or smaller tends to lead to poor results. This 
statement is borne out by the linear regression which is drawn in on Figure 
3.2.3. It passes very near the point (1,1). 

These observations are sufficiently important to restate. A 
given set of discrete charges from an avalanche has electric fields associated 
with it which can be computed exactly using Coulomb's law. It is desired 
to describe this set of charges in a continuum approximation. To do this 
requires a choice of distance scale for describing a "charge density". 
Equivalently, a choice of distance scale is required for averaging the electric 
field. The choice of distance scale must be consistent with the scale of the 
physics of interest. More fine-grained distance scales always lead to larger 
local fields, and larger fluctuations in the local fields. A choice of grid 
size in the CASCAD algorithm is a choice of the distance scale over which 
averages are to be performed. The computed fields will depend on the choice 
of grid size. 


The behavior of the peak computed local field in CASCAD and its 
dependence on grid size is entirely intelligible. The physical distance scale 
of importance is the avalanche length. Choice of a grid spacing significantly 
larger than this will lead to an averaging over positive and negative charges 
which will tend to obscure the effect which is being studied. Choosing grid 
size much smaller than this will lead to a very noisy (non smooth) local field 
description with larger values of the local field. The optimum choice of 
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distance scale is therefore very nearly equal to the ionization length. 

For this choice of gridding, the analytical results of Section 3.2.1 provide 
a reasonable representation of the CASCAD results. 

f' 

Of course, these results were implicit in the choice of algorithms 
in CASCAD. The very strong dependence on grid size shown in Figure 3.2.3 was 
not anticipated. This strong dependence is due to the relatively small number 
of particles involved in the avalanche. Inspection of the analytically pre- 
dicted number of generations in Table 3.2.1 shows it to be a very weak function 
of the problem parameters. The number of electrons is typically 2^^ = 1024. 

For this small number of electrons, the averaging will be highly sensitive to 
the choice of scale used for the averaging. The strong dependence on gridding 
for a relatively small variation in the grid size simply reflects the small 
number of particles on the avalanche. 

We conclude that the charge densities computed by CASCAD may be 
adequately represented with Equations (3.2.4) and (3.2.7). The parameter t' 
is taken to have the value given by Equation (3.2.6) where N_® is unity and 
N is the total number of electrons in the avalanche. These formulae have 
been used to provide the initial conditions for the ACORN computations of the 
next section. 


PAGE I? 
quality 


Because of the small number of particles in the avalanche when 
self-fields become important, it will be necessary in future investigations 
to determine more carefully the statistical fluctuations about the ensemble 
mean given by the above procedure. There may be significant features of the 
process contained in these fluctuations. 
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3.3 ACORN STUDIi 


The computer code ACORN is designed to compute the self-consistent 
evolution of the electron avalanche. A series of computations has been performed 
using this code for various values of the impressed field and several different 
material models. These computations are discussed below. 


3.3.1 ACORN RELATED ANALYSIS AND SCALING 


The computations to be performed were determined by the following 
considerations. The equations which are solved by the ACORN package depend on 
the three transport coefficients p» D, B (mobility, diffusion coefficient, ion- 
ization rate) and the impressed electric field Fq. The ionization rate B is re- 
lated to the inverse of the ionization length (denoted a) by the equation 
B = aVp where Vj^ is the drift velocity in the ambient field. The quantity a is 
often referred to as the first Townsend coefficient. Specifying a,p and the field 
provides a value of B. We will treat u, D and a as the set of independent transport 
coefficients. 


The transport coefficients are functions of the local field. From 
the results of Section 3.1, it was determined that the mobility and diffusion 
are very nearly independent of material model for the range of models and electric 
fields considered. The quantity displaying the greatest variation with material 
model (and field strength) is the inverse ionization length, a. This parameter 
was used as the discrimitiation variable. That is, it was decided to choose three 
material models which had first Townsend coefficients a which bracketed the 
apparent range of variation of the m's computed for the ten distinct models. 



Material model 10 (see Table 3.1.1) was specifically chosen as 
the nominal parameter set. Variations in material model parameters were chosen 
to bracket the parameters of this set. By referring to Figure 3.1.12, it may be 
seen that the ionization coefficient a for scattering set 10 is bracketed by 
ionization coefficients for the other scattering sets. The values of a on 
the high side (slow decrease with decreasing field) are clearly bracketed by 
those of scattering set number 9. This set was chosen as the set having the 
weakest dependence on field in the regime of interest. The decision on the 
low side was much more difficult. Scattering set number 1 had such a sharp 
dependence on field that values of a for low field were not obtained in the 
SEMC computations. It was decided to eliminate this scattering set on these 
grounds. The ionization is simply too low at many field strengths of interest. 
Scattering sets 2 and 4 have very similar slopes, but different absolute values. 
It was decided to use scatter*ng set number 2 to represent the lower bound on 
ionization behavior. Its behavior is most different in that it crosses some 
of the other ionization curves. 

The final choice of scattering parameter sets dictated that ACORN 
computations be performed using SEMC data from scattering sets 2, 9 and 10. 

Set 10 represents the nominal case; set 9 represents the high ionization limit, 
and set 2 represents the low ionization limit for the computations. The specific 
transport coefficients jsed for the computations were chosen as follows. The 
existing version of ACORN has the requirement that the diffusion coefficient be 
a constant, independent of location, field strength, or orientation. The SEMC 
computations presented in Section 3 above demonstrate that this constancy is a 
reasonable approximation. We have used the corresponding values from the SEMC 
calculations in each of the ACORN computations. These values are presented in 
Table 3.3.1 below. 


56 



ORIGINAL PAGE S3 
OF POOR QUALITY 


Electric Field 

Scattering Set V/m Diffusion Coefficient 


2 

2 

X 

10® 

1.18 

X 

10"^ 

2 

1 

X 

10® 

9.77 

X 

10"® 

2 

8 

X 

10» 

8.19 

X 

10‘® 

9 

4 

X 

10® 

1.42 

X 

►-» 

O 

1 

9 

1 

X 

10® 

1.28 

X 

^ ' 
1 

o 

9 

8 

X 

10® 

1.24 

X 

10"^ 

10 

4 

X 

10® 

1.73 

X 

o 

1 

10 

2 

X 

10® 

1.60 

X 

o 

1 

10 

1 

X 

10® 

1.59 

X 

1 

o 

10 

8 

X 

10® 

1.55 

X 

10"^ 

10 

6 

X 

10® 

1.49 

X 

10"^ 


Table 3.3.1 

Diffusion Coefficients from SEMC 
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The existing version of ACORN calls for a specific functional 
fit to the electronic mobility m as a function of the electric field. The 
results of the SEMC computations of Section 3.1 demonstrated that the dependence 
of the mobility on scattering parameter set was quite weak. A universal field 
dependent parameterization was therefore chosen for all ACORN calculations. The 
equation for the mobility p Is given by 

V » UM ' { 3 . 3 . 1 ) 


where 


=1.29 m^/v-s , 
A = 0.44 


(3.3.2) 


and the local electric field F Is expressed In Volts/meter. 

As evidenced by Figure 3.1.12, the dependence of the Ionization co- 
efficient a on local electric field F Is reasonably well represented by an ex- 
ponential. This type of parameterization has been chosen for all ACORN calcula- 
tions, with the specific parameters depending on the material model (scattering 
set). The functional form of the Ionization coefficient a was taken to be 

Fm 

a = Ojj(| exp [ - 1 ^] » (3.3.3) 

where a and have units of m’^ and the local electric field F has units of 
V/m. Table 3.3.2 gives the values of the parameters otj^ and F^ for the various 
scattering sets. 

The computations presented In Section 3.2 below were performed 
using the transport coefficients noted above. The values of the ambient field 
were chosen to span the regime of Interest, and are reported specifically in 
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Scattering Set 


(V/m) 

2 

1.37 X 10® 

3.55 x 10® 

9 

4.35 X 10® 

7.04 X 10® 

10 

5.10 X 10® 

1.63 X 10® 


a * exp [- -f ] 



Table 3.3.2 

Coefficients for the Equation of the Ionization Coefficient 
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the results section. As a result of these calculations, certain similarities 
in the computed solutions were noted. These similarities prompted an investi* 
gation into the anticipated scaling dependence of solutions of self-consistent 
avalanche equations. An investigation of scaling relationships sheds further 
light on the expected behavior of the solutions. A discussion of this feature 
of the problem follows. 

The equations which describe the evolution of the avalanche as 
solved by ACORN are; 


Ij- • 6 n. + ?(n.uf) + D?*n. 


(3.3.4) 

an. . n- 


(3.3.5) 

. ejnt -..nj , 

e 


(3.3.6) 

where n. is the mobile electron density, n^. the hole density, ? the electric 
field, u the electronic mobility, D the diffusion coefficient, e the electronic 
charge, c the dielectric yermitivity of the medium, and 3 = is the ionization 

rate, a being the inverse ionization length discussed above and Vq = - yF being the 
electronic drift velocity. Denote the ambient electric field by Fq, and the values 
of 3, and y evaluated at Fq as 3q, yg. Let the functions g^(F) and 92 (F) be 
defined by the equations 

u(F) ■ Uo 8 l<''> 


(3.3.7) 

B(F) = 92 (F) . 


(3.3.8) 

SO that gj and g 2 both have values of unity at F = Fq. Choose the unit of 
time Tq to be Sq~^, the unit of velocity Vq = MqFq, with the unit distance Xq 
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given by Xq * VqTq (ambient ionization length). If the unit of electric field 
is chosen to be the ambient field Fq, and all quantities are replaced with their 
non-dimensional analogs (which we denote with the " symbol) (e.g., t * ^Qt), 
Equations (3. 3.4 - 3.3.6) may be brought to the form 


8n- 

3t 


m 


■A /*«. A « A 

92". ^ * (9i"F) + n 




(3.3.9) 


3n+ 

— 

3t 



(3.3.10) 


V • F = ^2^"+ - "J ' . (3.3.11) 


where two dimensionless variables, and ^2 ^^^ve been introduced, and 
are given by the relations 


% 
ic 2 
^0 ^0 


^ _ o'^o 
^2 ' £ 


(3.3.12) 


(3.3.13) 


The variable Ci is the dimensionless diffusion coefficient which 

^ A 

was introduced in Section 3.2.1- The variable ^2 depends on the choice of 
the unit of density ti^. A free choice of Oq is possible. In particular, 
nQ may be chosen so that the variable has a value of unity. The choice 
leads to the following expression for ng: 
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"0 


^ 0^0 

eun 


e 


(3.3.14) 


The second expression in Equation (3.3.14) for Oq may be compared to an expression 
previously obtained for a "criticar charge density. If the critical charge density 
has the form en^, where n^ is a critical number density, then n^ is identical to 
the value of n^ given by Equation (3.3.14). Thus, choosing so that Cg ^^as a value 
of unity is equivalent to measuring number densities in terms of the critical 
number density n^. Recall that the critical charge density (and associated critical 
number density) occurs as a saturation charge density near the streamer tip. 

Equations (3.3.9 - 3.3.11) explictly display the scaling properties of 
the avalanche equation. Because the parameter may be chosen equal to unity, 
different problems are specified by different values of the parameter fTj,, and 
different forms of the functions g^ and g2» In the limit that g^ and g£ are 
weakly varying functions of the field (and hence nearly unity) different problems 
will be specified by the single parameter Cp This parameter is the same para- 
meter introduced in the discussion of the electron continuity equation in Section 
3.2.1. Different problems are specified by different values of Cj. Other 
parameter variations may be obtained by dimensional scaling. 

Of course, the functions g^ and gg introduce other parametric de- 
pendence in general. If the functional form for the mobility given by Equation 

A 

(3.3.1) is used, then the function gj depends only on the dimensionless field F, 
and introduces no further parametric dependence. Using the fits of Equations 
(3.3.1) and (3.3.3) gives rise to a form of the function gg which has the 
expression 

9o(F) = exp { ^ (1 - i )} F (1 - A) • <3.3.15 

^ '^0 F 

where A = 0.44 was derived in the above mentioned equations. 
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The function 92 depends on the dimensionless parameter (F|^/Fq). The range 
of variation and for the computed results Is noted below. Within 

the existing models used for ACORN, these two parameters are sufficient to 
describe all possible distinct solutions. 

Several other interesting features of the scaling should be noted. 

Let subscript "naught" again denote values of parameters in an ambient electric 
field Fq. Equation (3.3.12) gives an expression for a dimensionless constant 
which can be formed from some of these parameters. The set of parameters de- 
fining a problem is (Dq, Bq, Uq, Fq, e/t), together with the subsidiary relations 
Vq = UqEq, and Bq * This entire set of parameters may be characterized in 

a slightly different manner. One chooses three Independent dimensional parameters 
to provide the fundamental dimensions. The dimensions of any of the given 
quantities may then be written as powers of these basic dimensional units. Any 
given quantity may be rewritten as a dimensionless constant multiplied by the 
appropriate algebraic expression of fundamental units. 

For the above set of quantities, a convenient choice of the dimensional 
parameters is Bq, Vq = ^lo^O’ ^0* remaining parameters (Dq, e/t;) may be 
re-expressed as dimensionless constants times various powers of the dimensional 
units. With the indicated choice of units, the dimensionless diffusion coef- 
ficient is given by of Equation (3.3.12). The dimensionless coupling con- 
stant corresponding to e/t (denoted by ^ 3 ) is given by 
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The reason for the choice of notation, and the precise relation of ^2 
^2 Mill emerge below . 

The quantities and ^2 characterize the dimensionless parameters 
which can be formed from the given problem parameters. Shown In Figure 3.3.1 
Is a scatter plot of the values of and ^ 2 - constructed as 

follows. For each SEMC computation (10 materials models, 5 field levels) the 
values of and C 2 Mere computed. Each SEMC computation gives rise to a pair 

these pairs was plotted on the figure. 

From th<» plot, It may be seen that the two parameters show a 
very wide range of variation, (as much as five orders of magn1tud#» for ^ 2 )* 

The plot Is rather remarkable because a definite covariance Is shown between 
the values of these presumably Independent variables. A linear regression of 
the observed covariance yields the equation 


Cl = (0.379) • • (3.3.17) 

2 

This equation suggests that the quantity (^ 1 /^ 2 ) Is very nearly a constant over 
a wide range of parameter values. This quantity has the expression 


4 


C 


2 



(3.3.18) 


This equation, together with Equation (3.3.17) Implies that the transport 
coefficients Uq, Dq obtained from SEMC computations satisfy the rough equality 


°0 

s (0.379)- 

^0 



(3.3.19) 
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For the choice of models which have been computed. Equation (3.3.19) is an 
Einstein type relation which gives the field dependence of the mean electron 
energy. Though the relation (3.3.19) appears very reasonable, there is no 
a p riori explanation of this behavior. It remains a problem for future 
understanding. 

The net result of the above discussion is that a SEMC output is 
characterized by a single parameter, The second parameter '•s obtained 
from Cj by the relation (3.3.17). Two things should be noted. The first re- 
lates to the regression Equation (3.3.17) and the second relates to scaling of 
ACORN results. In a preliminary report of this work [Beers, et al . , 1981]; 
(reproduced here as Appendix 3) it was reported that nearly 

a constant for the parameters of interest. It was noted that the quantity 
(^2 /? 2 ) also occurs in Equation (3.2.25), the expression for the self-field 
of the undistorted avalanche. It was speculated that there was some deeper 
meaning associated with this coincidence. This speculation is now believed 
to be incorrect. 

The constancy of the quantity was obtained from an 

inadequate (eyeball) fit to the data. The properly weighted linear regression 
of the data represented by Equation (3.3.17) yields the approximate constancy 
of (Cj /C 2 ^- ® different result. Though the resulting expression 

for the Einstein relation is not understood a priori it does have a reasonable 
form. Were (C2^^^/^2^ constant, the resulting Einstein relation would 

contain an expotr* tial of the field on the right-hand side, a very difficult 
expression to reconcile. The constant value of (^ 2 /^ 2 ) ^Iso makes it clear 
that this relation is completely independent of the self-field effects of un- 
distorted avalanches. There is no relation. 

In the same publication [Beers, et al., 1981] it was further sug- 
gested that the covariance of the values of C 2 and exhibited by Figure 
3.3.1 provides an explanation of the similarity of the form of the solutions 
of the ACORN computations. While this statement is superficially correct, it 
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does not provide the proper interpretation of the phenomenon. As shown in 
the discussion on the scaling of the underlying equations, the physical de- 
scription may be cast in a dimensionless form which depends on the dimension- 
less parameters ^2* quantity ^2 given by Equation (3.3.13) de- 

pends on a choice of the unit density n^. If the unit density is chosen to 
be Oq = aQ*^, the quantity ^2 takes on the value of ^2* With this choice, the 
equations depend on the two parameters and ^2* 

However, as shown in the previous discussion of the scaling, the 

/V 

un't density may be chosen such that ^2 takes on a value of unity. Even if 
there were no covariance of and ^2* the problem may be brought to a form 
which does not depend on ^2* reason that a single parameter defines a 
problem is now apparent. There is no predefined natural unit of density for the 
problem. The only naturally occuring density is the critical density, a quant- 
ity which is dynamically determined. The conclusion to be drawn from these facts 
is that the dependence of the solution on the parameter Cp ''s trivial, and can 
be obtained by dimensional arguments alone. In the limit that g^ and gg are con- 
stant, problems solved by ACORN depend on a single parameter, i;^. 

The characterization of the scaling laws may now be completed. 
Dimensional analysis is used to bring the equations to the form (3. 3. 9-3. 3. II). 

A choice of Oq is made such that ^2 is unity. Explicit use is made of the ACORN 
functional fits for mobility and ionization rate to provide explicit functional 
forms for gj and g2- The two functions g^ and g2 depend on the dimensionless 
field F and a single dimensionless parameter (which we denote c^). The quantity 
is simply the ratio Fn^Fg of Equation (3.3.15). Thus, all problems solved by 
ACORN in this series of calculations depend on two parameters, Cj and ^4* There 
are no other parameters. 

The dependence on two parameters is only apparent. Shown in Figure 
3.3.2 is a scatter plot of the values of (Cj, c^) the plot was constructed by 
computing the values of and C4 which occur for the different SEMC computations. 

Each point represents a separate SEMC calculation. Again, a strong covariance 
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with very little scatter can be seen. The relation between c.j and dis- 
played by Figure 3.3.2 is easily understood. The quantity r.j depends on r.^ 
through the quantity Because the dependence is exponential the small 
variation in other problem parameters is marked. 

The conclusion is that for practical purposes, all ACORN 
computations may be shown to depend on a single paran»eter. The parameter 
(or may be chosen as the discriminating parameter. All other changes 
of problem parameters which do not vary the value of lead to solutions 
which can be obtained by dimensional scaling. 
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3.3.2 ACORN RESULTS AND DISCUSSION 


A serifiS of eleven computations was performed with the ACORN 

code. Transport coefficients were obtained for scattering sets (10, 9, 2) 

as described In Section 3.3.1. Scattering set 10 was used for ccmputatlons 

at 5 field values (4 x 10^, 2 x 10^, 8 x 10®, 6 x 10® volts/meter), 

while scattering sets 9 and 2 were used for computation at three field values 
g 9 8 

(4 X 10 , 2 X 10 , 8 X 10 volts/meter). The rationale vor the choice was 
discussed 1n Section 3.3.1. 

The plots which were generated as diagnostics of the calculations 
were quite remarkable. Visual Inspection of these plots leads to a very strong 
Impression that the form of the solution Is the same, Independent of problem 
parameters. That Is, the plots are all qualitatively the same as the plots 
which were previously presented for a single computation [Beers, et al., 1979]. 
This "sameness" of the solutions was regarded as highly significant. A sub- 
stantial Investigation was initiated to understand this occurrence. The results 
of that investigation will now be discussed. 

The investigation has been based on the notion of similarity 
solutions of mathematical models. The concept of similarity solutions of the 
hydrodynamical equations has a relatively well-defined meaning [Landau and 
Lifschitz, 1959]. In the simplest sense. It refers to solutions which depend 
on certain ratios of problem parameters. The most elementary examples in- 
volve one-dimensional time-dependent flows which are not characterized with 
either a time or distance scale. The solutions of the equations can be shown 
to depend on the ratio (x/t) only, and not x and t separately, where x is the 
spatial variable, and t the time variable. The use of the term similarity is 
related to the same word in geometry. 

After completing the work reported herein, an analytic similarity 
solution to the one-dimensional avalanche problem was discovered. It is a 
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similarity solution in the ordinary hydrodynamic sense. The details of 
this solution will be reported elsewhere. Its properties have not been 
fully Investigated as of this date. In the remainder of this section, a 
discussion will be given which presents some of the results which led to 
the belief that this solution existed. The importance of this type of 
solution is that all possible parameter variations can be accounted for by 
properly scaling and interpreting a single computed solution. 

The first step in a search for similarity solutions is to perform 
a detailed dimensional analysis. The results of that analysis were presented 
in Section 3.3.1 above. It was shown that the dependence of the solutions on 
the multiplicity of problem parameters could be represented by a single non- 
dimensional parameter, together with dimensional scaling. 

The second step in attempting to discover a universal similarity 
solution is to construct simple, easily viewed diagnostics of the solution. The 
plotted outputs of the quantities in two dimensions as a function of time are 
too complicated to provide simple measures of the similarity of the solutions 
(except in the pattern recognition sense that the visual appearance of the 
plots is the same). Simple quantitive measures are required. The remainder 
of this section provides some of the simple diagnostics of the solution which 
were investigated. 

The diagnostics presented below were investigated specifically 
to determine if the ACORN computation could be interpreted in terms of a 
single similarity solution. For the purposes of this investigation, a general- 
ized sense of the words "similarity solution" was felt to be important. In 
this context, we need a generalized idea of when two solutions are the "same". 
The most general point of view is topological. We refer the reader to the 
book by Thom [1975] as an example of this generalized point of view. The im- 
pact on the investigation reported below was to allow a wider latitude in 
the investigation. In all further discussions in this section the word 
similarity will have a deliberately vague meaning to allow for these possible 
general izations. 
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The result of the Investigations Is that all the ACORN computations 
appear to be Interpretable In terms of a single similarity solution. 

The simplest diagnostic of the evolution of the avalanche Is the 
total number of electrons Involved In the process. By Integrating Equation 
(3.3.4) over all space, the following equation Is obtained for the total 
number of electrons N 


dN 

M 


<3> N » 


where 

f d^x 3n_ 

< 3 > = 


(3.3.20) 


(3.3.21) 


is the density weighted average ionization parameter. The solution of 
Equation (3.3.20) is simply 


N = Nq exp 


r .t 

I j dt‘ <B(t' )> dt' • 


(3.3.22) 


For early times, when self-fields are unimportant <3> 3q where 3q Is the 
value of 3 in the ambient field Fq. It is sensible to plot the number N versus 
the dimensionless time 3Qt. Shown in Figure 3.3.3 Is a composite graph on 
which the computed results for all eleven computations are shown. From this 
plot (in dimensionless variables) it may be seen that the solutions have a 
general trend. As the self-fields build up, the rate of avalanching decreases, 
giving a decreasing slope to the graph. The number of electrons involved shows 
a remarkably small scatter. The mean fit shown on the plot represents the 
average behavior. It is apparent that the boundaries of the envelope of points 
have a very similar appearance. The composite plot was generated to display this 
similarity. Shown in Figure 3.3.4 is a plot of points from four separate runs, 
with a distinct symbol for each computation. In this figure It is easy to discern 
the distinct smooth curves of the individual outputs. The curves are clearly similar. 
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What 1s meant by this similarity may be made explicit by the 
following recourse to a simple example. From Equation (3.3.20) It Is ex- 
pected that N Is a unique function of time. Assume that the function Is 
Invertible and that the inverse Is analytic. Then <3>, which depends ex- 
plicitly on t, may be expressed as a function of the density N (chain rule). 
Assuming <B> to be an analytic function permits a Taylor series expansion 
about N = 0. The first term In the expansion Is clearly Sq. By retaining 
powers through the first order, It Is found that 



where Np Is dcifined by 


^0 _ d<6> 

Np = dN 


N ^ 0 


# 


The minus sign has been Introduced In Equation (3.3.23) to Indicate that the 
curve In Figure 3.3.3 is turning down. With the approximation of Equation 
(3.3.23) (which Is valid for sufficiently small N), Equation (3.3.20) becomes 


^ = ft N . !!- 
dt 


which win be recognized as the classical equation of population growth In the 
presence of competition [Davis, 1962]. The solution of Equation (3.3.25) Is 


N = 


1 + 


Nq / Sgt 

Nf r 





(3.3.23) 


(3.3.24) 


(3.3.25) 


(3.3.26) 
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Figure 3.3.3. Ccwnposite Scatter Plot Displaying 
Separate Computations of the Num- 
ber of Electrons Versus the Dimen- 
sionless Time 3gt. See text for 
explanation. 
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.3.4. Plot of the Nunter of Electrons Versus 
Bgt for Several Different ACORN Computa 
tions. See text for explanation* 
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where Nq is the population at t » 0. Figure 3.3.5 is a plot of Equation 
(3.3.26) for the condition that Nq » 1, and the three cases Np ■ 10^, 10^, 10^. 
The solution saturates at a value Np making ‘ihe notation transparent. These 
curves clearly have the same form and are only shifted with respect to one 
another. Shifting can be accomplished by plotting N/Nq versus BQt, and 
choosing Nq such that Nq/Np is a constant. 

Equation (3.3.25) is only an approximate expression, valid for 
small N. Higher order terms will very likely alter the behavior. In fact, 
there are reasons to believe that N increases linearly with t for long times 
[Beers, et al., 1978], so that the expansion of Equation (3.3.23) is at best 
asymptotic. The illustration serves to clarify the ideas of similarity and 
does have some physical content for small enough N. 


A second simple diagnostic also shows behavior suggestive of 
similarity. Let z denote the distance along the field measured from some 
reference position. Let the symbol < > denote density weighted average values 
of quantities. For instance, the mean location of the negative charge is 
given by 



(3.3.27) 


It follows from the equation of continuity for the electron density (Equation 
3.3.4) that the following equation holds: 


d <z> 

“Jt~ 


= <v^> + (<6z> - <3> <z>} 


(3.3.28) 


For early times, when self-fields are unimportant, <0z> ar <B> <z> because 
of the constant value of 6. Thus, for the undistorted avalanche, with 


<v 
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Figure 3.3.5. Plot of Equation (3.2.26 
Several Values of Np. 
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<z> « <Z>Q + (uQ^Qt) (3.3.29) 

or 

ttQ (<z> - <z>q) • $Qt • (3.3.30) 

Departures from a straight line on a plot of «q <z> versus 3Qt will depict 

the Integrated effect of the dynamics on the right-hand side of Equation (3.3.28). 

Shown 1n Figure 3.3.6 Is a composite plot constructed from all the 
computed solutions of the motion of the centroid of the electrons. Two branches 
of the solution are apparent on this graph, but most of the computed points are 
along a single curve. Shown on the figure Is a constant acceleration fit of the 
data (I.e.. quadratic In Bt). The value of the acceleration In dimensionless 
units BVq Is 2.0 x 10 . The most Interesting feature Is that the velocity Is 
Increasing, not decreasing. The enhancement of the field due to space charge Is 
generating a faster drift with the field, faster than drift In the ambient field 
alone. The possibility of this occurrence was suggested In a previous publication 
[Beers et al . , 1979]. 

Introduce the quantity 63 by the equation 


63 = ^ - <3> • (3.3.31) 

Assuming that the Initial value of <z> vanishes for t = 0, the solution of 
Equation (3.3.28) Is 


where 


<z> (t) = exp [B(t)3 


f df <v > (f) exp [- B(t')] , 

•'o 


B(f) 



63 (f) 


# 


(3.3.32) 


(3.3.33) 
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Figure 3.3.6. Location of the Centroid of Electrons 
as a Function of Time Plotted in Dimen- 
sionless Units. Composite of all ACORN 
Computations. Refer to text for explanation. 
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This form of the solution shows that the behavior seen in Figure 3.3.6 is 
rather remarkable. That is, the fit to the equation shows a quadratic de- 
pendence on time, while solution (3.3.32) shows explicitly an exponential. 

This demonstrates that either B(t) is very small, or that some other delicate 
cancellations occur. Of course, the solution (3.3.32) is implicit because 66 
depends on <z'^. This dependence may explain in general the weak time de- 
pendence displayed by the computed solutions. 

Other evidence suggests however, that the delicate cancellation 
is real. Referring to Equation (3.3.23), it may be seen that <6> depends on 
N. For short enough times, N depends exponentially on time as given in 
Equation (3.3.26). If this explicit time dependence for <B> is introduced 
into Equation (3.3.28), then the right-hand side contains a term with ex- 
ponential dependence, while the left-hand side is linear in time (from the 
computed solution). Thus, the other terms on the right-hand side of Equation 
(3.3.28) must also introduce exponential terms which provide a delicate cancel- 
lation to provide a linear behav^ir in t. 

The significance of this point is that it is suggestive of a 
similarity solution. Because <v^> is not expected to depend exponentially 
on time, the simplest solution to the cancellation problem is obtained by 
demanding that 66 be small, i.e., 

<6z'‘ - <z> 0. (3.3.34) 

If Equation (3.3.34) is approximately true, it implies a remarkable constraint 
on the solution. Recall that 6 is a very strong function of the electric field 
(and hence the charge density) . Relatively small variations in the field distri- 
bution could make large differences in the moments occurring in Equation (3.3.34). 

It appears that the negative feedback in the dynamical equations adjusts the charge 
and field distribution so that Equation (3.3.34) is approximately true. Because 
6 is such a strong function of E, it seems unlikely that there are many field 
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distributions which provide the minimization associated with this constraint. 

A natural interpretation is that the solution has moved to a stable equil- 
ibrium (in some sense) and continues in that state. 

Further diagnostics of the evolution may be obtained from the 
charge distribution. Inspection of contour plots of the distribution of 
charge gives a strong impression of similarity. Visually, there is little 
to distinguish the contours for the individual computations from those pre- 
viously presented [Beers, et al . , 1979]. An attempt to quantify this similar- 
ity may be made by considering the multipole moments of the distribution. 

These are scalar diagnostics of how the charge is distributed in space. ACORN 
computes these quantities as part of the numerical package. It must be men- 
tioned that they are computed in a moving frame. The procedure is as follows. 
The centroid of riegative charge is first computed (<z> in the above discussion). 
A similar computation is performed for the positive charge. The origin for 
the multipole expansion is then chosen to be at the midpoint between the two 
charge centroids. 

In Section 3.2.1, a discussion of the dipole moment p of the un- 
distorted avalanche was given. It was shown that p quickly tended to the limit 
p = 2Qa, where Q was the total charge of the electrons in the avalanche, and 2a 
was the avalanche length. An inspection of the dipole moment as computed by 
ACORN shows that this behavior persists throughout the self-consistent phase of 
the evolution. A plot of the dipole moment in time shows an identical behavior 
to the behavior of the total number of electrons involved in the process (Figure 
3.3.4). The dipole length shows a very small variation. 

The behavior of the quadrupole moment has also been investigated. 
Utilizing the techniques of Section 3.2.1, it is possible to derive an exact 
expression for the quadrupole moment of the undistorted avalanche. Denoting 
this moment as q, it is found that 


q 




+ SQt - 1 1 


(3.3.35) 
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where Q Is the total electronic charge, ag Is the Inverse Ionization length, 
and Bg Is the Ionization rate. After a few generations, the term in brackets 
increases linearly with time. An Interesting quantity to plot Is the dimension- 
less ratio (qag /4Q) as a function of the dimensionless time Bgt. 

Shown in Figure 3.3.7 is a plot of .x s dimensionless ratio as 
a function of Bgt. The origin of time for each individual computational set 
was chosen to correspond to the origin for the number of electrons In the distrl 
bution. At first sight it may seem that there Is little order In this figure. 

It must be remembered, however, that the quantity which is plotted is the ratio 
of two exponentially growing quantities. Small absolute errors in either of 
the quantities can lead to relatively large percentage errors In the ratio. A 
reasonable Interpretation of this plot Is that the quadrupole moments of the 
charge distributions are also very similar for different parameter choices. 
Generally they have the same behavior as given by Equation (3.3,35). It 
should also be mentioned that each Individual set of points for a given compu- 
tation varies smoothly with time. 

No higher order moments were computated in the ACORN calculations, 
as they were not required for the convergence of the numerical algorithm. 

All of the diagnostics discussed above support the statement 
made at the beginning of this section which noted that a visual inspection 
of the plotted output strongly suggested a great similarity of all solutions, 

A major portion of this similarUy was deduced in Section 3.3.1 where the 
scaling properties of the solutions were discussed. It was determined that 
the proMem could be reduced to solving a set of dimensionless equations 
depending on a single dimensionless parameter. The results of this section 
have shown that simple diagnostics go further than this and strongly imply 
that solutions having different values of the dimensionless problem parameter 
are also similar. It would appear that when properly Interpreted, there is 
but a single problem described by ACORN. 
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Figure 3.3.7. Dimensionless Quadrupole Moment qaQ^/4Q 
Versus Dimensionless Time BqL Composite 
for all ACORN Computations. 
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4. INITIATION MECHANISMS 


Because of mounting evidence that the majority of spacecraft 
dielectric discharges in space occur for low differential voltages [Stevens, 
1980], it was decided to expend the majority of effort for this task in- 
vestigating the effects of internal charge distribution. The results of 
this investigation were published. This paper, [Beers et al., 1981] reproduced 
here as Appendix 4, adequately discusses the effort and may be referred to for 
details. It is believed the results demonstrate that the internal charge dis- 
tribution has an important bearing on the occurrence of discharge. 

The behavior of avalanching and streamer formation as a function 
of the manifold material parameters follows from the computed results of section 
3. Local field enhancement increases the probability of occurrence, as does a 
decrease in scattering probability. The quantitative change is provided by the 
computed results. No quantification of the decrease in scattering probability 
near various defects has been given. Field enhancement due to the internal 
charge distribution has been given in Appendix 4. Field enhancement near de- 
fects has not been quantified. 

During the course of the research, it was noted that the classical 
literature on thermal breakdown failed to account for field distortion in the 
medium. A proper formulation of the problem was developed to account for the 
effect. A brief account of this work was published [O' Dwyer and Beers, 1980] 
and is reproduced herein as Appendix 5. The inclusion of space-charge distortion 
does not change the fact that thermal breakdown is a high temperature phenomenon. 
Consequently it is of little interest to the spacecraft discharge problem. 

The overall conclusions from the initiation mechanisms study are 
that the conditions in space are conducive to low differential voltage break- 
down and that huried change is probably involved. 
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5. QUENCHING AND BLOWOFF MECHANISMS 


Laboratory charging of spacecraft dielectrics with electron 
beams has led to an extensive data base on the subsequent breakdown. Two 
features of the observed process are: (i) the amount of charge cleaned off 

in the process is approximately proportional to the area of the sample 
[Balmain, 1979] for sample sizes up to 1 m^; and, (ii) the discharges are al- 
most always accompanied by the release of electrons from the front of the 
dielectric into the vacuum. The electrons released quickly escape to the ex- 
perin»ent chamber wall and manifest themselves as a return current to the 
sample substrate. The magnitude of these blowoff currents is large, involving 
as much as 50% of the total charge released there (the remainder existing via 
other current paths). In addition, the peak current scales approximately with 
a characteristic dimension of the sample [Balmain, 1979]. The first of these 
observations is nonnally called surface cleanoff. The second is described as 
the blowoff current. Certain theoretical aspects of thess two phenomena are 
discussed ^elow. 

The fact that laboratory discharges appear to involve charge 
from most of the sample surface suggests that once a discharge has initiated, 
it is self-sustaining, and does not cease until the source of energy for the 
process (within the field due to the trapped charge) has been relaxed. It is 
of interest to ascertain whether the negative tip streamer described by ACORN 
has properties which are consistent with this behavior. 

Before addressing this subject, it should be remembered that 
the negative streamer is not envisioned to be the discharge. Rather, a 
Markovian sequence of small streamers, initiated near the positive ends, 
is expected to provide the mechanism for discharge and propagation. Because 
a pre-'.ise mathematical formulation of this process has not yet been given, 
the discussion of the quenching of the discharge given below must be viewed 
as speculative. Explicitly, the discussion is concerned with possible quench- 
ing mechanisms for the negative tip streamer described by ACORN. A conclusion 
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about the behavior of the negative streamer can become a conclusion about 
the behavior of the entire discharge only If It Is the quenching behavior 
of the negative streamer which controls the propagation. If there are slower 
processes which control the stochastic Initiation of subsequent streamers, 
then the above conclusion does not follow. Of course, If any mechanism 
exists which tends to quench the streamer Itself, It does follow that the 
entire process Is quenched. 

Within a streamer, there are a number of processes which could 
be expected to lead to a quenching of the propagation. These processes are: 
diffusion, trapping, recombination, encountering a region of field reversal, 
encountering a macroscopic or microscopic material boundary, and global 
field decay. Each of these mechanisms Is discussed In turn. 

Diffusion Is already Incorporated Into the numerical nodel Im- 
plemented In ACORN. With the Inclusion of this effect, an avalanche appears 
to evolve Into a propagating, self-sustaining, Ionization front. While the 
effects of avalanche are smeared out due to diffusion, the effect does not 
appear to lead to any mechanism which would stop the avalanche process. 

Trapping Is not expected to be a viable mechanism for quenching 

the avalanche process. Shown in Figure 5.1 are the ionization times for all 

the models which have been considered in earlier sections of the text. Typical 

•ft 

trapping times are of the order of 10 sec. [Hayashi , et al . , 1975] or more, 
depending on whether shallow or deep trapping Is considered. Thus, In order 
for trapping to be a competing mechanism with avalanche, the avalanche times 
must increase to values which are as large as the trapping time. From the fig- 
ure, it is seen that this occurs for much smaller fields than those of interest. 
In fact, one may turn the argument around. If trapping times are sufficiently 
small as to compete with avalanche, then the avalanche will not initiate at all. 
Once it initiates trapping cannot be of concern. 

This conclusion may be reached from a slightly different viewpoint. 
It was seen In earlier sections that once an avalanche Initiates, It quickly 
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reaches a critical charge density given by Assuming no 

cancellation by positive charge, a purely negative charge leads to a 
critical number density = cpQaQ/e. Shown in Figure 5.2 is a plot of 
this quantity for all the models considered in earlier sections. Typical 
trap densities are of the order of '0^®/cm^, with a maximum of perhaps 10^^/ 

3 

cm . With the densities shown in Figure 5.2, it is clear that even if the 
trapping rate were very fast, all traps would quickly fill, and there would 
still be adequate charge in the conduction band to continue the avalanche 
orocess. 


Recombination is somewhat more difficult to quantify. There are 
several reasons for this. The first is that recombination with the parent atom 
(so-called geminate recombination of Onsager [1938]) is somewhat difficult to 
observe. Second, most observations of recombination have been performed for 
relatively small carrier concentrations {« 10 /cm ). In this regime, re- 
combination is dominated by trapping processes. That is, recombination occurs 
between a mobile electron and a trapped hole, or vice versa. This trap control- 
led recombination tends to give rise to a kinetic theory term which is linear in 
the mobile electron density. At higher carrier densities, typical of those noted 
above, direct two body recombination is expected to dominate (the kinetic theory 
term is proportional to the product of the mobile electron and hole densities and 
increases quadratically with the carrier density). It is possible at the very 
high densities associated with Figure 3.2, that even higher order processes may 
become important. This has been observed in high dose irradiations of semi-con- 
ductors [Van Lint, et al . , 1980]. Finally, all the recombination processes are 
sensitive to the electron distribution function. For the high fields considered, 
significant distortion of the distribution function takes place. This shift to 
a heated electron distribution function tends to make recombination less likely. 

Even with all these caveats, it is probably still safe to say tSat 
recombination will not quench an avalanche once initiated. Assuming a term in 
the rate equation of the fonn r*n*m, where r is a rate per unit density, n is 
the mobile electron density, and m Is some hole density (whether trapped or 
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free), the rate r is related to a capture cross-section o through a 

formula r ■ <ov>, where < > indicates an ensemble average, and v is the 

electronic velocity. The largest quoted value for o [A. Rose, 1963] is 
-13 2 

of the order 10 cm . This large value of the cross-section occurs 
at very low energies. At higher energies, o falls off rapidly. The maxi- 
mum value for the recombination rate r occurs at low energy. For electrons 
of near thermal velocity at room temperature, the mean speed is about 10^ cm/ 

sec. With these estimates of a and v, the recombination rate has a maximum 

-63 18 3 

value of 10 cm /sec. Assuming a hole density of 10 /cm , gives a minimum 

-12 

recombination time of 10 sec. Even this worst case estimate gives a time 
which is larger than the characteristic avalanche times as noted in Figure 5.1. 
In practice, it is expected that the maximum recombination time is substantially 
larger. For example Gross [1978] quotes a value of r = 10 cm /sec for Teflon. 

Therefore recombination cannot quench the processes which are 
occuring in the active region of the avalanche, the ionization front. Once 
initiated, it will continue. In the body of the streamer, however, it might be 
possible for recombination to reduce the conductivity. In this region, the 
field drops substantially from ambient levels, and there are significant free 
carrier densities present. If nothing further happened, except for the propaga- 
tion of the single ionization front, then the recombination process will take 
over and decrease the f r 'e carrier concentration. Further quantification of the 
importance of this process requires a precise formulation of the discharge in 
terms of a succession of small streamers. 

The question of whether a region or field reversal will impact 
the negative streamer formation has a strai<jhtforward answer. It is topolo- 
gically impossible for the field reversal region to impact the evolution once 
initiated. The reader is referred to Figure 3 of Appendix 4 which shows the 
field configuration near the front surface of a dielectric irradiated by a 
space-like electron environment. It explicitly shows the field reversal of 
the postulated mechanism [Meulenberg, 1976]. To the right of the zero of the 
field, the electric fie^d intensity is negative. Electrons and subsequent 
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streamers In this region will therefore drl^t to the right and never encounter 
the zero crossing. Similarly, to the left of crossing, the field Is positive 
and electrons and streamers will move to the left In this region and never 
encounter the zero crossing. The point is simple. Electrons are trying to 
escape the region of trapped electron density. There Is no way to configure 
the boundary conditions so that drifting electrons will encounter a region of 
field reversal. 

The most difficult processes to quantify are those associated with 

encountering a physical boundary, whether macroscopic or microscopic. For this 

reason, the following discussion is primarily intuitive and heuristic. Clearly 

If the negative streamer encounters a conductor, then the avalanching process 

will cease. Subsequently, the region of space which has been left with mobile 

carriers will rearrange itself so that the surface of this region Is equi potential 

with the conductor which was encountered. For the conductivities computed in 

-13 

Section 3, this happens very rapidly, 10 sec. The primary process has 
been stopped, but the net effect has been to project an "effective" needle-like 
conductor into the dielectric. Assuming the potential has been maintained, and 
the positive tip is nearer the cathode, the field continues to be enhanced at 
the positive tip. If nothing further happens, the process has ended. A repeti- 
tion of the avalanche mechanism must be invoked for further spreading of the 
process. Further statements must await a formulation of this mechanism of 
breakdown. 


On reaching a vacuum boundary, :t is also clear that the avalanch- 
ing process must cease. The disposition of the avalanching electron front 
thereafter Is not quite so clear. Two things could happen. The front could 
reach the surface and encounter a potential barrier which stops the electrons 
at the surface. Alternatively, the potential barrier could be sufficiently 
small that the electrons are energetic enough to penetrate the barrier and 
escape to free space. Of course, the experimental evidence suggests that 
the la^’ter occurs. Theoretical arguments to support this conclusion are more 
dif T*’ to muster. The electrons have at most a few eV of energy above the 



conduction band edge of the solid. The relevant barrier Is that of the work 
function, or electron affinity. It typically has a value much less than the 
work function of metals. Values of 1-2 eV at most may be expected. The major 
uncertainty concerns the effect of the deeply trapped space charge (the internal 
electric field) on this value. It is not absolutely apparent that all aspects 
of this field have been accounted for by merely computing the bulk electronic 
distribution function In the field. Indeed enhanced secondary emission of 
charged dielectrics has been reported from time to time suggesting that some 
lowering of the work function barrier may occur. 

In any case, the barrier is at most about the same energy as a 
characteristic energy of the bulk electron distribution function. One there- 
fore expects a significant fraction of the electrons in the front to escape. 

The exact escap? ratio depends on the details of the distribution function 
and the electron affinity. A detailed formulation and calculation is outside 
the scope of this effort. 

The effect of the above Is that the avalanche ceases, a short 
burst of electrons Is emitted into the vacuum, and an enhanced field conduct- 
ing region Is left behind. Further evolution depends on a successirn mecha- 
nism for avalanches. 

The effects of microscopic discontinuities have not yet been 

Investigated. 


The dependence of the solutions on global field decay follows 
from the discussion of Section 3 and the trapping and recombination discussion 
of this section. From the scaling laws of Section 3, the characteristic time 
for avalanching is 6 q“^. The quantity increases exponantially with a decrease 
in field. If sufficient charge has been released so that the bulk field is de- 
creased significantly below the initiation field, then the trapping and re- 
combination times can begin to dominate, thus providing the required quenching. 
It should be mentioned that an exact quantitative relation cannot be given 
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within the present format f r the field at which this occurs. As the field 
decreaseSi and trapping and recombination become dominanti the assumptions 
which were used to solve the transport equation for electrons become invalid. 
That is, it was explicitly a^isumed in the formulation that the field was suf- 
ficiently large that the quasi-free electron approximation could be used. 
Transitions to localized states were to be treated as a perturbation. The 
limit of interest requires that trapping and recombination become dominant. 

In this limit, it is the conduction band motion which is a perturbation on 
the trapped states (sometimes called hopping conduction, or trap modulated 
drift). The effect is that significant ionization shuts off even faster than 
implied by the exponential behavior. What remains clear is that this 
mechanism will turn off the discharge, but only after a significant fraction 
of the charge has been lost from the insulator. 

A review of the paper by Leadon and Wilkenfeld [1978] leads to 
the belief that magnetic forces are probably not important in the blowoff pro- 
cess. A review of the numerical values presented therein is sufficient to 
come to this conclusion. In order to obtain a magnetic force even remotely 
close to the ambient electric field, the following parameters were chosen; 

Q 

electron drift velocity in the dielectric, 10 cm/sec; peak discharge current, 

3 3 

6 X 10 Amps, and distance from punch through channel, iO meters. With 

these parameters, the electric field equivalent of the force is 10^ V/cm. 

These values are probably not appropriate. The drift velocity 

used is approximately two orders of magn'tude too high when compared to those 

computed frcmi SEMC, and measured under a very wide variety of conditions 

The discharge time (10 ns) assumed is significantly shorter than observed 

for the assumed sample size (10 cm). A value closer to 10^ ns [Balmain and Hirt, 

1980] is more reasonable. Thus, the assumed peak current is probably two orders 

of magnitude too high. Finally, the magnetic force falls off inversely with the 

distance from the discharge channel, so that 1 cm from the channel, the force 

is reduced another order of magnitude. Since it is Impossible to emit the 

amount of the blowoff current which has been measured ('v< 50 A) for samples of 

-3 

tnis size from a spot 10 m in radius (see discussion below about space-charge 
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limiting), it is necessary to look at larger distances from the punch- through. 

These considerations imply that the importance of the effect is an estimate 
of at least five orders of magnitude too high for the example presented. It is 
concluded that the mechanism is unimportant. 

More generally, if these forces were important, then their effect 
would be strongly dependent on sample siz? (since the discharge current is). It 
is difficult t'' conceive of circumstances which would allow scaling of this 
force (it decreases with current) in such a fashion as to always dominate the 
electric field force. More generally, magnetic forces are seldom of importance 
in transport processes unless both electric fields are shorted and inertial 
effects dominate. 

Two other mechanisms can be envisioned for the blowoff process. 

The first of these is the termination of an avalanche at the surface, or, more 
generally, the punch through of a discharge to the surface, with the sub- 
sequent emission of electrons. It is important to distinguish between large 
scale or macro-discharges and small scale or micro-discharges in considering 
this mechanism. While it is possible that this mechanism is important for 
micro-discharges, it seems very unlikely that it is responsible for macro- 
discharges. 


The reasoning is based on space-charge limiting. The large dis- 
charges which have been observed cannot possibly have been emitted from a small 
area characteristic of a surface punch through. Under space-charge limited 
emission conditions, the maximum current which can flow (I) for an applied 
voltage drop V is given by the Langmuir-Child law I = where G is the 

perveance of the emitting system. For the purposes here, a detailed discussion 

2 

of G is not required. Instead, we assume a 10 cm sample which has a grounded 

grid 10 cm from the front surface of the dielectric. The perveance in this 

-2 -3/2 

case is G = 33 mA cm kV ' . For a 20 kV potential drop, this gives a maxi- 
mum space charge limited current of 30 Amps. This estimate is based on the 
assumption of uniform emission from the entire surface. The total allowed 
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current will scale roughly as the actual emitting area. Thus, a small punch 
through with an area of 10’^ cm^ (which Is large for a punch through), can emit 

nearly four orders of magnitude less current than is observed. We do not pursue the 
quantitative aspects of this line of thought further. It is apparent that a 
few small punch through channels cannot emit enough current to account for the 
observations of current scaling for macro-discharges. 

By contrast, micro-discharges with very small total blowoff cur- 
rents could be sustained by this very same mechanism. A surface directed aval- 
anche/streamer can arrive at the surface, and the free electrons can be ejected 
from the surface. This was discussed earlier In this section. It seems quite 
likely that this type of small scale charge emission occurs routinely In charged 
dielectrics. The details of this mechanism should be developed for describing 
small scale scintillation (probably without damage) of charged spacecraft di- 
electrics. 


The final mechanism which has been considered Is the diffusion 
and release of free electrons generated In a subsurface propagating arc discharge. 
The Idea Is quite simple. Based on the conceptual model of a stochastic sequence 
of avalanches, It Is expected that streamers of the type described by ACORN will 
propagate through the subsurface trapped charge layer, with the direction of 
propagation parallel to surface. Tvldence of this type of subsurface propagation 
has been given by Balmain [1979] from observing damage paths In the material 
subsequent to breakdown. These discharge channels are expected to branch out 
into a tree, which is similar In form to the damage patterns. Filamentary 
processes are therefore expected to lie under the entire surface area. 

Each of the filaments may be described by a conducting channel of 
the type computed by ACORN. In these channels, a large density (see Figure 5.2) 
of free electron carriers is generated. These electrons will diffuse outward 
away from the channel and a portion of them will further drift to the surface, 
where they will subsequently escape. 

It Is to be determined whether this mechanism Is quantitatively 
reasonable. The mechanisms which could prevent the electrons from reaching 
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the surface are trapping and recombination. Because the holes are relatively 
Immobile, and the diffusion Is away from the regions of large concentrations 
of holes, recombination Is not expected to be Important. Trapping Is Included 
In the following estimate. 

Assume that a quantity nQ/cm electrons are released at depth d 
below the surface In a planar configuration. The process of escape to the 
surface may then be followed with elementary diffusion In the presence of 
absorption (trapping). Let t be the time constant for trapping, and D be the 
diffusion coefficient. Define a characteristic distance L by the equation 


L = >/Dt • (5.1) 

For each electron released at a depth d below the surface, the fraction (n/r^) 
which reach the surface Is 


^ = exp [- f ] . (5.2) 

"o 

Some numbers are instructive. From the high field computations of 
Section 3.1, D has a value of about 10 m /sec. Assuming a trapping time of 
10"® secs [Hayashi et al., 1975] gives a characteristic length of 10"® meters. 
Typical depths of the charge layer are a few microns, so that the exponential 
In Equation (5.2) is not extremely small. Since a copious supply of electrons 
will be generated in the discharge channels, the fraction reaching the surface 
will be more than adequate to account for the observed amount of charge released 
In the discharge. 

The time T required for electrons released a distance d below the 
surface to reach the surface Is given by 


T = d^/D • 


(5.3) 
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Choosing d » 10 m gives T = 10 secs. Compared to the time scales of 
macro-discharges, this is quite fast. The rate of formation of discharge 
channels (which would control the release rate of diffusing electrons) 
would therefore control the current reaching the surface. The model is 
consistent with the blowoff current having the same pulse shape as the 
propagating subsurface current. The scaling with area, which has been 
observed, is consistent with this process if the subsurface discharge has 
the anticipated scaling behavior (peak current proportional to length). 

The conclusion is that diffusion to the surface of electrons 
freed from a subsurface discharge is a very likely mechanism for the blowoff 
current production. Note that emission is expected over most of the surface 
area in this mechanism. The released current is not so severely limited by 
space-charge considerations, and is comparable with observations. 


6. SUMMARY 


A description of continuing work on the development of a 
detailed model of arc discharge dynamics has been given. The work de- 
scribed Is centered around an extensive parameter study which was per- 
formed with the SEMC/CA3CAD/AC0RN family of computer codes. The salient 
results of this effort are summarized as follows: 

• Code modifications were Incorporated Into the SEMC/CASCAU/ 

ACORN family which permit reasonable code to code Interface. 
Adequate output routines were also included In the package. 

• Dimensional analysis and an extensive confutation?.! study 
with the SEMC code have provided a characterization of the 
single electron distribution function and associated transport 
coefficients over a wide range of material models. 

• Analytical approximations for single electron induced un- 
distorted avalanches were compared to an extensive set of 
computations of the same phenomenon using the CASCAD code. 

The problem parameters were those obtained from SEMC compu- 
tations. It was noted that self-fields become important 
when only a small number of particles are involved in the 
avalanche. The transition to a continuum description for 
the self-consistent further evolution must therefore be 
interpreted as an ensemble averaged description. 

• An extensive series of computations of the self-consistant 
evolution of electron avalanche/streamer formation suggests 
that this evolution can be described by a self-similar 
solution of the macroscopic equations. The scaling dependence 
of the solutions on almost all of the parameters has been given. 
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• A model for low voltage differential breakdown in 'ii- 
electrics was given. The electrons from the environment 
penetrate the dielectric and give rise to a trapped charge 
layer in the dielectric. This layer gives rise to large 
internal electric fields near the surface which cannot be 
reduced by enhanced secondary emission. These fields can 
become large enough to induce breakdown. The space environ- 
ment gives rise to the field configuration of the proposed 
mechanism. This is in strong contrast to monoenergetic 
laboratory simulations. 

• The likely quenching mechanisms for terminating a discharge 
appear to be associated with encountering a vacuum boundary 
and/or global field collapse due to a release of sufficient 
charge. 

• The most likely mechanism for large blowoff current appears 
to be diffusion to the surface of free carriers released in a 
subsurface discharge. 
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Reproduced From: 


Lin, D.L., Electron Multiplication in Solids, Physical Review 
B20, 5238, 1979 
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Electron multiplication in 5olid» 

Dofif L Lifl 

icirnrt liff’luunont tnforpu'ainL iJJO OiJ C^tiUtiouu HmJ. iiuiit iiO. IVnnu. I'lryinHi 2JIH0 

(Rwn»«l « July 1979) 

Trrdlititi 4» Maikot pinuMO. ■ (brrir* u drvelutvil lot llir c.tli uliiiion ol rl'il'im 

nmliit'licaiinn in v'liii^ 11i( rlft'iion mean fiw ptiKj t>3t tcousiii (<lnmon Kai.rimy. iipiual pliomm 
rinis^iofi. ind pair pitklm non pK r\tc* tit tUumeO to Ik indcptndtni of llit tnti|t> of the tirciiori tlir 
nuM.I»i c.f loriitiM'ii ftrnit |itr metn fiat pttb ii i unotrtol fumiini ,if dinirn\ioiite»i torryt ijiuk 
S rttitl i\(<n of pliof.nii viiir'inyt art iluditd In ibt caw of laiitrnpic uttltiin|i, cnitip.tnu>n of ili>- 
pirt.rnt uofk uilli Iht itNuli' ohiaintd by numrnrall) lultinf the lloli/niaiin ri|UNiiiiii indicairs lliai yi*«l 
aiiirrmrnl ii louivl only m f.itl-inuliipliuilion ca«A llie diwiepanc) noted in klou-mnlnplKaliiiii ..isrN it 
lK.lirtrd 10 be due lo ruber the Irvt iinnfrnl leuuurmtnu or ihr mild tingulaniy of ilie iiaiiMiioii 
liioliabiliiin in the nunieru^l appiuacb to aol»in| iht Ouiumann tuuaiion. 
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I IMfUJIlUniUN 

Under (be stresK of a verv strong electric field, 
gt'iniroiiduclnrs and dicloi-irirs lend to exhibit ava. 
lanchc breakdown. A qu.witiUtivc nicacurenicnt of 
the (liarge niul(i|}|icatir>n in l>-» junr’ ons of 5i and 
Cif was earned out and it.s implicam,i on the ion- 
izulion rale (icr unit length, as a function of the 
.itipltfd electric field, w a.s analyzed Iv-odecadea 
ago.' I he fundamental pmcesscs for rharge mul- 
tiplii’aiion in snlids were asbumed to be analogous 
to titii fl aw nsend's fi niechanisni*} tif gas dis> 
charge. The ionization rate is an important para* 
meter which brings nut the details of the micro* 
scopic solid-state properties from the macro- 
scopic characterestics of breakdown nieasure- 
nients. Methods have been developed to calculate 
tilts parameter, most of them Involving the solu- 
tion of the Onitzmann equation in a high electric 
field. Thus, neglecting the then- unknown band 
structure of silicon, Wolff’ expanded the electron 
distribution function in terms of l>egendre polyno- 
mials, kept the first two terms, and solved the 
IinUrniann equation in a steady state. Uaraff em- 
ployed the concept of collision density and derived 
an integral equation by Uiplace-traiuforming the 
Holtzmann equation which he then solved numer- 
ically. To solveUhe Doltzmann equation at high 
fields IS generally difficult, as pointed out by Wan- 
nicr.* On the other hand, Shockley* considered 
the collision procesbcs as probabilistic processes 
With exponential probability distribution and, 
treating only the electrons which survive any col- 
lision, directly obtained the innization rates in 
llic low -field and high-fuld limits. Although the 
validity of some assumptions made by iihockley 
was questioned at that titne, the concept of expo 
nenlial proti.ibi tity is partirul.irly attractive in its 
simplicity. It appears that for the calculation of 
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Ionization r=tles, one may avoid solving the Holtz- 
mann equation by treating all prncessc's as Markov 
processes; and following the tnrrfion of the eleetrou 
in a completely stochastic way. 

In this paper, we will addrcs.s Ihe c.i Icul.ii hui of 
ionization rales in bcmicondjciors or dielfctrics 
where band-Blrurlure^delails can lie neglected. A 
strong electric field S is :ip|ilied to a (hick i l.iti of 
material. The clecirnn whose li aek we w ill fof. 
lou starts wi' zero energy in the l•oMducllnll imul, 
Wc Xisume Uul there are three energy raiit’e.^^ in 
each ol which only one sinelp.pU-rirno process is 
possible. Thus from zero to /iw Itiie energy ol the 
only optical phonon), only aeoustic- phonon scat- 
tering IB possible; between and £« (the lomra. 
(ion threshold to produce an clcrlron-hole pair), 
the electron can only emit the optical phonon, as- 
summg no optical phonons to be present so that 
the absorption of optical phonons is negligible: 
above £,, electron excitation from the valence 
band la the only interaction the elociron has with 
Us environment. We also assume that the mean 
free pathu of these three fundanieni.il processes 
ac* independent of the cleclcon eiieri'v Further- 
more, these throe mean free paths are assumed 
for simplirily to have the same value a. These 
assumptions are in line with previous works. *•*■* 
Moat of them can be relaxed il one so wishes: the 
involved modifications ol the equations, to be 
derived below, either are trivial or can be car- 
ried out with minor effort. 

The main quantity we are interested in is the 
mean Ionization distanee lor an electron starting 
with zero kinetic »*ncrgy at z = 0 (r is Ihe depth co- 
ordinate). An electron, emitting » opiical pimoons, 
will reach £, at location / = (ft, ♦ nfl-’)'riJ. Since z 
and n are two random variables. Uie mean dis- 
tance I per mean free path is related lo tin* mean 
number of optical- phonon emissions, ii, l>y 
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OF POOR QUALITY 

M. I I I k U > M I, I I n M «; A T lo N I > S i» I. M» ^ 


» K ^ J , fS A * «(/)* ,>A' I ) (I) 

itcUinc .V' /•,'•’(?» and /< fu' l\ . ili#' mean nuiiw 
luT n( i'liM , rut ■ hn|t< pain- (rnduc*^ li> a 
t'lc'Clroii III lr<urlin»; a iiii'aii fret* path t iv civcii 
•>> 

(v.fifa « U'' . (2) 

1 III' I Mil lilt- ri^hUliaiid Hide Uikcii ijrv ol lie 
• lira nii’.tii Jiiilarre n, ltii> I'iiiiron liastnlravcl, 
UK' I rracliinc in nrder in (.rnduiii 3 pair. Tli«‘ 
i’».*rl value ol lliib tn!ra ii e.in iunUiii'p Ih not iin- 
porlaiil , ’■ * hiir evi- r Heir wr .usmiMiC a,» k 

l,(l.).iliiii| (2' intiicjlrii llkii if> a uiuverFai 
luiK iiMit (i| ,V and R fiince if, ai> vi ill U> Ehown la- 
ter , drj«'iHlh only on ,V and /?. 

Tlir criirr(.l drsrriljed hri" could l»' a()pltrd to 
t*ie calrulalinn ol the avcrai;r enerpy n( hot elcc- 
Irniib rmerttuii’ frnii. a KUb ol civvn Ihicknebk 
The prqjer equation reads 

r fifnt (3) 

and >7 could be oiriained exactly in the saiiie uay, 
an described in See. 11. The only change in (n re. 
pUcr y - dv in the ariruntentft |jy », Though Monte 
Carlo rrsultK* do exist, the .suggested approach 
here gives an exact answer. 

II THCOKV 

Let * be the mean (ree path o( a collision (a Mar- 
kov process). The no-cnllision probabibty for an 
electron in travelim; a distance f is e'*^ and the 
prolubilit) ol one collision in tU ihtiz/K. The 
mean collision distance is given by 



the m*»an (ree path of this prnci'.ss, 

At some reference level, an electron is released 
with both kinetic and potential iiicrgies zero. Un- 
der the influence of a constant electric field S, 
this electron bcluves like a particle of some in- 
ertial mass In a gravitational field Now, if this 
electron cxperleAces an acoustic (elastic)coliision, 
It simply adjusts the direction of its velocity with 
magnitude unchanged. If an optical -phonon cmis- 
aicn occurs, the electron not only changes the 
velocity direction but also loses some quantum 
kinetic energy (the quantum energy of the op- 
tical phonon). Since in a conservative field, there 
is a one-to-one eorrespo.ndence between kinetic 
energy and potential energy, once the reference 
level of zero potential is chosen, the (act that the 
electron loses flu> kinetic energy ivii) Oc equiva. 
lent to pushing this electron .a otsunce r,^,rf clo- 
ser to the reference level, as shown m Fig I 
The probability of having a rolhsion ix controlled 




»l» 


••i 

K 




1, 

till t lllri Iron ai /I eiiiiiii (III oi'lK'.il I )>Mik>n aixl Iin 
locatiuo It rxiiwequently put al A, da • h^/nk. 


by the oxponniiial (unclur, ol the lenr:ih of the ira- 
jertory thia electron travels after the latest col- 
lisioo, 

For a Markov process, the history of the ctec- 
tron before a rnllision is completely wiped nut 
when a col'.ismn occurs. In oilier vvordb. thi angu- 
lar di'itributinn of the cross section ol a rollinnn 
in independent of the dirrctirio o( this inculeni elec- 
tion A few examples arc m order: 

(a) Forward scatter, ng, Tlie electron can only 
travel in Die direction of the electric force. There 
IS 00 elastic collision in this case because the di- 
rectinn is lixed. 

(b) Two-direction scattering. The electron can 
move either along or against the electric field. 
Movement in cither direction can cause collision 
v.ilh the conb(<quence of, say, lulf chance ol the 
electron going along and half chance of going 
against the electric field. 

(c) Isotropic scattering. The cross section is 
isotropic and thus of Markov type. 

in general, we can assume the angular disiribu- 
tion of the probability to bo P[0) no matter what 
direction the electron is moving in immediately be- 
fore a collision. 6 is the angle the velocity vector 
makes with the horizontal line as shown in Fig. 2. 
The electron is at a location y, (v xa/X) below the 
reference level. I.ol 7*(y,, v,)dVj Ik* the probability 
distribution that the electron, suffering a rolliaion 
at y,, will not have any collision before reaching 
the height V| below the reference level and having 
a cotlislm in dv, at v,. The collision here could 
be acoustic nr optiral. This probability is compu- 
table from classical mechanics, as can be done 
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fnllrtwinc (hr rrsutu of (hr A|i(M'fufi«. (( in «ai>y In 
priMf llwit f,T7(\ ,, v,V/v,» 1 if rcali/«*>> ilui, 
for j jsivfii projorlilc anr.lt, iht* pioiutiilii\ <t| hav- 
ing’ a ciiilr'iun r dloNb n( hrm Ur th<> ('li cti <n 
Iravrlb II. iinilv, i t I . Oner P(fti m 

aubit'iii'rf and 7(v,,\,l is rntnpulid, (hr loniralum 
rau> fan lx oUaini."l aH fnilou ». 

I.i'i 

- f ,ki\ ,* f T(v,. v,)«/v, 

■'o •'»,) 

•if l(u‘ prolMbi lify than an elfciron suffcrini' a enU 
liMort ai I, uiM not havr any rnllibian btfort rc.uh- 
*'H5 ''o ^ inniaatinn Ihrtsliold for prorfuc- 

in*! an I'lcclrnn-holc pair. Lot Ay bo Ifio width of 
the uc oust ic -t>r at I erint; region, i.e., 

Av - f5Aa/.’5» - fiui/eSx^RX. 

We Can llion rnmpute, by tfic technique of reRroup- 
inK, the probability of having w optical -phonon 
eiiiib.->innh , repurdlosN of the nuiiil>or of bcattcr- 

1 


1 = \(0)» f .V(v,)7'(v',, OV/v, * { </V| r ilV;T{v,,\' 

/g ''O ^0 

♦ J ’ ’ V , *® d V , r tv „ ,V , - A V ) r(v , . 0 ) . 


if (»••"•) denote! the number! of aroustic and op- 
tical Rcaltcrings. then the first term in Eq. (6) is 
f0,0), the second (1,0), and the fourth (0, 1). For 
the thiixl and the last, we need to break up the in- 
tegration /pO and (hen aubettlute the appropriate 
equaimii (nr 1 as before. This procedure could be 
repeated indefinitely. Collecting the terms with 
(n,. 0> for It, aO, 1 , 2, 3, * gives 

At(0)* f*’dv,Ntv.)r(v,.0) 

Jo 

♦ j *dv,J dVjhffVjinvp.y.jriVj.O)* • •• 

-AT(0)4 r‘dv,W(v,)r(v,.0) = .Vf0) , (7) 

where A'(v,) satisfies the lollowing Inhomogeneoub 
integral equatio.n of the second kind: 

pV^V(v,)rtv,.v,). (8) 

Jo 

We CJU 1 Intrrpretr at Hie rrnorni.ilirod ttiy 
intf rarlmr. with acoustic phnnnnr.' vctmou o( the 


ings due to acoustic phonons. 

Tlie electron releaKrsl at lero kinetic energv at 
the poti-ntui reference level reachci. »■., cilhtr 
uilhnul any i nlliiiion or with at leasi one cfilli‘.ion 
I acoustic or nplic.il). In other words, 

1 • MOW J I X r(y,,0Wv, . 

After bre.ikmg up the inlf):i >il inn luiiiis ii.tolhc 
aioubiic ao'! of)(n.il rivoos. the I oo th>' iight- 
h.ind suit (Mil tic iiplai«‘d tiv 

Vfv,) ♦ r ” 7(v j, v,l/v , 

Jo 

when V, IK in the acoustic ficattering region, or liy 
.V(v,.. Ay)* f**T(v,.V, - AvWy, 
when V, is in the o()tical cmissi'Mi region: 


)r(v,, 0 ). 



f/y,.V(y, 


A\ )7(v,. 0) 


( 6 ) 


I 

bare quantity A'(v). 

For one optical -phonon emission, all the terms 
of the type (x,, I), «, » 0, 1 , 2 need to be collec- 

ted and regrouped. After two regroupings and one 
relabeling, the probability of one optical cmiHsion 
Is given by 


f’dvVCCv, -AvlflVi.O), ( 9 ) 

where akr.o satisfiefi an intoftral equation, 

rtv,.>'i)-7‘tv,. v,)* f dv, nv„y,)f(v,,y,). (lO) 

Jf 

Again, due to 4 he presence of the acoustic -scat- 
tering region, the bare transition prnliability den- 
sity 7"{v,,v,) IS rcnnrmatir.ed to /(v,,\,f. 

A similar procedure leads to the following prolu- 
bility for two optical- phonon emissions. 





Avinv,. y, 


-^v)7■•(v,.0) 


(in 


Ti e rule for writing the probabi li tv of anv nu:i' • 
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bor nf al- can tic 

by b() Ul) /')>,, 0) Hidti airs tikil lltr 

t'irciron is i rli'aM'd at i - 1>, nutting; alt pnsHibli* 
inti’iMCliniis U'llh tlir »•' ustir pliomtiiN, anil rrai li 
fs I, Ml till' optical t'lnissmii itv,mn Kniiltint; an 
rtplii.il |ili<>uon ut I , am! llii'ii liv < IwtnttM^; ils In. 
t.ilumiiii,. Jti , till' I’b'i 1 1 'ki I lu oliiiicrs all pns 
sil'li'. .ncluilinp i'lMd, oiiii.lMTS oI acouMic-pliuooo 
tiCaliiM'inph liHotc tlic iicM optiial (iIioikui pmiis. 
smn at in the (ijitic.il uoissino rcpinn, a*, dcs. 
cnl>«‘il by .ill. It (■ (i|>lical-|)hnmi'i cmib- 

i.um at V. puiilicb tlic liKiili'm ol tbr ^■lc(tr<*t In 
V,- Av. nclui’in tl.is ncM tsilion anl Ibc pair- 

) 
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|irn(luclinn lbrc^^lol.t no •■•jilical cmihMon occui't« 
while any iiuiiibcr r.i acoaslic scattcrmitK iii.iy take 
pLici'. Tliib lb rcpicsi'nted tv \\\,- &v). ficoiJif^; 
t'i|uatioiis like yt|. [1 1 1 ih Completely equivalent to 
deseritiint’ the physical processes they represent 
II IS t;aurai.li ut, lion; Uica.i. tiie .ibovcproh.i- 
bilil les a 1 e ilcnveii. llsil the suin ol all proljabili. 
Ill /• equ .1 Is ijiiily . 

The avera);e liuiiitv r ol optical phonon ennssiom 
IS piven by llii' snin o( the piodurts nl the iiumlier 
ol iiptii.il ptiunon.s emitted and the ('oircsp<>iidinp 
protsitiility , 


ORlQ^At HAtifc idr 
OF POOR QUALITY 

tllt.lHtih V I t I M* 1. 1 1, a I Ml h I '• '1*1 III'- 


ri- 


0N(0)elj”V,.V('.- Av)7\v,.0)*2 j *</v, A'tv, - A>> 7 (> ,. V, - Ay)ftv„ 0) 
♦ 3 f , - Av) f *</v, f(v,, V,- Av) ( *</v, 7’lv,,y, - ^y)ti\^,0)* • •• . 

efc •'at •'si 

I 


( 12 ) 


l.ayiiip down the protiubi I um'S with uiMphls (the 
number ol phnnnos emitiivl) at. in Fip. 3, and suni- 
miiii: veitically tor each i>.|umn. one obbilns a 
(unction C salistymp 

6'(v,) = 7*(v,.0)4 r®dv. 7'(v,,y,- Av)G(v,) (13) 

for e.ach entumn When all the columns are 
Burnmed, the liiv’l iapn?.sr.ion (or the averaRC 
number o( np'>cal- phonnr. emissions is given by 

r®dyA\v_ Ay)//Cv), (14) 

*'4f 

where 

W(v,)-G(v,)t r’®dv,f(\,,v,- Avl7/(v,). (15) 

9 ■ «- Osf (JMICAi IntVSICN 

• • «- Us Orric*!. IniisiONS 

• • • uc. 

* • • • • 

• • • • • 


1 11 . 

; i ; 

U (• - - • 

rii’i I Snmrrii"' a tolnrin cues il Si.e iiiinj' -all t' f’ 

iMven H. 


Equation (14) indicale.s II is a function of v„, 
whicb IS equivalent to the energy ratio Xr t ^/cSk, 
and Av. another energs’ raiin, (e W.V) 

lliveii H and X, is uniquely dclerimned by 
I'Ai (I,. 


Ill RESULTS ANU bISCTISSIUN 

As mentioned in See, II, there are three interes- 
ting cases of icatterings: forward scattering, 
two-direclion scattering, and isotropic scattering. 
The iirst rase provides a check of the formula- 
tion described above against an algebraic approach 
to be published elsewhere.' The second case gives 
a hint of hou- the backward scattering (agaujst the 
electric force) would change the mean itmiaation 
distance, while the ronipbcation of calculating 
the trajectory length is kept to a mmumuro. The 
third case is that of the real physical situation 
which las been discussed by others in the litera- 
ture. 

Case (i): forward scattering. 

jO for y,<y, or y,<y,«Ay 

nv„y,)-j (16) 

^ exp - |y,- max(y,, Av)| otherwise. 

Backward scattering (v,< v,) is forbidden In this 
rase. The \(v) corresponding to this transition 
protnlulily is given by rxp|max(v, Av}~ Using 
the lormubition of Sec. 11, one obtains -1, and 

(I -«A)r''-2 lor Uic average number ol 
optiral phonmi eniiSMnos in the respeciivr rases 
of W = 7 and | These results agree with that of 
the Simple algebraic appmarh.* Kor smaller R's 
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of ilii‘ U •»;, luirruTicnl 

rv’ultulKMis in a ut<lf-ci way ol flio ii t{-f;ra) 

«'qu.tlions in Sor It yiold llti' baiiii' fi lluii of Iho 
at^rbrau' .snlulion." y or Ihr iinalyf ir cunfit matiun 
fil llii> ociuivali'tirp ol (heso iwu (lilU'ronl upprcurh. 
Oh, a conltfciui (' It, |K(>n ♦ I ) - oilr"*’ |<>r H I, 
{hi I I), nti’dij to bi> proven, the A(v) is (lo|in*‘il as 
//(v)i h(v) r’’'*’ and IS f/(v) (or *-v The 
di'finitirin nl h'{>n) can iie (oimJ in lb f. B 'I liat 
lh.:i conjt dure is true Im in - 1 and 2 has 
cl.rt'keiJ 

Casetii) Iwn-dircriirin BoatiiTini; with iTjual 
piolxbilities, 

7'(v,,y,)= \ t»p(- !> ,- V, l)»: fxp|- (V, ♦>•,)!. (17) 

A'(y)* <’‘*®cosh(y) in thin case, Shown in Fig. 4 are 
the transition probabilities (or various Vt > Note that 
an electron, alter having a colbsion at a given lo- 
cation, may move against the external electric 
force (backward scattering' ind have another col- 
lision somewhere upstream in die electric field, 
further delaying the time at which it will reach the 
ionization threshold. Consequently, the mean ion- 
ization distance for this case is soiTifovh.it longer 
tlian for case (i) above. In Fig. 5, Ihe ionization 
events per mean free path of these two 

cases (or /f = 0.05 are compared in a scmilog 
scale. The one-half chancr to scatter baclward 
docs change substantially lor low electric 
fields (targe £,/eS\). 



iV'i. -1. fi'i" i\io-.liri(Utin i-rattenn,: willt 

rqu.i I !'• cl> i': 11110“. . II* n' A ,',/fr) I''. 





KIG i. I'.nization evcnli- per menir (lef path for vari- 
ous tvp.'s of scriltiTirv; h^; 0.0.1. 


Cat c (111): general scattering, 
nv,. v.Wv., ^P(0)d{f **| ‘*z.»i“»' 


(16) 

The functions S,(v,,.v,,B) and S,(y,,v,,B), as can be 
found in the Appendix, are the arc lengths of the 
electron from v, to the points on the 'rajectory at 
which the height is v,. h ii> the projectile angle 
and /’(6) the corresponding prolxibility |/'(f<)= const 
for isotropic scattermi;] Trivially, 

Af(y,)- (i9) 

There are three important observations con- 
cerning Eq. (16). The first is that the angular in- 
tegration appears as a multiplicative factor where- 
as in Daraff's formulation* it shows up in the ex- 
ponent. The second is about the behavior of 
T (, Vj , v,)at some seemingly singular imints. For 
example, at co.s'e<= v,..v,, both terms on the right- 
hand side nf Kq. (18) arc divergent. They sum up 
to a (mile value becau.se 5, - S vv hen cns' iii v.. 'v,. 
Another example is when cns(»= O m that the de- 
nr.rniiiaier in ft»e Irpantliiiiic fimclinn o( ,y,fi ^ . v , . M 
in Kq. (A2j hei onus zero. 7 he eos't' la el or iq 
iron* ol thr. |n|.,.irilhmir (uoriion eli mm.ites l*iis 
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t 


iJIt 


Jti 

h.iv.or I »»• ivl n.h.iviir n! T\v, <,) 

IS Mtfi «'( |'h)Mi aHv rl* ar Ih'Iu 0»*' drlimiinn oi 
7 o . . V , t ilM’JJ i:i vrr a bf.iMi rinj.' ,»t \ , , thi* pt'cUi 
hiliU fi! liK' III At Si ' lUft HI,' '*iltiii( a !iri|> il\' ■ :tl t, 
j.Kii'uIrt lx- liitili* tx'raust' Ihf .."fiiiiiil o( Ihc Irajt't- 
lui V .Ilf Hit III*' 'J.I-. Ml ip si. rinks 'i. /.fiu ul Uvi.'-i 
.1 (a-- 1 as «/\ (1 -M Uii .Ul's .i|)(ii‘-)iiiiial«' thprPN- 

n I'ti i ili'i" I (.1.1.1111 .1 .siiii’uL.1 ity at V, ‘ \ i. 

'Ill iliir'l rvtisnrv.it inn iiiviilvi's till normaliratirm 
( f 7'h V ; |(,.ul’' n (I 111 .lul/iiiia'H .tlly f.!ilii.ftcs 

!i.»' ( 1 1 | i; t •■i.H'iil fil fr>rn.ali/ i' iiiii it ib 

l'ii II ul.itfii II, a iM'I'i |)i (ilf.ii'ilistir wav. 

An r ,'ii v,.!i';v va> of i .ir r uii|', 'lUt I*., talcula- 
tidiis .iitl.ratcj If. Sf‘v. II IS to actiullj «" inputc the 
pr'’lulHiity of III ‘iptiral |ii.('iii'i) I'liiissin a for 
Ilf r 0 , 1 .2 . Oncti ,s and f arc knr.-.vn, these 

prolvtl'ililici. arc jusl siriplc iiitc|',rals. The irr.- 
p.ii l.mt .idvantaiti of thi.s is Ic t’hcck and see if the 
sum ol all protiatiili'ies equals uiuly, a very etrin* 
nmt constraint on ihc arcuracy of both the most 
lui.rtanii lUa! lum tin:. itself and all the nu- 

mcriral prociAlun s ifivnivci. In oilier wnids, 
once 7 (v^ , V ,) l^ ; cl, a probability diM rilmljnn is 
olitamcd which t.is to U Miinircd up to unity to be 
correct, This ri“,ui rciiicnt l'„i;- licin obeyed to 
gu.i r.i nil e the fr- t iciness r f • h.o ren, putt’d i esulLs. 

In actual ( alnJatioiis ,i i ti ,idy si, ate will be 
reached, as correrMy observed by ll.irall,* in 
which (he ratio of adjacent prohabililies of emit' 



VTit •• 1'rith.iHlitv flisti iliul'iiiis of two itiffiTrnt Beal' 
trrinc f :ii'lvr.i for Pime P aiv! V '1 ne sun of eaci. 
pros.'bi >dv di .sir ibt. lion cqu.-il.*: uniiv 


liitil III ibd If • I optical tiliniv'is appro.ic'ics a ceii- 
MnM, Thi constant, called >, ib strictly con- 
I relied Ir, 7(\ . 1 ,) Itself and ri.i, l«e computed hy 
usilit either the eigenvalue nieihod or iterative 
procedures. Using this constant . the protiability 
.it the o.iset r.l (he ste.uly M.tte, and (be sum of the 
prolj.ibitil.es l eforc rcachm,: ib" ste.idy sia't, one 
c.tn compuie the total proliability wluth should be 
unity if everything is right 

In I ig. (i we I liou the prolutulity distributioii;. ol 
i .'use (ill and case (ml (nr /( = b.o;) and A = ij. I lvi 
steady Mate ii. detei miiufl uhen the fraeluin.tl 
i baru:e o( die proli.dnlity ratio y •t' than lO’’' 

(in Sonic cases SO 

The result of the present work for the isotropic- 
•scattcring case is presetited m Fig. 5 for rt = 0.()!>. 
Also shown in the figure is a curve of Uarail for 
ident.cal a.B sumptions. The dillerenec is about a 
factor of 2 at X = 10 where, on die avcr.age, an 
ionization event occurs every thousand mean free 
paths (slow cluirge tnuUiplicalinn). The discrep- 
ancy gradu-tlls disappears as X is decreased. 
Comparison o( HaridCs work and the present ap- 
proach (or R^O.Oi sho'J S excellent agreement (or 
X (lom 2 to 14: the mean ionization distanre in 
this rase ts generally Ic.ss than HW mean (rec 
paths ((.iKt cliarge muUtplication). It is noted that 
in the sinw-niultiplication cases the steady states 
always possess probability ratios very close to 
I, i.c., Since the equation for the mean loni. 



no ~ Ncnilicr ol | .nri. jir >l•.^r^•^l ncr mean Inc path 
as ,a funtlion of h, for v.vri.ju.n vbIucs of P ‘-A»7F,. 
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y.Hion il ••uiiHi inviilvfs tfrtits prr^,uilin«.al to 
t (1 - >1 di)il I (I ->*•’, Ihf »e',utt i;» \«ry sensi- 
(ivp In 1 . ) u'lfn > i*< nrar 1. The y l)V itu' pri's- 
I til (aliulatini) lor /( ■ U.O'i atnl A ■■ 10, (or 
iti 0 '.'90 ft I ai)it till' sum ol all pi oUiIiiIiIil'H is a 
hatisiacini > 1.!'0|2. 

■J'l'i’ n'sulis n| Die laU'iilalmns of ar»' r'.'iwn 
III Fit’,. V. T.iblf I pivi's till* cnmpuU'd v.ilui's TIu' 
luiwi'i'ii llii-so rtssulis and H.iraO 's vi m k 
1 s I’.i oi-i .1 II'. K'llrrlor l.iri;i' ,V^,, (la.st u'uUiplKU- 
lini.) iti.ui I'lr siii.ill .V^i, (slo'Ji' mull iplicat u.o). 

A.s di'Siril* <1 at ivi‘, ili.s may I*' iluf In llii- di((.T- 

I III n qairi.'.'i.i I'l.s in ip.'iliii;; ) 's nr tlif iiiitd binp.u- 
larity in Il.ir.iff’s nxprt'.s.sion of llm approximati' 
transition prntialii lily density T{h^,h,). Tli'‘ over- 
all atireemuil i.s considered I'.rnfl in vipwintt the 
very diflcrcnt a|)proai'ln‘S used. 

Shockley’s for the high- field limit concep- 
tuiilly corresponds to nur h -0, In his simple 
model, the high-field Imiil corresponds to the 
rase of very tint clerironb viliosc eiirrpies are 
always grialer than H,. In car model, those hot 
electrons cannot produce any optical plionnii, thus 

II - 0. The argument for olitaininp, Inw-field lirmi 
m his mnde! is not apjilic.ilde to the speriat situa- 
tion of »,>= 0. Lven if we choose X, X‘l, the 
formula for .V,,,, is iliir. simple model reads 

c‘*/rt.Y, which does not corre.spond to any 
of our curves. 


tv CONCLUSION 

A theory is developed for the calculation of the 
mean ionization distance in solids. The founda 
tinn of this formalism is the same as that of the 
Doltzmann equation, namely, concepts of prolu* 
bility. This shortcut approach clearly relates the 
phyjics to the equations and has versatility in ap- 
plications. A simple .'ipplicaiion to the isotropic- 
scattering case produces results that generally 
agree with that obtained liy nunicncally solving 
the Doltaniar.n equation. 
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Am si>i\ 

Thi.s ap;x>ndii uiv*'S thi* I'mihuU lor Ihr l('ni;th 
o( the ^•U•cl^nn I ra}«-L'tnrv uivl>‘r tlio aclinii n( a 
i-i/itKiaiit ••klornal ck'tlnr fon i- Tlii' rcilinK 
u'Ih'It till' kiitrtic cnf’ri’v I'l Ih'- cU'flrrtn ts 7.cto ih 
as i •■■0. Tl>e cr)orclin.kti i iis iiii-asuied 

. 1 


downward ihe dirrrltnn of iho forcr) and 

lliiTcloM* rSr Tl>v eUriron at i, is pro- 

Jci'IihI at anpleb t $ {6 positivr) with rreiwci to thr 
horizon, as .shown in Fip B, Py t'lastiral nuch- 
anirb, the traji'ctni^ Icnitths from thr pnini a, to 
till* inierMTiimiK with a iMrizontal line ai hnttht 
Zj arc |;ivca by 


Si(z,. i() a 


I 1 , / I « sinf< \ 

,|Sinf(«cof:'olii ; \rrT—,~j ,ir: 

*,r'‘ ♦ (i,/a, - I- / 


,M<„ . , , W . < , [s,n 1. . . l„ ) 




(Al) 


(A2) 


If Ihcsr IcntMlis and tlic ^ coordinate arc (.rated 
by * and v ie defined as r/ A, ttic diniensinnlcbs 


I — — — 

Irajeclory lengths .S,(v,,v,.P) and S,(y,.\,,C) are 
obtained. They are used in Eq. |1B). 
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Dy treating the emivvinn of optical phonon* as a Markov process, a simple analytic method it 
developed for calculating the electronic ionization rale per unit length Tor dielectrics. Theeffecit 
of scattering from acoustic and optical phonons are neglected. The treatment obtains universal 
functions in recursive form, the theory depending on only two dimensionless energy ratios. A 
comparison of the present work with other numerical approaches indicates that the effect of 
scattering becomes important only when the electric potential energy drop in a mean free path for 
optical -phonon emission is less than about 23% of the ionization potential. A comparison with 
Monte Carlo results is also given for Teflon. 

PACS numbers: 77.30. -f p, 32.80. — s 


I. INTROOUCTIOM 

Interest in calculating ionization coefficients in semi- 
conductors and dielectrics dates back about three decades. 
McKay' derived a curve of a, versus applied electric field 
from experimental data. WollT developed a method for 
computing a, by expanding the electron distribution func- 
tion in terms of Legendre polynomials, keeping the first two 
terms in solving the Boltzmann transport equation in the 
steady state. This approximation is ii'stified at high electric 
field strengths. Shockley^ conjectured that the steady-state 
electron distribution function must have a spike in the direc- 
tion of the electric field By neglecting all other electrons not 
in this spike, he wu able to obtain an expreuion for such 
that the logarithm of a, was inversely proportional to the 
electric field. This approach is reasonable for low electric 
fields. Baraff,* in an effort to urify the work of Wolfl' and 
Shockley,* numerically solved the Boltzmann equation. He 
obtained an important plot of the quantity a, A versus the 
quantity E, /ttk for various values of the pammeter . 
In these expressions, the symbol A is the optical mean free 
path, E, is the energy gap between the valence and conduc- 
tion bands, Au is the energy of the optical phonon, and t is 
the external electric field. Further work by Baraff * concen- 
trated on the high-field case. A spherical harmonics expan- 
sion was introduced and a truncation procedure utilized to 
obtain a closed system of equations. A consideration of tem- 
perature effects was later given by Okuto and Crowell,* and 
Crowell and Sze.* An alternate approach, involving the con- 
cept of Markov processes, was studied by Lin.* The results of 
this approach agree with those of Baraff* for the quantity a, 
A at small values of E,/cfA, but disagree by a factor of two 
for larger values of E,/eeA. Monte Carlo studies ofhot elec- 
tron distributions in thin insulating films were reported by 
Baidyaroy el al.^ From these investigations an interesting 
scaling law for the average electron energy of the steady- 
state distribution was numerically obtained. 

To date, the calculations of ionization coefficients have 
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centered around the solution of the Boltzmann equation and 
concomitant electron distribution function. Though the cf- 
feett of scattering are properly included in this approach, the 
numerical methods required to solve the Boltzmann equa- 
tion may hinder one from gaining a direct insight into the 
problem. I n this paper, the problem will be approached from 
a different point of view. 

Just as the Boltzmann equation is derived by the con- 
cept of probability (via crou sections), so it will be assumed 
that the motion of an electron is a stochastic process. The 
two approaches are connected by the common concept of the 
mean free path for the underlying interactions. The crou 
section multiplied by the number density of scatterers gives 
the inverse mean free path of transport theory. For a stochas- 
tic proceu, the mean free path is defined in the usual fashion 
as the mean distance of free trajectories. The probability dis- 
tribution of free paths is aaeuined to be of exponential form 
u described in Sw. II. 

To simplify the problem, the effect of scattering is ne- 
glected, an assumption which is justified for large electric 
fields. This simplification, however, does not eliminate the 
essential aspect of electron multiplication in the phenomena 
of dielectric breakdown. By comparing the present results 
with those of Baraff*** (in which scattering is included), one 
can asseu the importance of scattering at various strengths 
of the external electric field. 

II. THEORY 

If A is the constant mean free path between collisions, 
the probability of no collision occuring for an electron trav- 
eling a distance x is exp( — x/A ) and the probability of one 
collision in dx is dx/A. The mean distance x between colli- 
sions is given by the expression 

jf= jT*xe-*^Vx//l->l, (1) 

thus providing a consistent definition of 

In the following, we will treat the collisions as stochas- 
tic events and calculate the ionization coefficient (for the 
case that the scattering effect is neglected) by regarding the 
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edtxl of the t'oIlKioHN 4% ihr evrtil% o( a Maikov prr^-m. 

Thu« an elecirnn i« rrleaMd with lero energy, gaint en> 
ergy from (hecnn«lunt electric Held t, and Ivnet an amount of 
energy Aru by einiuing an optical phonon of energy Aru aAer 
tranaveraingaconMant mean-free-path/i. The question to be 
addretaed analytically is, “What is the mean distance re> 
quired for the electron to ' :ach the ionization energy £,?” 

Shockley's* conjecture about the aforementioned spike 
is based on the Markov concept in a trivial way, as follows. 
The number of electrons, per unit volume, which have sur> 
vived the transport to energy E, without collisions, is pro- 
portional to expi — £. /egA. |, where the quantity £, /eg is the 
distance traveled by the electron to reach energy £,. These 
are the electrons which can "ionize" electrons in the valence 
band, i.e., promote them to the conduction band. The num- 
ber of ionizations in a unit length is therefore proportional to 
exp( — E,/egA ): 

a, « c esp( - E,/ggX ). (2) 

The results of the present analysis, on the other hand, 
indicate a more complicated dependence on [Et/egA ). 

To simplify the analysis, choose the ionization energy 
E, such that the relation £, *■ (n -f l)Aai holds, and let the 
quantity AZ, denote the mean distance that an electron will 
need to travel in order to reach the energy E , , given that the 
electron was released at the origin with an energy £ — Mu. 

There are three basic quantities involved in this prob- 
lem: ( 1 ) the quantum energy of the optical phonon Aiu, (2) the 
ionization energy £,, and (3) the electric potential energy 
drop in a mean-free-path eeA. Two independent dimension- 
less energy ratios can be constructed: 

R mftu/E,gnd X g»Et/egA. 

An important distance parameter of relevance to the prob- 
lem is the distance that an electron must travel in order to 
gain from the electric field an increment of energy equal to 
that of an optical phonon. This distance is given by the ex- 
pression ^/ee. Moreover, by choosing the mean-free-path 
A as the unit of length, as is done in the following, this dis- 
tance may be expressed as a dimensionless ratio RX, also 
denoted by AZ. The probability that no collision occurs for 
an electron traveling a distance dZ is exp( — RX ). The nota- 
tion avexpt - RX ) and^ ^RX is used in the following de- 
velopement. 

Figure 1 illustrates the evolution of the state of an elec- 
tron beginning with energy Mu at the origin as it travels a 
distance AZ in the applied field. If no phonon is emitted (an 
event with probability a), it will move to the state (/ -f 1). 
Emitting only one phonon, with probability afi, it will return 
to the state /. To end up in the state J\j the electron must 

emit (/- /-*- I) phonons. The probability for this event is 
aP{IJ\, v,herc P[IJ]:^' ^ * '/(/ — y 4 1)!. The probabilities 
for niultif )e emission processes are calculated via simple in- 
tegrations using the underlying probability assumption stat- 
ed al the beginning of this section. Finally, the electron has 
the icri.aining probability to emit / phonons and drop down 
to sia le I The probability for this event is not aP\l,\) because 
in the elastic region op:ical-phonon emission is forbidden 
The distance factor Z;|/<n ^ I ) is connected by a sequence 
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of appropriate probabilities to the distance factors + AZ 

where k takes on the values k ^ [I + 1), /,—,l. Thus Zg is 
related to Z, by the equation Zg ■■ Z, -h AZ. The general 
relation is given by the equation 

Z, - a{Z,^ ,+AZ)A‘a2P{IJ) (Z, -h AZ) 

/-* 

+ [l-a-aJ;^/»(/y)](Z,-|-dZ). (3) 

Because the goal is to calculate Zg, it is natural to eliminaw 
the intermediate quantities Z, in the above relations by in- 
troducing a series of kemals K (/ ) which satisfy the equation 

Z,^Z^-K[l\AZ. (4) 

It is apparent that the fint of these, £ ( 1 ) has a value of unity, 
£(1) » 1. By substituting the defining relation of Eq. (4) into 
Eq. (3), a recursive equation for the kernel is obtained. This 
equation is 

a£(/-»- D- I -h£(/)-o ^ P[lj\K\D 

l_a-o ^ F(/^)£(l). (5) 

y » 3 P 

This equation, together with the initialization condition 
£(1) — 1, determines all of the kemals £ (/ ). By impection. 
the first few are given by 

£(2) = o-'-»- 1; 

Jk(3) = a-*-Hl 4- 1. 

It should be noted that Eq. (5) depends only on the two di- 
mensionless energy ratios R and X because a and P[!J\ 
functions of R and X. 

Even though the kemals £ (/ ) are known by solving Eq- 
(5), the various values of Z,|/ = l,2—.n -f 1) cannot immedi- 
ately be obtained from Eq. (4) because the quantity Z„ is still 
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FlU 2. loniMtion coctncient m • function of rlctiric field for vinous K Y 
« iFe number of cicciron-holc piirt produced in « mein free pith; X ii the 
mcrf> fOK* •oniuiion energy to electric po«entiil-ener|y drop in i mein 
fire pith 

unknown. It is powible, however, to compute the special 
distance Z, , by other means, as is described below. The 
quantity Zg may then be obtained by an application of Eq. 
|4). For electron energies greater than the ionization thresh- 
old K (n -f I|^, the electron of interest can raise a va- 
lence-band electron into the conduction band (collision ion- 
ization). This occurs in travelling a mean free ionization path 
which is assumed to be much smaller than <1. With this 
assumption, the relation Z, , *A,/A is immediately ob- 
tained. Setting / B + I in Eq. (4) then determines Zg as a 
function of R and X. ' 

With the mean-free-path chosen as the unit of length, 
the ionization coefficient (denoted by F ) is given by 
Y m I/Zg. Thus, the number of ionization events F in a 
mean-free-path is a function of the two energy ratios 

R kAu/Z, 
and 

XmmE,/e€l. 

III. DISCUSSIONS AND CONCLUSION 

Calculations have been carried out for various values of 
the parameter R These results are shown in Fig. 2 

fortheconditionthat^, = /I. The cun c for/? - 0(n <»)is 

obtained by special considerations In this limit (zero opti- 
cal-phonon energy), the electron docs not lose any energy via 
the emission process so that the relation — E,/cf + A, 


holds. 1 his reialkw implies that F is given by the eapression 
F -• 1/(1 -f Z) for 4, * 4. |6| 

llic curves of Fig 2 as generated by the present aiiJ) sis 
do not have the simple behavior implied by Eq. (2) (a straight 
line on a semilog plot), except for the case when R a 0.5. 
This complicated behavior is expected because there are 
terms of the type exp(m/?4r ) (m « l,2",« Min the equa- 
tion determining A|n -f I). It is only for the condition lliat 
m B /I -I I lhatexp(/n/?Jf ) >- exp(.V|, the behavior of Eq. (2). 

Comparison with BaraflTs* and Lin's* results, in which 
an isotropic scattering effect is included, shows that the pre- 
sent results agree with theirs only for small values of 
X{XS 4). In this regime, the electric field is so high that the 
direction of the electron after an isotropic scattering reverts 
almost immediately to the direction of the applied Held. On 
the other hand, for large X (small electric field), the scatter- 
ing effect should be very important. Indeed, the present re- 
sults deviate significantly from theirs under these condi- 
tions. The difference is shown in Fig. 3 for a value of 
R - 0.02. 

The curves in Fig. 2 are the same in general form as 
those from more realistic calculations* *; thus increasing the 
ratio of the optical-phonon energy (Ara) to the ionization 
threshold energy (£, ) rapidly decreases the number of ioniza- 
tions per mean free path. Though the isotropic scattering 
effect must be included in describing ihi: artual physical situ- 
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FIO 4 lonixaiion length t> a function of uicn'al electric Ack* for CF, 

ation, the etsential feitureft of electron multiplication in die- 
lectrics is contained in. and described by, the present simple 
approach. 

The above approach has been applied to a crude model 
of Teflon. The opncal-phonon mean free path is chosen to be 
^ « 26 the ionization energy is taken to be 

6.S eV,'" while the optical-phonon energy is given the 
value 0. 1 1 eV. The calculated ionization length for 
this model is plotted in Fig. 4 as a function of the applied 
electric field In the region above 10“' V/m (the high-field 
region), the reliability of the present model is greatest. Also 
shown in the same figure is the result of a Monte Carlo calcu- 
lation of the same quantity." The Monte Carlo model in- 
cludes the effects ol'an energy-dependent mean free path" 
(/’ = 26 A/( I - ttu)/E )"’) and elastic scattering. The two 
approaches give results which are in reasonable agreement 
(within a factor of 2) in the high-field region above lO"* V/m 
(see Fig 4). 


In summary, a method has been derived for solving ni,' 
problem of two competing processes in which the cfact of 
one of the processes is to produce gain in a continuous fasl^' 
ion, while the effect of the other is to dissipate that gain ia a 
stochastic way. An application of this method to the study of ' 
a simplified model of electron multiplication in diclectria ■ 
yields the essential features of the electron ionization codk i 
cient. This coefficient is a central feature of all theorim of ' 
dicicctnc breakdown which caplicilly invoke impact ionita- ' 
tion. The relative merits of these various theories, at well aa 
theories which do not capiicitly utilize impact ioniution (for 
ctamplc, the theory of Frdhiich'’) is best ascertained from 
the literature on the subject.'* 

Other applicatiofu of the method include the analysis of 
nose currents in device physics; this may produce us^ul 
information about the process occuring in the device in the 
prebreakdown stage. The steady-state energy distribution of 
hot electrons in thin insulating Alms subjected to high-clrr- 
trK fields could also be obtained by the present laet^. 
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lltKOKKlll..! !‘k l|’; RUt '»i srNlVViK 

l.\ ;U)I iP sili.i.i i I Kit;.* 


h , 1. . iSi'i r !> > V , Vi . I I !!•' , .trii! ^ . I . I /i . 

HutTb AsMi.o l ltOH, lUi . 

i’OKl (Uflii- Hex J '>•»') 

Kuulon, Vltk.lKl.i. J.’O'iU 

INTRODliCriON 

A "lodrl h.v« prr.’lotiHlv hit n ilos. itli.l O.-) whlih 
pr 'viilfb for the sill ti'ii at stent .nimerti.il doatrlptl n 
ot tin- spjtl.tlly ninl tompor.iMv ioIumoiu i v. liii Icn of 
ill. tti'f. .IV. 1 1 .Indies tn sol Id if lo 1 t r i • . V.irlcuia 

tt. It'll o. .'t the prop.i.vit iiu: !U>,.itl'e tip <>t I )u i v.i I .iin he 
i.ive hi on pteert..! li), .in.! .jiiilitt i^i Wi 1 1 olteleil 
i!i..iit the pissihli I'^por t .iiKi' I't ri'i'uirl.i^ jr 'ves'.es .it 
t'l. positive tip le.iillnn ti> i.itt.’de illre. ti.l Mre.imer !.'i 
.iinl subse'pient bre.iktlown ( 2 ). .\n extensive set of 

niinei tcal e.ileul.it Iona fi.is now been peril r.sed within this 
m.i.lel i riimework. The purpose of the cjiKiil.it Ions w.is to 
i-.pl.ie the depet'Ji n. e of Import. int i .i 1 on 1 1 1 f oivi 1 fe.ilnre 
on the Input par.imeieis tif »fie model. A prollmln.irv de- 
s riptli'n .'f the results of this p.irameter study Ki p.l'i-n 
1 1 this note. The research 1 •; p.irt of an onKoitiR proitrar. 
to urderatand the vharRlnR and discharKlnR properties of 
elect r.’ii beam Irractl.itcd solid dielectrics. 

M.iIlklAl TAh.VM! n.kS A'.T' IKAN.sPOKi’ UiM iUil.MS 

I hi tr.insp.irt c.-i t t Ic lent s t" I’C .i.M-d .is p.irarict it.s 
• .r t h. p..!i I oscop ii di's.rlptlon of f fie .u.il mchc rust 
tir.-.t hi .ibt. lined. In t lu- reported ippro.iil-. (1), tliese 

piri.melcrs iie of't.ilncd bv solvlnc the Kolt’m.iiin tr. ins- 
port eijiiat leii for the siiiRle electron distribution fuiic- 
tl.m. The Holuli.'ii la obtained f>v .i s Inple-si .itter inR 
Monte t.irlo appri'.ith as lmplement.il tn tin' computer code 


• Work supported bv tlie Natlon.il Aeron.iut Irs and Spi.'o 
Ad" 1 n 1 .St r.it Ion , fowls Roscarc!i ('eiitei under Contract 
.‘.un'bor NAS '-if 1; f'.K 
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PROEM 

QUMJTV 


I'M 


s. '''i (i*. III! r><,uiii‘ti ii (Mil (i.it.i .irt- •! ri.'lfl ■'! t lio 

s> .It t fi liip (>i >ii rs-u-H (.1 lutcrlil rk«1i t ( anJ I lit' 
i, l<lli <l iloiirK f ii lii It rcitpl li . 

} t till. I'lirpKHf «.i| ililu |i.ir.irii ti'i htuiiVt u vvrv 
Kl-'iilr mitirt^l m.'vltl w/im tliokcn. Thrrp Hint icr Inc, 
rr wfii' ciiiiii Idcrc'd ; rc at tt^r Iiik t ron .irouitt 1 l pli>>iitiini ; 
<>nib' K'li ol opt > 1.^1 p'umims; and t. hc ilnt ion of iiucrliciJ 
ti .-insit tons I)..)! I Ik loll lonir.il ton; . Thi' mji. free p«th 
I li r .iroiistii' plimion ncdl ti-r iiip v»u tokru to be con- 
stant. riiif •pltcnl phonon spcctiuai w«» ropr«R«nttd with 
.1 sinitK- optli iil plion<'n of energy It . The energy de- 
pt i.di-n. i* of the mc.ir. free p.ith for opt leal phonon einle- 
Hlon in tiMt Riven by SeitP (3), while the ttnlc iv ect 
with j viitldhle pnrdneicr Interbai.d tranaitionn <irc 

dvbcritod with j thrcKhold for r luet r.m-hole creation 
tcgftlier with d connt.mt me.ii) free path above the 
thieshnld. Scattering from acuunt Ic phonona waa taken to 
be Isotropic, whl!e both inelastic proccnaob were aasumud 
to occur without angular deflection. 

A material model la therefore dearrlbed with the 
following act of parametera: 

\.r *o‘ !•'" • 

whore i( lb the eU' tronic charge, m* the effective luab 
of the electron, and t the relative dielectric conatant 
of the material. The Fcouilnlng problem parameter ia the 
impresacd electric field The aet of qnantltlea 
im*,A .iCf) vai> chosen ea the fundam/fital dimenalonel aet, 
and tfie Boltzm.in equation waa brought to non-dimcnaional 
form. The reaultlng equation ia totally Independent of 
r . The following five dlmenalonleaa parametera charac- 
terize the reaulting equation: A. - a_/'> .; R •Hc/E.; X* 

the linear problem* la*lndependent of .'I. For this pareme- 
te‘r atudy, the dependence on the variable waa aupprea- 
acd. A value of X. was chosen suf f iclently*small that 
ionization occurs practically Inatantancously once the 
threshold is reached. This ia a good physical assumption', 
and is in accord with other rcaearchers (4,3,6). T|ie 
p.nramcter X is taken to be a measure of the electric 
field. I'p to scaling by the fundamental dimensional aet, 
d iwiterlal model is therefore specified by a poir of 
values 

Computations of the transport coefficients were 
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p«rfv»rni.J vlih viliif. if ‘ in the nii/,. f ron (I, 2'>J 
.ind v;ilucH of K In th*» rjti^c from I /• x 10"^, 2 x I0*l)» 

Both thi* difiuf'iin cowf fie lent .iiiJ the mobility were 
found to bi* rnlher Insensitive to llu vj’ues K nin! ■. in 
the high field limit (X in the nn,;.* of (2. 10 1 ) wliich ih 
tile parameter rttgine of concern for tiiib n.)io. Both U 
and Vi vary aa a function of the field atrength. In 
<1 iT.ena lonloMM unltHtO Increasta from 0.1 to 0.2 as X le 
dveraased from 30 to 2 • cor rcapund Ing to electron iMtat' 
ing for the higher field values. Similarly, u decreases 
from 2.2 to 0.6 over the same parameter range, exhibiting 
the anticipated trend toward saturation of the drift 
Velocity at high fields. The mean energy (effective 
temperature) shows a slightly greater variation with 
the parameters i'-q.K. indicating that the coefficient in 
the Einstein relation depends mure strongly on these 
parameters th.in v or D individually. The mean energy 
varlcb by .ibout a factor of five from the lowest fields 
to the hi|;hest, wi ih il : t.i !\ s for 1 .. v..liieh in th-. r.tti,;e 
i t X 10*2, 1,8 X 10“ n depend in;* on the values of ■ Q»h- 
The quantity most sensitive to the paranietcr 
values (and the field) is the ionization length, or its 
inverse i (the first Townsend coefficient). Shown in 
figure 1 are plots of m (in units of a func- 

tion of X for 10 different scattering sets. Olven the 
rough independence of u and D on the parameters and 
R, the quantity '« Is the quantity which characterizes a 
material model for the macroscopic description discussed 
below. The general features of the curves in figure 1 
may be understood from pre 'ous seni-analyt leal work (4, 
5.6). 


SrRtA,MER EVOI-UTION 

The nvicroscopic model of the scl f-co«is i stent evolu- 
tion of electron avalanche requires specification of the 
parameters p, D and a (1,2). A universal value of D 
was chosen, given Its weak dependence on both material 
model and field. A univeraal fit of the mobility aa a 
function of field strength was chosen because of its 
weak dependet c ! on material model. Finally, the ionlza-. 
tlon coefficient i was fit with the usual form 

• - *0 (- 

where 'i^ and are the fitting parameters which depend 
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I'll 'll .iml K. Mils fiiiutliii.il li'iri pi lu ;>li'i. <1 r.itlu'i 
pi'iiil I rpi i-si'iit .It li'ii Ilf tilt' ilitii I'l !U;iilt I. 

Till' lull till I I'liit 1 1 ioi|i tn li.tipl It lni till 
tliiiilts iipi.itliiiN wi'ii ill ti'i mini d .it> I ilIiiWN. A ptlnuiiv 
ili\tiin .iv.il nil hi ui*. ilrvi- 1 1 'prit In i Mi nit t.iiln l.isluin 
II. iii^ till' il I St I lltiil lt>n til It'll 1 7.it iiiiiN iilil . I lilt'd I I I’ln 
M,.'K . (1). Ilif ivnlulli'ii w.is Hli'pi'i'tl will n till' hi'lf'llild 
■ <f till* .'IV. I I .nil III- I'ht .1 llll-tl .1 .'.lllli' I'f SO' ii| till- .ippi ll*il 
fli-l'l. lilt' ri suiting; iIuiiki iI I i.l i ll.iil IiMI w..-. usi-d .it. 
lilt- ht.iilln^ print li-t till- m.ii rtiKi'i'p It' di-si i i pt 1 1 'li . 
liit.pi-t t tni t'f I tu-'u- il I !.t I lluit loiiH lnl^ sill wn t lu-n. to hv 
riitlti'ly I I'tiNlHtnnt with hlm|li* ainilytli.il 1 1 pi nm nta- 
I lull). ii| till ptlir.ny v.iid 1 si m t i-d . .ivn I .in. In- (7). 

Till- nwii I I'Si up li' I ijin. -ili ili-ricill IliK t l.f I'Vnl nl Irn 
i>! till' .i\ . 1 1 niclir r iv hi' I’l.-np.lil tr i I' .i lit. l.'i It s.s l.itmtl), 
I. Intli-iif. tli'p. ml i'll twr il I r.i I '■ 1 1 n 1 1 .s . p.it ,i: .ri ri t wh i li 
li.ivi- 111 1 11 i'liii).i-n t-' Ilf ■ • (lti) f) 111 . I , • (i) Ti- 

lin' I ll'.t p.i I .rni' 1 1 I ir.iv l>i i I'i'xpi .'-..si'd .ii- lilt* l.it It' i>l 
till nil-. Ill I'li-i.tii'n I 111 t ny ti- tin iiii'tpN ixti.iiti-il linn 
tin- f li-ld In .in It-ii I .Ml inn lii.i-tli. Tlif hf t I'lid is t lif 
j.itln t’f twn i-li’i'tilf fli'ltls; tli.it dn. .i sliip.lf ili'i'irnn 
Itu.iti-d one Inn 1 ..It Inn li-iip.tli illst.mi. .iiid tin- .imhliiit 
lli'ld. 1 arp,iT valiii-s nf nirropt'iid In dl(fui.inn 
hflnp • f 1 .It 1 vi- 1 V m.iri Imp.'i t .int lli.in drift, wlilli •r..illt-r 
v.i 1 III"; I'l't 1 i'spi>ii.l !n lilt- nppi'silf r I'li.l 1 1 Inn . Till' iju. 'ti- 
ll I , Is .1 mi'.isnti' nf tin- slit-nplli .-t tin- i-xti-rn.il 
111 Id ‘ill ti-rms nf till' I iindiiiii'iUa 1 tnnlnnh M'puiatinn 
field. For the nodel p.irumiti'is .ikh-'i l.it ed with t I pure 
I, f. . v.it U'k In the i.ingi- .>f (<i x I0'^, 1 x 10“' 1 while 
v.irles In the r.inpe |S x I 0 *h, 7 \ 10“'). 

Until p.iranietels show wldi* vaii.itlnn fni the nmili 1 
v.iliier. I hnwen fnr tlilr rtudv. Fnr the mn.lels .ind field 

levi-l« .ixHOclated with tlpnie 1, ,i nc.ittet pint nf tin- 
actual values of f .iiid ’ , w.ih peimr.ileil . On a Inp-lnp 
plot, the valueu rijspl.i\til .i rcm.irkah I e irciid. F.ssfn- 
tfally all the values lie within a narrnw hand of a 
airtight line on the plot. The ntralght line has a 
slope of rnughlv 1.'. cor resp.ind Inp to a conatani value 
nf the qu.intlty 3^*). Conirut.it Ion of for 

the vnrlniia model values allows that it varies hv less 
than one order of nvignltude (corresponding to the hand 
width) over this wide range of model values. It is 
poNslblc to examine the solution* of the equations ns 
a function of a single pnr.ineter (s.iy *j) with the re- 
n.ilning parameter determined h> thi- const. nicy of f..j. 


120 



ORIQINAL PAGE 18 
OF POOR QUAUTY 


)>4 

■ ■ 'j j',n .4t 1) tl»f iiu.nlier *»f com 

pul.il IiuIh 1 rc'l to cNpInrc the I’jraD.ctoi :>p.u'v. Thin 
i.k.iliiip l.iw in not i.tidt-rntuod at proHiii*.. It alioul'J be 
U'leJ, li.'Vi\.r, tli.it '.j Is the ll:uiij|.:il. , p.i i .i:. i-l er 
whlcli occur:. wlu :i Otit Ir.il iiig the Ktrcri|(th >'f llio sclf- 
I'liltl In the unJlni or'cl iv.il ui.i’Ik ol lulled t* uhovo (7). 
Ihe s iKn 1 f It aiice nf thlb colncldriice of paraneters i« 

•I 1 Bo m>t undet stood. 

lh.>ii£li stri'Uiily lndlcitiv«, the .ihove si.iliii^ dis- 
cussion is overs in.pl i f led as it only cuncernu imolf 
with the overall dimension.il .nspects of the problem. In 
fact, the transport coefficients .. and • have a rather 
strong dependence on the field, and vary significantly In 
s; .ice inJ tine during the iv.ilut 1 n of the stre.iruT. 

Ihe pnr.ir.oters determining this viilitlo'i add furllier 
d liuMis i o;' 1 • s- I'.iiliMt; to the pr.'hlem. 11. e pos..l)illty 
th.it the d« pendence issue i.ited with this coefficient 
v.iriillon in a single pri'Meir. mig'it h.ive a malar iTip.ict 
on I lu' ‘oim ol till' soluti'in le.id to an initi.il choice of 
o-niul.it ions h.i'ej prim.irily on distinctive function.il 
toims fir till' transport coef f ic ieiits . Referring to 
figure I, it was decided to compute the evolution lor 
ific ioni/.ition curves .iss.'ciated with model numbers 2, 9 
and 10, conput at ions being performed over the span of 
field i..ilues I;’ tlu- flgiir.. I'nfortunalely , this choice 
li.id t.' a scries oi comput.it ions In which the Jiniaision- 
l";,'. p.ir meter ’ , only v.irlid by t'wo {'rders of magnitude, 
.-j-U s lO-‘, A x‘10-»), ^ 110-2, 8 X 10-1) rather than 

ll’•••r tlie full r.ing»« implied by the previous discussion, 
f. plot. It ion of some of the smaller valui-s of awaits 
♦viture work. 

It folliws from the defining eijuat ions that the total 
1.1 diet of ele'.lT-ns in the stre.imer grows exponent ial ly 
;i.. "'p • ^ Indicates .in electron density 

welgfited .ivi't.ige of the lonl/.itlon r.ite. lln* f;ua“tlly 
• provides the simplest diagnostic on the distribution 
of electrons .inJ fields during the avalinche. Shown in 
f igiKe .’ .ire the results of eleven separate computations 
of the growth of the number of electrons as a function 
I'f I t. The solid line is a best fit to the set of all 
p. 'lilts In this composite tlgure. Individual data points 
.ire indicated to show the spread, Tlic figure auggests 
Very strongly that seemingly dlfletent compunt ions are 
in i.icl ve-y much the same when viewed in .ipproprlale 
d iru tit l>>n I lu s var l.ibl es . 
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Aiuitlit i I'lti.imott I vliuli I'tiviJis .1 tl i>l 

llu' i'voliit ii<ii th till* |>i.ik ut'Mi.Uv (k'l ihf nutxl- 

nium ttoli! i>nti tnci'i.ii ii t ) . !ii •! i' u<ii>. i i >r, 1 1 s;> uiiils, t la 
•ml It I'ciik thdi^i di'imlly ih ^Ivt-n by In* 

titl'd, loi f.ith i'( I be si'^'.itati t itn^tul .il liMm iit‘>lt-d abovt , 
llii»- v.iliii i.' ati.ilni'd tiK a s,.t ui at Ikd v.iIuc. I'liiiilly, 
llu- posit Inn ol the tonirold ot iu't;at ivi- chaigi- as a 
luntllon ol I tmi' rviy hi lollowotl In il Inu ns i«>it 1 « ns units. 
Thesf points also lie loiighly on thi- sdnit cuivi , with 
the same degiee ot scattei as In I ip.urr i. Thu resulting 
dtlft veil', ilv l;i helnp, .ift. t lei .it til al a tale of l.sxlO' 
tin illmenslonlesN units ^p>). 

The above diagnostles, topethir with del. tiled ex- 
amlnat ion I'l the evolution of the \ai ious computations 
veiv strongly suggest that the ev'>liitloii of negative- 
tip • tte.imeis tends to a similarily solution ovei a 
large range of problem parameter values. This type of 
behavior oiay be expected from systems exhibiting negative 
feedback. 


fl) P.L. Beers. V.W. Pine, II. T. Ilw.tng , and 

II. W. Bl.U'mbeig, 1, 1 . h . JR.XN:? . 

NS-26, 5127 (Decenrber’, 197d). 

(2) B.t. . Beers, H.C. Hwang, and H.W. Bloomberg, 

(1979). p. 212. 

(3) F. Scltr, PJh^s. Rev. 73, 549 (1948). 

(4) i;.A. Barnff. Plos. Rev. ]28. 2507 (1962), 

D'X"' n3." A2b ~~(1964y.' 

(5) D.L. Mil, Phjfs. Rev. ^0, 5238 (1979). 

(6) D.I.. I.ln and B.I.. Beers, J. Appl. Phys. . 

(to be publ Ished) . 

% 

(7) H. Rsuthcr, Fl eet ron Avlsnches and Breakdown 
In Oases . BuMerwurth, Washington, (196^). 
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Ficiure 1 

i’lot of tJiu inverse ionization 
length i (in units of versus 

tiio dimensionless fielti variiibles 
Ei(q^iel) -or ten different nuternl 
nxxJels (see te-xt for discussion of 
nuterial model ). 
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Figure 2 

Plot of the miflber of electrons in 
the streoner siiailation as a function 
of the dimensionless time Pt. Itie 
plot is a ocmposite of eleven distinct 
ociifxitations for different problem 
Iivametcrs. 
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ABSTRACT 

Tht primary purpost of this paptr It to 
tuggatt that low difftrcntlal voltaga, tawll tnergy 
ditchargat attoc1«tad with Internal burled charge may 
ba a dominant mechanism by wnich ttorad alactrottatic 
energy Is released from dielectrics on board orbiting 
spacecraft. The evidence from space at given by 
Stevens^ It noted. Laboratory experimental evidence 
of Freder1ckson< Is Invoked to demonstrate that d1s> 
charges occur under circumstances with no external 
potential drop. Previous calculations^'^ are then 
reviewed which Indicate that significant Internal elec- 
tric fields can exist In dielectrics charged with multi- 
ple-kilovolt elcrtron beams under conditions Involving 
little or no external potential drop. The Internal 
discharge mechanism of Neulenberg' Is '-ecalled, and 
new calculations suggesting that the space environment 
Is conducive to the formation of the conditions re- 
quired by this mechanism are presented. Experimental 
procedures for checking the suggestions of the paper 
are developed. 

BACKGROUND 

Early In the history of concern with the 
spacecraft charging problem, Meulenberg^ proposed a 
mechanism for break-down in electron irradiated dielec- 
trics which Is critically dependent of the character of 
the electric fields very near to the surface of the 
material. For a variety of reasons, subsequent experi- 
mental work did not concentrate on the conditions re- 
quired for this mechanism. Most work was directed to 
configurations resulting In large surface to substrate 
differential potentials. For these charging conditions, 
very severe and distinctive discharges were observ- 
ed.**' •8,9,10 conditions of the reoorted experi- 
ments (referred to as "Standard" herein) are those 
associated with a monoenergetic electron beam Incident 
on a sample with a grounded substrate. For sufficient- 
ly large beam energy and/or current, discharges were 
universally observed. The area scaling relationships 
that were expected and observed^ suggested extremely 
large amplitude, long pulse width discharges for the 
large areas typical of spacecraft. Reports^^ of the 
observation of these monster discharges under "standard" 
conditions lead to the realization that If this type 
of experiment were Indicative of space conditions, then 
the spacecraft charoino problem was extremt^ly severe. 
Direct measurements *0,11 of the electric and magnetic 
fields associated with these large discharges directly 
supported the Idea developed earlier that these large 
current discharges were associated with a space-charge 
"blow-off" current. 

Within the past year, the idea that the 
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discharges associated with these "standard" simulations 
may be significantly different from those which occur 
In space began to evolve. 

Recent experimental work by Frcderlckson^ 
has shown that most types of coamwn spacKraft dielec- 
trics undergo spontaneous internal breakdown when ir- 
radiated by h,gh-energy electrons (n 1 MeV). By In- 
ternal breakdown, we mean that the experimental con- 
figuration Is such that there are no external potential 
drops across the dielectric sample -- breakdown must be 
caused by field Internal to the dielectric. While 
these results are no particular surprise to workers In 
the area, and earlier worki? has shown similar effects 
In more complicated configurations, Frederickson's re- 
sults do establish the presence of the effect beyond 
doubt. The results also establish the fact that any 
program designed to ascertain breakdown criteria for 
spacecraft dielectrics must consider far r«re detail 
than has heretofore been assimmd. In particular, the 
measurement of external potential differences under a 
variety of charging conditions Is clearly sufficient 
to characterize the onset of breakdown. 

In his paper, Frederlckson^ explicitly 
notes that the generic features which characterize the 
accumulation of charge for NeV electron Irradiations 
are no different than those which characterize multiple 
keV electron Irradiations. The Implication of this 
statement is clear. Irradiations with multiplc-keV 
electron spectra characteristics of the spacecraft 
charging environment can be expected to lead to break- 
down even when external potential differences are small. 
Supporting evidence for this viewpolrt Is also becoming 
available from space. Stevens^ has pointed out that 
SCATHA data Indicate that breakdowns are occurring on- 
orbit. despite the fact that only small differential 
voltages are being observed. The appearance of small, 
relatively benign discharges, under the conditions of 
small differential potentials suggests that the "stan- 
dard" simulations are giving rise to qualitatively dis- 
similar discharges caused by a different mechanism. 

The internal breakdown observed by Frederickson^ and 
proposed by Meulenberg^ Is certainly a candidate mecha- 
nism. The conditions required for this mechanism In 
a space-like environment are explored below. 

CONCEPTUAL MODEL 

Previous papers3,4 on the subject of In- 
ternal charge distributions In dielectrics subject to 
spacecraft charging type environments have provided the 
computational frame work for understanding the dielec- 
tric conditions of Importance for internal breakdown. 
Indeed, It was explicitly noted in these publications 
that low-differential voltage Internal breakdowns were 
a real possibility. 

The non-charging beam Introduced In Refer- 
ence 3 Is a convenient device for understsnding the 
essential features of Internal field buildup. A 
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non-ch«ry bt««i for t i^arti'.ular mitertdl i. 4 ooji.. 
for whUh the surface emits precisely one electron for 
eecb electron incident on the > iterldl. A»^uninii thet 
the emitted electrons ere prluoi ily scconder1es«ri>in1t> 
ted froei a region very near tfe surface f < 00 X)], and 
that the primary electrons per.ctrjtc suOstantlally dicp- 
er than the secondary emission region, it It apparent 
that the prlemry deposition leads to a charge »epara- 
tion In the material. This region of sepsrated charge 
contains an Internal electric field (even though there 
Is no external potential). For conceptual purposes, 
assume all secondary charge is emitted from the precise 
geometric surface, and that all primary charge is de- 
posited at a single depth d. Let o be the conductivity 
In the region from the surface to a depth d (assumed 
for simplicity to be uniform). The Internal electric 
field Is given by the well-known expression: 

E • - (1 - exp (- f t)) ^ . 

where J is the beam current density, and c the permit- 
tivity of the med1ue>. 


tiavi'i'i chosen a surface neutral Izur , the 
basic experimental program Involves varying the energy 
and current density of the penetratinj beam until break- 
down occurs. The location of the breakdown threshold 
as a function of deposited charge and beam energy should 
provide an excellent empericjl criterion for Internal 
breakdown. Combining such experiments with conputa- 
tlons for the conf iguratton will permit the Identifica- 
tion of Internal conditions at breakdown, as has al- 
ready been shown. 1 It might be speculated that break- 
down will occur when Internal fields reach levels of a 
few times lO' V/m - the breakdown levels found by Fred- 
erickson.l 

The point of this section Is very simple. 
Theia exist free floating front surface chargina con- 
figurations such that no potential drop Is developed 
between the surface and substrate for which significant 
Internal flalds art developed. These large fields are 
near tha front surface. If larca enough, they can Itad 
to breakdown In this region. 

RESULTS 


Essentially all qualitative features of the 
Internal charge and field buildup may be ascertained 
from this expression. The emgnltude of the Internal 
field Is limited only by the conductivity of the Irrad- 
iated region. Writing the conductivity as Oq ^ Oo • 
where Oq Is the ambient conductivity In the unlrradla- 
ted utirlal, and Is radiation-induced conductivity 
(RIC), the various regions may be explored depending 
on whether Oq or Op Is dominant. We will only considtr 
the circumstances "in which op Is dominant, because 
most materials of practical Interest exhibit this domi- 
nance at levels of current characteristic of exoatmos- 
pherlc environment. 

If Op Is linear In dose rate, the Internal 
field will saturate at a maximum material dependent 
value, which may or may not be above the breakdown 
strength for the given material. If o» 1$ sub-linear 
In dose (as It almost surely Is for high enough dose 
rates and/or doses), then the field *ay be driven arbi- 
trarily high by increasing the driving current. Break- 
down Is guaranteed for sufficiently high current level. 

Examination of the range-energy relationship 
for electrons shows that the mean energy loss over the 
range of the electron decreases with Increasing primary 
electron energy. It follows that the mean volumetric 
dose throughout the deposition region decreases for In- 
creasing electron ener^. Thus, for a given material, 
and given current level , an Increase In the energy of 
the primary electrons will lead to an Increase In the 
Internal field. 

Both the above statements are easily quanti- 
fied within a computational model of the type prevli^js- 
ly given. 2. 3 We prefer to explore the possibilities for 
achieving the above conditions. A single monoenergetic 
beam can be non-charging as noted In Reference 3. For 
a given material, only a single Incident energy can 
satisfy the conditions, so the possibility for varying 
the second parameter (depth) Is not available In this 
configuration. The simplest way to achieve non-charg- 
ing with a given monoenergetic electron beam is to pro- 
vide a second source of Irradiation which produces a 
copious supply of secondaries. Possibilities for the 
second source are: low energy electron beam. Ion berrii, 

photons. The thiro choice Is less desirable for mater- 
ials which are photo conductive for the wavelengths of 
Interest, as it Introduces yet another variable Into 
the problem. Once the basic effect Is quantified, of 
course, it will be Important to Include sunlight for 
full environmental simulation. 


A number of calculations have been perform- 
ed to compare the Internal charge and field conditions 
In a typical dielectric due to a monoenergetic electron 
beam to those due to a space-like electron distribution. 
The space-like electron environment which was chosen 
for simulation was an Isotropic Maxwellian with a temp- 
erature of 11.5 keV. This environment was taken to re- 
present the hot electron portion of a double Maxwellian 
fit to the electron distribution function durina a 
severe substorm. 13 The electron transport model for 
determining conditions Interior to the dielectric has 
previously been described. 3 The model Is one-dimension- 
al . To specify the charging, the Incident electron 
distribution function Is required as a function of the 
surface voltage. For true space-like conditions, this 
specification reouirts t!ie solution of the coupled 
Polsson-Vlasov problem for the dielectric/space craft 
configuration 0 / Interest. This type of solution can 
be obtained using the NASCAPlA computer code. 

For the purposes of this study a much 
simpler appraoch was chosen. Two distinct realizable 
configurations admitting analytic solutions were chosen. 
The variation from one solution to the other Is believed 
adequate to be Indicative of the type of change which 
might be expected were the details of the actual dis- 
tribution function to be used. The first configura- 
tion Is that of a semi-infinite surface, a one dimen- 
sional problem In which the pertinent dimension Is the 
distance from the surface. The equations describing 
this simple situation rre well-known. For an Isotropic 
Maxwellian at Infinity, the Incident distribution Is 
due to an unperturbed Maxwellian In the directions 
parallel to the surface, and a shifted cutoff Maxwellian 
In the direction perpendicular to the surface. This 
represents a shift In both energy and angle as a func- 
tion of surface voltage. For larger surface voltages, 
the Incident flux is more nearly mono-directlonil . The 
second configuration 1$ that which Is familiar from 
spherical probe theory. That Is, the surface ;s uni- 
form and spherical and the object size is small compared 
to the plasma Debye length. The incident flux In this 
case is that due to a shifted and cutoff Maxwellian In 
the total energy. The flux remains Isotropic as a func- 
tion of the surface voltage. 


The sample chosen for computation was a 2 
mil thick Kapton. The computations were performed with 
the following normal materials parameters: ralative 

dielectric constant, 3.45; density 1.23 g/^^; prompt 
conductivity coefficient, Kp * 1.16 x lO"*' mho/m/rad/ 
sec; dark field dependent conductivity given by the 
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dtu of AdIK). MconOtry yitid proporttorul 

to turfoct doM according to tho nodtl of lurk*, Wall 
and f rater Ickton.lA uting th« data cooptlatlon of 
lurkt.^' Tfirtt co«Tut4t1ont Mere parformed, one u>. 1 ng 
a 20 ktV aono-tnargattc noraally tncitent bain (labal* 
lad MONO); ona using tha surface flui taken from the 
ona*d1aont1ona1 aodal noted above (labelUd MAX), and 
one using the surface flux taken from spherical probe 
theory noted above (labelled SfT). For the latter 'mo 
cases, the distribution at infinity is the 11. b keV 
Isotropic Naxivelllan noted above. All threo coe 4 >uta- 
tions Mere performed with a unperturbed current density 
of 0.8 nA/cm, representing a severe charging environ- 
ment. 

As anticipated, the substantially higher 
secondary yield associated with the distribution of 
electrons (particularly the angular distribution) 
leads to substantially lower surface voltages for the 
space-like environment than for the mono-energetic 
charging situation. It should be noted that at sat- 
uration, the rather high field dependent conductivity 
of Kapton makes the bulk leakage current the dominant 
contribution to the equilibrium cuirent. Because of 
this fact, the equilibrium voltage for a monoenergetic 
source Is a very weak function of the beam energy at 
the given current densitv. Similar charging curves 
may be expected for smaller beam energies, and the 
choice of 20 keV My be viewed as representative. The 
relatively small difference between the two cases with 
space-like environments My be attributed to the dif- 
ferences In secondary yteld for the two. The secondary 
yield obtained from the SPT flux Is higher, leading to 
a lower surface voltage. Overall, the figure displays 
the well-known fact that space-Hke distributions give 
I'ise to substantially lower charging voltages than those 
obtained from mono-energetic charging simulations. In- 
clusion of an Ionic current Mkes this difference even 
larger. 

Shown In Figure 2 Is the location of the 
Man penetration of the charge In the dielectric (the 
charge centroid) as a function of tiM. The monoener- 
getic results are similar to those previously present- 
ed.^ The results for the space-Hke environment *iiow 
a very weak dependence on both model and tiM. The 
major feature of the figure Is the substantially deeper 
penetration associated with the space-Hke environment, 
showing the Importance of the higher energy electrons. 
Shown In Figure 3 are the fields Internal to the di- 
electric for the three charginn conditions at a time 
near saturation of the external potential. The bulk 
fields (to the right) show the variation (from model 
to model) expected due to the variation In surface 
voltage shown In Figure 1. The higher the surface volt- 
aM, the higher the bulk fields (distances In the di- 
electric beyond the penetration of the electrons). 

Nhat might not have been anticipated Is that the space- 
Hke environment leads to higher fields near the front 
surface of the dielectric, with the dependence on sur- 
face voltage reversed. The higher the surface voltage, 
the lower the front surface fields. The differences In 
this region are quite strong, being almost a factor of 
three between the monoenergetic and spherical probe 
theory slwlatlons. 

Consideration of the charging In the frame- 
work of the Meulenberg^ bilayer suggests that the be- 
havior of Figure 3 Is precisely what one would expect. 
Layer secondary yields are In the same order as the 
order of th.' front face fields. With similar amounts 
of burleo charge, the field In this region must vary 
with the surface positive charm. Shown in Figure 4 
Is a comparison of the internal enarge density at tne 
saM time as Figure 3. The Internal charge distribution 


for the two space-Hke slawlatlons Is virtually Identi- 
cal and not distinguished on the figure. On the dist- 
ance scale of the figure, the surface positive charge 
would appear as a aelta function and Is not Included 
In the display. The space-Hke environment shows the 
substantially deeper penetration suggested by Figure 2, 
and gives rise to the greater depth over which the 
electric field varies as shown In Figure 3. Over most 
• ,f ih,* rangi» of depth*., the dopo.tttd cnargt densities 
are very similar, being a few thousand couiombs/md (few 
tii«s lOlo electrons/cmJ). The space-Hke envlronrent 
dues not give rise to the strong charge enhancmxent 
very rear the front face seen for the monoenergetic 
simulation (which Is In part artificial, and reduced, 
by the front positive charge density). The toul nega- 
tive charge deposited In both cases Is very similar. 

Several further calculations were performed 
to examine the dependence of the space-Hke siMlatlon 
results on the problem paraMters. The transient con- 
duct1v1:y coefficient K was allowed to be *9 tiMS 
bigger and 10 tiras sMller than the nominal value, 
representing the range of typical values of K . The 
dark field dependent conductivity was al1owed''to he 10 
tiMS larger, representing a typical ‘‘fix* to the charg- 
ing problem. The secondary yield was allowed to be a 
factor of 2.7 times larger than nominal, representing 
both a possible "fix”, and also a large contribution 
to secondary yield due to Ionic envIronMnt. (For much 
larger values of , the Mterlal will charge positively) . 

All calculations were perfonaed using the spherical 
probe theory Incident envIronMnt discussed for the pre- 
vlcus results. Only one parameter Is varied from nom* 

Inal for each calculation. 

Shown In Figure S Is the surface charging 
voltage versus tiM for the nominal case (MOM) and the 
four variations noted above. Variation In It give rise 
to only minor variations In the final charging voltage 
as shown. Rearrangament of the charge In the deposition 
region has little effect on this variable. Of course, 
an Increase In Oq, or 6 gives rise to a substantial re- 
duction In surface voltage. That an Increase of Oq by 
a factor of 10 only gives rise to a voltaM reduction 
of a factor of two reflects the sharp field dependence 
of the dark conductivity In Kapton. The much slower 
charglno for the case of a larger secondary yield Is 
precisely as expected. However, even In this case. It 
Is the sharp field dependence of the conductivity which 
turns off the charging. 

In Figure 6 are shown the Internal electric 
fields at saturation for the saM paraMter variation. 

The variations In K give rise to variations In the 
field only In the deposition region. While It might 
seem surprising that the bulk fields are Identical when 
there are variations In surface voltage as shown In 
Figure 5, a little reflection should convince the reader 
that Figures S and 6 are both correct. Rearrangement of 
the charge In the deposition does effect the surface 
voltage but does not affect the bulk field. Within the 
deposition region the variations with IL are those ex- 
pected from previous work. 3 SMller ^values of IL 
give rise to larger values of the front face field. 

Increasing the bulk dark conductivity Is 
effective In decreasing both the bulk fields and fields 
In the deposition region. This Is expected. By contrast. 
Increasing the secondary yield leads to a decrease In the 
bulk fields and an Increase In the field near the front 
surface. This situation was discussed above, and Is 
expected within the frame work of Keulenberg bilayer.® 

One should also note the behavior of bulk 
fields for the large Oq and large f cases. Note that 
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tilt bulk fluid for Uit Itter 1i lar^tr dttplU tho 8. 

fact that tht aurfact voltas* It lowtr. This It a 
conittutnct of tht fact that a tuhttantlal auount of 
charft It dtpotittd at dtptht which art comparablt to 
tht taapit Uilcknttt (SO alcront on tht ftturt). Thut. 
for thin taa^ltt In a tpact>11kt tnv1roi«tnt, tht tur* 
fact vo 1 tt 9 t It not an accuratt diagnottic fo* tht bulk f. 
fitidt. Shown ;n Figurt 7 It tht chargt dtniity aitocl- 
attd with tht fitidt of Figurt 6. Thtit proflitt art 
contltUnt with tht abovt ditcuttlon. 


R.C. Naititon, R.J. Churchill, and E.J. Vadlow* 
tky. *NBaturtaitnta of Farticit Calttlen froa 
OltcharM Sittt In Ttflon Irradlattd by Nlth 
Intrty Citctron Ittat.* JUK Trtnt . [ic. kl . 
hS-», S141 (1879). 


J.C. Nantvicz and R.C. Adaao, "Spacocraft Charg* 
Ing Studitt of Voltagt Srtakdown Procttttt on 
Spactcraft Thtrwal Control Mlrrort,” rtf. S, 
Idta. p. 22S. 


OlWVSSlOh 

Tht Ntultnbtrg^ Idta that an Inttrnal dlt- 
chargt way bt Inltlatad In tht dtpotitlon laytr of , 
titctron chargtd dititctric hat bttn Invoktd. Ev1dtnct‘ 
hat bttn givtn which tnowt that Inttrnal ditchargtt do 
occur. S^ct obttrvatlont havt bttn nottd which tug* 
wtt that brtakdown conditions In tpact arc different 
froa thota In "ttandard" laboratory tlaulatlont. A con- 
ceptual aodtl wat prtttnttd which deaonttrattd that 
large Internal fitidt could develop with no differential 
potential. An tiptrlaental procedure for rtalliing 
theta conditlont wat nottd. Theta conditlont obvioutly 
favor an Internal brtakdown atchanita. New calculatlont 
were reported which deaonttrattd that tpace-likt tlaula- 
tlont give rite to conditlont which favor the Internal 
brtakdown atchanitfi; dtcrtattd bulk fitidt, inertattd 
dtpotitlon region fitidti timllcr chargt dcntUlct. 
Material t paraaetcr varlatlont were explored for the 
tpaca-llke tiaulatlon. The retulting changtt In Intern- 
al fitidt art tatlly underttandable In tcrat of the 
general theoretical aodtl. Inertattd conductivity It 
an effective deterrent to ditcharge. Enhanced tecondary 
eajttlon It not, and favort Internal ditcharge. 
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Figurt 1 . 



Surfic* PoUntUI v$. Tint for Thrt« CiMriIng CnvtronMntt 
(m* ttxt for conditions). 


Figurt 3. 


^turstlon Citctric Fitld Profllt Inttmol to Ololtctric 
for Thrtt Charging Cnvlroiwntt (sM Uit for cor.dltloni). 
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Figurt 2. 

Chargt Ctntrold vs. Tim for Thrtt Charging Envlronmnts 
(stt ttxt for conditions). 
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Figurt 4. 

Chargt Dtnslty vs. Dtpth Internal to Oltitctilc for Thrtt 
Charging Envlronmnts (stt text for conditions). 
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Figurt 

Surfict Charge Voltage vi. Tlw for Variations In Niterlal 
Faraaeteri (see test for conditions). 
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Figure 7. 

Charge Density vs. Oeoth Internal to Dielectric Near 
Saturation for Several Variations of Niterlal ParaMtars 
(see text for conditions). 


Figure 6. 

Electric Flclo v^. Depth Internal to Dielectric Near 
Saturation for Several Variations of Naterlal Paraneters 
(see text for conditions). 
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.1. f {rim..r 
»*ii>'U.N !•< [' irtiM'iit 
Sli'*V t)swi-^u, 1^1 '<> 

h. 1.. !.<•». rs 

liccr^ ‘.^sov ..itw' . 

1 ,u. B»ik JM' , 


!. I: I Hti! .V t 1011 . IIk hasit onu.it lonv j.uv.rn 'ii’. t i i 
liri.iiulj.iii .ir»‘ t tio fu- ii flow cnu.it inn 

L ' r • d . V i* HI .id ' J ■ Cl.‘ I i 


.ind tlic ri('i.it ii n of ..irriM v li.i : u 1 1 1 > 
.tiv I j!; i‘ 


wfuji tflO li.'ivo tlicir u-iUjI r.-amn^i T1 « .onNor 

VJt ion of energy is expressed by Lq I. .uid Jt ;• thorc- 
foio ru’ccssaiy. Equation 2 requires that the di^loc- 
tric rel.Tx.ition time l/c bo much vr..illir tinn t i-.o intir- 
■ivls Loniidtficd in the intcgr.it ion of l.q . 1. Tf. i s loiidi 
t ion holds for the examples worked lielow. 

If the electrical conductivity, o. r- not tin|Hra- 
lurc dependent, then the electric field, f. at inv point 
15 dttcmined only by the electrode geometry md I’o 
plied voltage. The problem then reduces t» th.it tro.itod 
in the clissic paper by i'u|ple, llartreo, iortri, .ind 
iyso'Ul) who o.il I. ul .itod the transient lumiit fir plino 
par ill el elect rode;. However, .i temper itiir dop. nd. i.t 
• Itv'rieil cor.diic t i\ 1 1 \ itself .i c.iuse oi '..n.it’iii r 
Mic electric field strength in order to '.|•l^fy Iq d. 
This < fleet .h.is been investigated by lVintli(») wh.i de- 
rived steady state solutions to the equ.it ions by .an.ilytic 
merhods for spciial functional forms of the temper.iture 
dependence of the electrical and thermal conductivities 
In the following wc develop transient solutions to the 
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t ion < for t!:ii.iil .n..l f 1 • r u ,i 1 iu.;.Uu t i > i ly fun; • 

f u»iii I'.iviMi fiy 

• / I (M 

u 

* n vxp f V^n (4) 

o 

rcspi c 1 1 vcl > , iktuTc <.’q, And A arc wonstant parameter!* 
of the dieirctru*. These functions are typical of the 
hchdvior of alkali halides. 


* • * I el LJ ee * rodes . For pi .me parallel c 1 cc - 

trodcs with t)iV : cooVXMiatc norm.il to the electrode sur- 
f.iccs, I and » bi .one with the n: e of o .iiid 1 


C M 
‘Tt’ 


I 


of , 
• • » 


K* ISO '-'.'IT' 


(3) 


' - j I. I \i) • .\ ' V r i 

fi ' 


toi 


whcri. ; IS the current lensity. Appl i c.it ion of a step 
function volt.ige it t » 0 i;ives the condition 


'o 


I 


d: 


\ for t > o 
o 


I’) 


where d is tiu' elcitiM.ie separ.it i »ii. AsMiitung electrodes 
of l.irjtc he.it c.ip.icity p. ives 

T » r for » o 

.iM»l : ' d 


win re T,^ is the .nh.en* t empi r 1 1 

Solution of the t'.irfial il i f f.rent i.i 1 eiput ion unJei , 
the coiiJiticiis «•, , .ind .'<;i >iel.ls the functions T(c,t) 

and Ff:,t) wiiich .ire both svirmetrie .I'lont the mid plane 
of the diel'Ctric slab he e.in therefore replace h he 

T « T for : - n 

” ' («.i) 

i>r/h; s n tor : = A/ 2 

lin.illv lit) • is net .i fiitu t loi of ’ , and 
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v.an hv calculated from 0. 

lor cox|Hitnt jons, the foll«>wing valuer were taken; 

Cj, » 2lM)t)K/m, Oq * cor.'*™'*, A » O.KbeV, C » 2x1 0*J/ra’r® , 
1^ < 80(lK, and d ■ Sxun. Fig- 1 shows the coiii|UJtrd cur- 
lent ilensity as a function of time for various applied 
voltages with the critical voltage apparently about 2.1kV. 
rig. 2 shows the distribution of temperature and electric 
field within the dielectric after steady conditions h-.-zc 
neer. attained. The field strength at the cent er U*2 . Smm) 
IS about half the average field strength, while near tlie 
electrodes it is about four times the average. This sug- 
gests the possibility that thermal effects might so dis- 
tort the field as to produce pure y electrical breakdown 
ni.ir the electrodes. The center temperature is about 
1 50C® above the ambient, and this is much higher than the 
temperature increase expected from the cor.stant field 
case 


2. Cyli ndrical Electrodes . For cylindrical geometVy 
with inner and outer electrode radii rj and r. respec- 
tively Eq. S is replaced by 

" 5V ■ F S ' T I? > ■ V' 

Application of a step function voltage at t > 0 gives 

/*■» E dr • for t > 0 (10) 

r, o 


if the outer electrode is m.ilQtained at ambient tempera- 
ture and the inner one allowed to take up the temperature 
of the adjacent portion of the dielectric, the boundary 
conditions become 


T ■ for r ■ r, 
o « 

3T/8r ■ 0 for r ■ r, 


( 11 ) 


Solution of the partial differential equation shows an 
instability in the current density similar to that shown 
in Fig. 1. The temperature distribution rises monoton- 
ically from ambient ne.nr the outer electrode to a maximum 
near the inner electrode. However, the electric field 
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1 st r ibut mil shows .in uit n st mu i!i v lopmoiit Jopoiul i n>j 
oil tlic -lOkti'J' r.ufii. ri>;i. ' mo J show t tu' iji.ti.^ 
fu'lJ .IS ,1 fiiiu'tioii ol' i.iili.'il Jl^tllU'c lor on .ip|>lui.l 
»olt.i»;o ut l.J.'lkV iiul I I I'v l roilo si| ,t i at ion of .<r'.fi In 
l>oth cas.s flic Jppliiil voltap.c is sli;’,l.tl\ in iso os • oi 
iritu'.il voltJRo. I i)j ^ illusti.itos t ho i h.ii'.p.v in t'lcld 
J I St r ihiit ion with time wh-ni t lio rit m r / r - At :ero 

t imo tho I'lcliJ non -uii i |oi i:;i t > is di t > ririnod hy ^oonotry 
only, after Ss thermal of foots pin, ^oom<'try nuiko tl»o 
field no.irly uniform, and after 10s thermal offoots Jomi- 
nato so that the field is larp,e near the niter i loot rode 
and small near the lontral ono. ( •(;. i illustr.itos tho 
different situation that prevails when f, • 6. In 
this ease, thermal offoets m.ike a oont r ihut ion to the 
field distortion, hut do not oveieomo tho •‘tronc, aoomot- 
rie contribution. 
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Kig. 1. Current tr«niiicntc foi pla:ic parallel 
electrodes at 5 nun spacinq and tor 
parameters quoted in the text. The 
various potential differences arc marked. 



Piq. 2. Variation of temperature and electric 

field between plane paraliel electrodes 
at 5 mm spacinq and 2.1 kV potential 
diffcrenco after attainment of steady 
c<^nditionn . 
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rx»|. Vatiatlon of lr“twocn 

I'yl Indrlcal c* Icct t'j n»n and 10 n«n 

radii) 10 a after .ippl iration of hV 

IHjtential dlfferanct'. 



Pig. 4 . Variation o( elactric- fl*ld hvtwran 

cylindrical alactrodat (1 nan and f> nan 
radii) 10!i aftar aj>rl of I.?*) kV 

I'otrntial dlfferoncr. 
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O'Dwyer, J.J., Space Charge Enhanced Emission into a Dielectric Solid, 
paper to be published in the Journal of Appl ied Physics . 
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by J. J. O' Dwyer 
S.U.N.Y. OswcRO, NY 13126 

A bstr act . Collision ionization may occur in a dielectric solid that 
is subjected to a large electric field, and the positive and negative 
species produced then constitute space charge within the dielectric. 
Emission from an electrode into the dielectric can be increased by 
substantial hcterocharge density adjacent to the electrode. It is the 
purpose of the current work to estimate both the dielectric conduction 
parameters and the electrode configurations that would be required for 
significantly enhanced emission. 

1. Introduction . We consider steady state current conduction by 
both electrons and holes in a dielectric in which a strong electric 
field is causing collision ionization. The conduction current density 
is givt'n by 


j « n e H ♦ p e Up ^ (1) 

where n and p are the electron and hole densities, and the corres- 
ponding mobilities. The electric field is and the electron charge 
is e. Tlie spatial variation of E is governed by Poisson's equation 

V . E « (p - n) e/£ ( 2 ) 

where r is the permittivity of the dielectric. The collision ionization 
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rate per unit time is assumed to be <(iven by 


3n 


I 



■ n a u E exp (-U/E) 
n 
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where a is an inverse length and H an electric field characteristic 
the dielectric. If diffusion and recombination are assumed to be 
negligible, the equations governing steady state electron and hole 
density are 

div (n u - n a u E exp (-H/E) ■ 0 
n n 

^ } (4) 

div (p u E) - n a u E exp (-H/E) ■ 0 
P P 

The basic problem is then the solution of the system of equations (1), 
(2), and (4) for given electrode geometry and electrode emission 
characterist ics . 

In a real dielectric, electrode emission probably begins from an 
asperity and the emitted charge begins to spread as it enters the 
dielectric. !t' the electrodes were, for example, parallel planes, then 
the emitted charge density would be large only at points within a 
reeion that spreads from the asperity towards the other electrode. In 
this case the ch.irge density pattern is not related to the electrode 
symmetry, and the problem is insoluble without additional assumptions. 
In order to propose a soluble problem that retains at least some ele- 
ments of the real situation, we write the equations in a form suitable 
for spherical geometry with emission taking place from the electrode 
of sm.iller radius. In these circumstances there are now two different 
reasons for enhancement of the electric field in the dielectric. In 
tlie first place the spherical geometry results in a higher electric 
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field near the electrode of smaller radius; secondly, the field near 
this electrode is further increased by local heterocharge. In the 
following calculations we separate these effects hy quoting a geometric 
field enhancement factor, and a total field enhancement factor. 


2 . Computatio ns for Spherical Geometry . Using the spherical fona 
for the divergence and writing th« results in the form of finite 
difference equations, we have from (1), (2), and (4j 

. . r — { j - j (1 ♦ HP) } ♦ Ar (5) 

cuc '^ n U r 

- P J 


exp (-K/E) - Ar (6) 

In addition the total current density obeys the equation 
Aj • - j- j Ar (71 

The boundary conditions used with these equations are: 

(1) At the anode j • j^, meaning that this electrode is blocking 
to hole emission 

(21 At the cathode j = j (E ) where the functional form used 

•^n •^n cath' 

in the following computations was that for Fowler-Nordheim 
emission (see e.g. O'Dwyer^. 

Numerical solution then consists of integrating the equations from the 
anode to the cathode and imposing the Fowler-Nordheim boundary condition. 
If the ratio of the electrode r.idii is g-eatly in excess of unity, then 
the first order expressions for the divergence in equations (5), (b), 

.ind l") lead to geometric factors on the field and current densities that 
are too low. he therefore drop the spherical geometry requirement and 
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specify the electrode system by its geometric enhancement factor and 
the i nt ere lect rode distance. 
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Thc following values were used for various piiramcters during the 
computations: 

e ■ 3 X 10 F/m (corresponding to * 3.4J 
H • 10® V/m 

7 8-1 

a Range 2 x 10 to 5 x 10 m 

Up Range 3 x 10 to 3 x lO’^^ m^/Vs 

U_ Does not enter computation if u ^ u 
n n p 

The electrode separation was taken at 90um unless otherwise noted, and 
the cathode work function at 2.5eV in the Fowler-Nordheim formula. 

Fig 1. shows the current -voltage characteristics for a geometric 
field enhancement factor of 22. In every case there is a current con- 
trolled instabilxcy with a certain critical voltage beyond which a 
negative resistance region corresponds to an unstable current -voltage 
relation. All curves are marked with the total field enhancement fac- 
tor that pertains at the critical voltage, and they show that the space 
charge enhances the field by an additional factor of about 2 to 3 over 
the geometric enhancement, the exact value depending on the collision 
ionization coefficient and the hole mobility. Fig 2. shows the same 
characteristics for an interelectrodc spacing of 290um and geometric 
field enhancement factor of 165; once again the space charge enhances 
the field by .in add it iona 1 f.ictor of al) 0 iit 2 to 4 . 

Fig .i. shows the effect caused by v.ariation of the geometric field 
enhancement factor while keeping the hole mobility, the electrode sep- 
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aration, and the collision ionization coefficient constant. The critical 
current density is relatively unchanged while the critical voltage is 
lower for larger field enhancement factors. The space charge enhances 
the field by an additional factor of about 2.5 to 5 over the geometric 
enhancement, the exact value depending on the electrode geometry. Fig 4. 
shows the effect of changir.g the electrode separation while maintaining 
the other parameters constant. Increased electrode separation results 
in larger space charge field enhancement, and larger critical voltage but 
does not change the critical current density. In both Figs 3 and 4 the 
current density is that near the cathode; this is in contrast to the pre- 
vious diagrams where the current density refers to the anode. 

3. Discus sion. Some general conclusions emerge from the computations, 
shewing the effect of varying certain parameters while keeping the others 
constant . 

(a) An increase in the collision ionization coefficient (all other 
parameters remaining constant) causes a marked reduction in the 
space-charge field enhancement factor. The critical current 
density is not substantially altered. 

(b) An increase in the hole mobility (all other parameters remaining 
constant) causes a marked increase in the critical current density, 
together with slight increases in both the critical voltage and 
the space-charge field enhancement factor. 

(c ) An increase in the geometric field enhancement factor (all other 
parameters remaining constant) caused a marked decrease in the 
critical voltage, while not substantially altering the critical 
curmit density. 
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(dj An im rcase in the electrudv sopurutiun (all uthcr para- 
meters remaining constant) causes n marked increase in the 
critical voltage but a decruasf in the mean field strength. 
The critical current is not substantially altered. 

In summary it appears that the most critical parameter for the 
possible existence of collision ionization space charge field enhance- 
ment is the hole mobility. For a cathode of radius lOum, the hole 
mobility should be of order 10*^^ m^/Vs or lower for the effect to 
occur at reasonable values of current density and voltage. Published 
results^’ ^ indicate values of this order for polyethylene, but it is 
not clear that a comparison can readily be made. The computations 
above consider the holes as the less mobile of the two species Involved 
in collision ionization, and the experimental results on charge carrier 
mobilities do not relate to the products of collision ionization. 

Given a siifficicntly low hole mobility, the existence of space charge 
enhancement at a given value of mean field strength is facilitated by 
the following factors: 

(.i) a large collision ionization coefficient 
( 1)1 a large geometric field enhancement factor 
(c) a large electrode separation 

In simple terms this means that avalanches of many generations from 
sharp asperities are most effective in promoting space charge enhanced 
emission - a fact that could have been seen without detailed computing. 
The major conclusion, therefore, remains as the estimate of the order 
of magnitude for the hole mobility (< lO”*’^ m^/Vs) required to produce 
*;ign I f leant space charge. 

i null,* a word sliould be said about the cxpcctoJ effect of tern- 
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peraturc. The paruiietors Must affected by a tcwpcratuii* change would 
probably be M und u^. both of which would be expected to increase with 
increasing tenper.iturc. This should cause higher values for both the 
critical current density and the critical field strength, with the 
exact details depending on the form of the temperature variation for 
the parameters in question. 
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Captions to Oidgraws 

Fig I. Cur rent -VO 1 1 age characteristics for a geometric field en- 
hancement factor of 22 and an electrode spacing of The total 

field enhancement factors at the critical point are ntarked on the 
curves. The collision Ionization coefficients refer to the curves 
above them. The hole mobilltes are denoted as fallows: 


-12 2 

- 3 X 10 ' mVVs 
P 

p ■ 3 X 10 *^m^/Vs 
P 

-lit 2 

p - 3 X 10 m^/Vs 
P 


Fig 2. Current -VO* tage character i st ics for a geometric field en- 
hancement factor of 165 and an electrode spacing of 290pm, The total 
field enhancement factors at the critical point are marked on the 
curves. The collision ionization coefficients refer to the curves 
above them. The hole mobilltes are denoted as follows: 


p ■ 3 X 10 *^m^/Vs 

P 

- 3k I0‘*^m^/Vs 

P 

-\k 2 

p “3x10 m /Vs 
P 

Fig 3« Current-voltage characteristics for a range of geometric 
field enhancement factors leading to the total field enhancement 
factors marked on the diagram. In all cases the electrode separa- 
tion was 90pm, the collision ionization coefficient lO^m'^, and the 
mobility 3 x lO'l^/Vs. 


Fig k. Current-vol tage characteristics for two different electrode 
separations and the sar>e geometric field enhancement factor of 22. 
The total field enhancement factors at the critical point are marked 
on the. curves. In both cases the collision ionization coefficient 


was 


I0®m"^, and the hole mobility 3 x I0“'^m^/Vs. 


148 



(A/m^) 


0 



-7 

4 6 8 10 

V(k V) 


149 



150 



6 8 10 
V (kV) 

151 








